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MEETING OF 1952 JANUARY 11 
Professor H. Dingle, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


John Anthony Bancroft-Wilson, M.A., 18 Severn Road, Weston-super- Mare, 
Somerset (proposed by D. W. Dewhirst); 

William George Barrett, Jun., 13 Heversham Road, Bexleyheath, Kent 
(proposed by A. C. Clarke); 

John W. Boyle, 9 Pearson Grove, Caulfield, Victoria, Australia (proposed by 
N. R. L. Blackmore); 

*Geoffrey Martin Brown, Department of Physics, University College of 
Swansea, Singleton Park, Swansea (proposed by F. J. M. Stratton); 
Phyllis Betty Comrie, 131 Maze Hill, Blackheath, London, S.E.3 (proposed 

by W. M. H. Greaves); 
*Phillip John David Gething, Ph.D., M.Sc., 6 Highmore Road, Blackheath, 
London, S.E.3 (proposed by G. J. Whitrow); 
Murdoch MacDonald, K.C.M.G., C.B.(Mil.), M.I.C.E., 72 Victoria 
Street, London, S.W.1. (proposed by F. J. Hargreaves); 
Edith Miller, D.Phil., University Observatory, St. Andrews, Scotland 
(proposed by R. O. Redman); 
Colin Andrew Murray, M.A., Royal Observatory, Greenwich, London, 
S.E.10 (proposed by R. d’E. Atkinson); 
C. Joy A. Penny, B.Sc., Norman House, St. George’s Road, Bickley, Kent 
(proposed by R. d’E. Atkinson); 
*Vincent Cartledge Reddish, B.Sc., 40 St. Helens Road, Leigh, Lancashire 
(proposed by A. Hunter); 
*Dennis William Sciama, M.A., Cavendish Laboratory, Cambridge (proposed 
by ‘I’. Gold); 
*Bryan Summers, 365 New King’s Road, London, S.W.6 (proposed by 
J. McG. Bruckshaw); 
Frank Cyril Drake ‘Thomson, 13 Akenside ‘Terrace, Jesmond, Newcastle-on- 
‘Tyne, 2 (proposed by H. Thomson); and 
John Victor ‘Thomson, Mirastelle, Woodland Way, Kingswood, Surrey 
(proposed by F. J. Hargreaves). 
I 
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The election by the Council of the following Junior Members was duly 
confirmed :— 


Kenneth Neilson Dodd, B.Sc., M.Sc., 36 Prospect Vale, Wallasey, Cheshire 
(proposed by W. H. McCrea); 

Charles Beresford Rayner, 3 Warminster Road, South Norwood, London, 
S.E.15 (proposed by G. J. Whitrow); and 

George Alan Wilkins, B.Sc., Ph.D., Royal Greenwich Observatory, Herst- 
monceux Castle, Sussex (proposed by D. H. Sadler). 


Sixty-eight presents were announced as having been received since the last 
meeting, including :— 
S. A. Mitchell, Eclipses of the Sun, 5th edition (presented by the author); and 


Third Cape Catalogue for the equinox of 1925.0 (presented by the Royal 
Observatory, Cape of Good Hope). 


The President announced that the Council had awarded the Gold Medal of 
the Society to Dr John Jackson, lately H.M. Astronomer at the Cape, for his work 
on stellar parallaxes and his contributions to the general problems of star positions 
and proper motions. 


* Transferred from Junior Membership. 





STAR-STREAMING FROM THE PROPER MOTIONS OF THE 
BOSS GENERAL CATALOGUE 


T. R. Tannahill 


(Communicated by the Director, University Observatory, Glasgow) 


(Received 1951 September 14) 


Summary 


1. The proper motions of the Boss General Catalogue have been analysed, 
according to the two-streams hypothesis, by Eddington’s method. A total of 
22 374 Stars of spectral types A5 to M have been used in the discussion. 

2. The galactic coordinates of the vertex are (343°°1, —1 °7). 

3. The equatorial coordinates of the solar apex are (273°°2, -+35°°4). 

4. There is a marked decrease in the solar speed and in the relative speed 
of the drifts with advancing spectral type, and a similar decrease with 
increasing galactic latitude. 





1. ‘(he proper motions used in the present work are those of the General 
Catalogue of 33 342 Stars by Benjamin Boss(1). Preliminary work—the forma- 
tion of a card-catalogue—was begun in 1939; various diversions, however, have 
combined to prevent the completion of the investigation until now. It is my 
pleasant duty to record my thanks to Dr R. A. Robb for help with the card- 
catalogue, and to Professor W. M. Smart for his encouragement throughout. 

2. The object of the investigation is the determination of the constants of the 
star-streams and the direction of the solar motion, and the variations of these 
quantities with the spectral type and galactic latitude of the stars concerned. In 
previous investigations on the Cape photographic proper motions (2, 3), it was 
found that the stars of spectral type A form an anomalous group as regards the 
phenomenon of star-streaming, and 8751 of these have been omitted, with 
1626 stars of types O and B, from the present discussion. Of the 33 342 stars 
in the catalogue as published, 476 have no spectral type assigned, and are also 
omitted, as are three stars designated as having peculiar spectra. Finally, 63 stars 
of the Taurus cluster (4) and 20 members and 28 doubtful members of the 
Ursa Major cluster (§)—not already omitted as being of early spectral type—were 
rejected. The Pleiades stars, being of early spectral type, are not concerned in the 
discussion; nor are the stars of the doubtful Scorpio—Centaurus group. ‘There 
remain for discussion 22 374 stars * of spectral types from A5 onwards. ‘These 
have been divided into three spectral groups designated A5-F5, F8—G5 and 
Ko-M, the last group containing also the relatively few stars with spectral types 
later than M. 

The Boss catalogue contains, besides the proper motions, estimates of the 
probable errors of the determinations in right ascension and declination. In the 
present work a rough system of weighting was adopted. A star is assigned 
half-weight if the probable errors of the centennial proper motions in both right 


* One star has been mislaid during the investigation. 
1* 
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ascension and declination exceed 0”-075, or if either of these quantities exceeds 
0”-085. ‘This system is by no means perfect, but it is rapidly and easily applied. 
The process reduced the total weight of the stars used in the investigation to 
17953. It accounts for the presence of “ half-stars” in Table II. 

3. Systematic corrections to the Boss proper motions have been calculated by 
Smart (6). In applying these to the data, a preliminary division of the sky into 
small areas has been made. Each area is bounded by two parallels of declination 
20° apart, and by two meridians 10° apart in right ascension, giving 324 areas in all. 
Equatorial areas (declination — 10° to + 10°) are indicated by the letter E, areas of 
north declination by A, B, C and D, and areas of south declination by «, f, y and 6, 
as shown in Fig. 1. ‘The right ascension of the centre of an area is indicated by one 
of the numbers 1-36. ‘Thus area 15C is bounded by the parallels +50° and 

+70°, and by the meridians 145° and 155° of right ascension. 

For the centre of each area the corrections in right ascension and declination 
have been obtained by interpolation from the tables given by Smart, and applied 
to the components of proper motion of each star within the area. ‘Taking into 
account the probable errors of the proper motions, and the uncertainties of the 
corrections themselves, this procedure is considered to be sufficiently exact. 

4. For the final analysis the sky was divided into 37 regions arranged, as far 
as possible, to lie in well-defined galactic zones. ‘These regions are built up by 
combining the small areas described above. ‘This division (of one hemisphere 
between 5° and 185° of right ascension) is shown in Fig. 1. Diametrically opposite 
regions of the sky are combined as explained in Section 5. During the analysis 
according to spectral type, it became apparent that certain of the regions near the 
pole of the celestial equator were too small in area and contained too few stars to 
permit of satisfactory analysis. Region 9 was accordingly combined with the 
region antipodal to 15, the composite region being designated 9-15; regions 27 
and 36 were also combined. ‘The original numbering of the regions has been 
retained throughout in order to prevent confusion. 

Details of the limits and geometric centres of the regions, with the total number 
of stars in each, are given in Table I.* The arrangement described gives, in effect, 
five galactic zones and two polar calottes. Zone I consists of regions 1-6, the 
mean galactic latitude being —0°-2. Two zones consisting of regions 7-12 with 
the region antipodal to 15 (mean galactic latitude +19°:5) and regions 13, 14, 
16-19 (mean galactic latitude —20°-7) form Zone II of the subsequent analysis. 
Regions 20-24 (mean galactic latitude +42°-5) and 25-30 and 36 (mean galactic 
latitude — 43°-4) similarly form Zone III. Zone IV consists of the regions 31-33 
{mean galactic latitude +64°-8) and 34, 35, 37 (mean galactic latitude —70°-1). 
Each region—and each zone as defined—is combined with the part of the sky 
diametrically opposite. For.most practical purposes, therefore, the galactic 
latitudes of the zones adopted may be considered as 0°, 20°, 43° and 67°. 

5. From the components of proper motion in right ascension and declination, 
corrected as described in Section 3, the total annual proper motion and position 
angle, to the nearest 0"-001 and 1°-0 respectively, were obtained by a graphical 
device. ‘The position angles only are used in Eddington’s method (7) of analysis 
into star-streams, the method adopted in this paper. ‘The reference meridian 
for position angle, for an individual star, is that through the star itself. No 
appreciable error is committed if, in dealing with a small area of the sky, the 


* Tables I to VI appear at the end of the paper. 
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reference meridian for all stars in the area is assumed to be that of the centre of 
the area. If, however, the regions are large (as in the present case), this error is only 
negligible where the meridians through the stars and that through the centre of the 
region are substantially parallel, i.e. near the celestial equator. For non-equatorial 
regions, therefore, the following procedure has been adopted to allow for the 
convergency of the meridians. Each small area—as described in Section 3—has 
been treated separately. Within each the position angles have been corrected 
by an angle $=(c—«,)sin8, 
(x, 5) being the coordinates of the centre of the region, and «, the right ascension 
of the centre of the small area. For the most part these corrections are very small. 
In regions of high declination, however, they may be greater than 50°, but such 
circumstances apply to very few stars. Regions in which the corrections for 
convergency have been applied are marked with an asterisk in Table I. 








Fic. 1.—-Division of the celestial sphere. The faint broken lines show the divisior into small 
areas used in applying catalogue corrections and the corrections for convergency of the meridians. 
The heavy broken line is the galactic equator, K the south galactic pole, O is the Drift I apex, 
@ the Drift II antapex, © the vertex of star-streaming (opposite to centre of galaxy), @ the 
antapex of the solar motion. The central meridian is at right ascension 95°. 

‘The regions shown in Fig. 1 have, as already stated, been combined with the 
diametrically opposite regions of the sky. If a star has position angle 6 referred to 
a certain meridian, then referred to the anti-meridian, its position angle becomes, 
for our purpose, 300° —@. 

6. ‘The position angles in each region, treated as described, have been grouped 
in 10° sectors, and the resulting distribution smoothed by taking running means of 
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three successive sectors. ‘The frequency curves so obtained were then analysed 
by Eddington’s “trial-and-error” method (7), fitting a curve of the form 
p =0°0492{N, G(AV,, 0 —8,) + NoG(AV 2, 0 —8,)}, 

where N,, Nz are the numbers of stars in the region belonging to Drift | and 
Drift II respectively, AV, and AV, are the projections of the space-velocities of the 
drifts on the tangent-plane at the centre of the region, #, and @, the position angles 
of the directions towards the respective apices, p the number of stars moving in 
direction 6. ‘The function ¢ has been tabulated by Smart (8). In all cases 
N, +N, is the total weight of stars in the region. Considerable difficulty was 
experienced in determining the constants in certain regions of the sky—parti- 
cularly those of Drift I] in the regions surrounding the Drift I apex, where the 
indeterminacy of @, and the low value of AV, make the determination of hV, and 0, 
very uncertain. Generally, however, the separation of the drifts was not difficult, 
the regions being well populated ; the best determinations required an estimate of 
#, to within 24°, of 0, to within 5°, and of AV, and AV, to within 0-05, in order to 
obtain a satisfactory fit. 

The results of the analysis, for the three spectral groups and for the total data, 
are shown in ‘Table II. 

7. From the analysis two equations of condition for each drift are obtained 
for each region: 
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X sina — Y cosz =—hVsin@, 
X cosasinéd + Ysinasind—Zcos 6= —hV cos. 


X, Y and Z are the linear components of the space-velocity of a drift relative to 
the Sun, referred to the customary equatorial axes, and «, 5 the coordinates of the 


centre of the region as given in Table I. In forming the normal equations equal 
weight was given to all regions. ‘The solution of the normal equations gives 
X, Y and Z, from which the coordinates A, D of the drift apices may be found from 
the equations 

tanA=Y/X, tan D=2Z2(X*+ Y?) 12, 
and the space-velocities hW of the drifts from 

hW =(X? + Y? + 2)". 

The results are shown in Table III, with the probable errors of the determinations. 
Separate solutions have been obtained for each galactic zone. ‘The total weights 
of stars in each drift are given by n, and ng. 

The results largely confirm those of previous work, as may be seen from 
Table A, in which the present results for the total data are compared with the 
corresponding results obtained by Eddington (7) from a discussion of the proper 
motions of the PGC, and with those obtained from the Cape photographic 
proper motions (3) ox stars within the same range of spectral type as is used in the 
present work, situated in the zone of south declination 40° to 52°. 


TABLE A 


PGC GC 


go°'8 
—14°°6 
1°516 


287°'8 
—64°°1 
o'855 


93°°4+0°'9 
— 9°9t0°9 
1*500 + 0'022 


274°°4+3°°6 
—7o°*1+1°°2 
0°797 + 0°021 


Cape 
go°-8+ 1°°7 
= Peers 
1°545 + 0°041 
304°°9 + 8°-4 
—76°*2+2°°0 
0°668 + 0029 
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A noteworthy feature of Table III is the decrease in the value of AW, with 
advancing spectral type. This variation with spectral type has already been 
noted in the previous work on the Cape data (3), but the probable errors of the 
determination were much higher in that case. The present results give a more 
“efinite indication that the effect is not accidental. No systematic trend is dis- 
cernible in hW,. There is a possibility of a slight increase in D, with advancing 
spectral type, but in view of the probable errors involved no definite conclusion 
may be drawn. 

The value of hW, shows also a definite decrease with increasing galactic 
latitude. A similar effect has been found by Jones (g) in discussing the PGC 
stars: namely, a decrease from 1-727 at the galactic equator to 1-415 at latitude 
43°. The present result—a decrease from 1°703 at the equator to 1-392 at 
latitude 43° and 1-364 at latitude 67°—amply confirms the conclusion reached 
by Jones. 

No other significant variations with galactic latitude are apparent. 

8. With the known values of the drift constants it is possible to calculate the 
theoretical values of AV and 6 for each drift. This has been done for each grouping 
of the data, and the results are shown in the columns headed AV, and @, in 
Tables IVa-IVe. ‘The residuals hV,—hV and 0,—@ are also shown. ‘These 
are small, generally, but high residuals occur—as expected—in regions near the 
drift apices, and—-unexpectedly—in a few regions where there seems to be no 
excuse for their occurrence, e.g. in the value of AV,.—AV, of region 28, Table IVB. 

A notable feature of Tables 1Va—-IVp is the preponderance of negative residuals 
in the upper half of the AV,—AV, column, and of positive in the lower. This is 
a consequence of the decrease in AV, with increasing galactic latitude, and the 
fact that this effect shows itself even in the residuals of the spectral groups is 
further confirmation of its reality. The effect does not, for obvious reasons, 
show in ‘Table IVe. 

A somewhat better indication of the convergence of the great circles defined 
by the observed position angles of the drifts may be obtained by calculating the 
distance £8 of the apex from each great circle. The appropriate formula is 


sin 8 =sin Asin (@,—@), 


where A is the distance of the apex from the region concerned. ‘The values of 
A and f are given in Tables 1Va-IVe. As regards £,, only seven values exceed 
10°, and all of these occur in the separate spectral groups, the highest value being 
19°-0. Larger variations occur in f,, a consequence of the lower AW value for 
Drift II and the greater difficulty of determining its direction in any particular 
region. 

g. From the values of X, Y and Z given in Table III we calculate the linear 
components €, », ¢ of the motion of Drift I relative to Drift II, and hence the 
right ascension «, and declination 5, of the vertex of star-streaming and the 
relative space-velocity {2 of the drifts. ‘The results are shown in Table V. There 
is a marked decrease in {2 with advancing spectral type—which might be expected 
as a consequence of the corresponding variation in AV,. The Cape proper 
motions (3) also show this variation quite clearly. A decrease of 2 with increasing 
galactic latitude is also present, but is somewhat masked by the relatively poor 
determination of this quantity in Zone I. The effect is completely obscured in 
Jones’s investigation of the PGC stars (g), owing doubtless to the higher probable 
errors involved. 
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Table B gives a comparison between the results obtained for the vertex of 
Star-streaming in the present paper and those obtained in some previous 
investigations. Column (a) contains Eddington’s results from the PGC stars (7), 
column (b) those from the Cape photographic proper motions (3), and column (c) 
those of the present paper. An analysis of the GC stars between magnitudes 
6-0 and 7°5 has recently been made, on the ellipsoidal hypothesis, by Delhaye (10). 
The grouping used by Delhaye is rather different from that in the present paper, 
the material being divided into three spectral groups B8-A5, Fo-F9, Ko-Ka, 
and into zones of high and low galactic latitude. The figures given in column (d) 
represent an average value, for the whole sky, of the Fo-F9 and Ko—K2 groups, 
comprising 7187 stars of the magnitude range mentioned. ‘The corrections 
applied to the proper motions in Delhaye’s investigation were those calculated 
by Blaauw and Delhaye (11), and represent a compromise between the GC 
and the FK3 systems. It is evident that, in spite of the marked difference in 
the data and in its treatment in Delhaye’s work and in the present paper, the two 
investigations agree with the earlier work in attributing to the vertex of star- 
streaming a direction significantly different from that towards the centre of the 
galaxy, namely 325° galactic longitude. For the brighter stars at least, therefore, 
the outstanding discrepancy in the theory of star-streaming remains. 


TABLE B 
(a) (6) (c) (d) 


PGC Cape GC GC 
ly 347° 342° 343° 337° 
'e —0°'5 2 —2 “oe: 
Q 1°87 1°75 1°82 


It is to be remarked that the values of €, y and { in Table V refer to the vertex 
opposite to the centre of the galaxy. 

10. From the analysis of Table II the ratio n,/n, of the populations of the 
drifts may be obtained. If hU, is the solar motion relative to the stars concerned, 
and (Ap, Dy) is the apex of the motion on the celestial sphere, we may calculate 
hU, and (Ap, Dy) from the equations 


(n, +n.)hU, cos Ay cos Dy +n, hW, cos A, cos D, +n,hW, cos A, cos D, =0, 
(n, +n,)hU, sin Ay cos D, +n, hW, sin A, cos D, +n, hW, sin A, cos D, =0, 
(n, + my)hU, sin Dy +n,hW, sin D, +n ,hW, sin D, =0, 


the quantities 4,, D,, AW, A,, D,, AW, and the ratio n,/n, being known from the 
previous analysis. 

The results of the solution are shown in Table VI, with the values of n,/n, 
for each grouping of the data. Consequent on the corresponding variations in 
AW, (Table III), there is a weil-marked decrease in hU, with advancing spectral 
type and with increasing galactic latitude. ‘The decrease in AU, with advancing 
spectral type is noticeable also in the results from the Cape data (3). 

Table C gives a comparison of the results for the solar motion as obtained 
(a) by Eddington (7) from the PGC stars, (6) by Smart and Tannahill (3) from the 
Cape photographic proper motions, (c) in the present investigation, and (d) by 
Delhaye (10) from the GC stars. In column (d), as in the corresponding column 
of Table B, the figures given represent a weighted mean of those given by Delhaye 
for the spectral groups Fo—-F9 and Ko-Ka. 
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TABLE C 
(a) (6) (c) (d) 
PGC Cape GC GC 
ly 29 20° 30° a5”. 
bo +26 +23° +21‘ +20° 
AU, o'9I 084 0°83 
n/N 1°48 I‘Io 1°19 


It is to be noted that Eddington’s investigation of the PGC stars included the 
type A stars; the results in column (a) are therefore not quite comparable with the 
other results quoted. With regard to the value of m,/m., numerous investigations 
since Eddington’s work have agreed in bringing the value nearer unity, parti- 
cularly for the fainter stars. ‘The Cape stars (2), however, when the early types 
are included, give the value 1-41, in good agreement with Eddington. 


University Observatory, 
Glasgow, W.2: 
1951 September 13. 
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TABLE | 
Limits Limits 
of of 5 
R.A. Dec. 


° ° ° ° ° ’ ° 
95 > 115 +10 + —30 — 9 24 190°6 
105 145 —30 —50 —3916 22574 
145 185 ——). = —5831 257°2 

;: . £8 +70 +5831 1003 
55 75 +30 +4603 124°6 
75 95 +28 o1 148°5 


1941 215°5 
185 —39 16 249°8 

75 +75 54 ge"o 

75 . +62 00 120°4 
95 115 +28 o1 156°7 
135 + 924 183°3 

85 —46 03 221°3 
—58 31 241°! 

5 : +39 16 98°3 
35 ‘ +28 o1 123°9 
55 + 924 15371 
75 + ‘ — 924 181°5 


135 + 9 24 194°7 
155 —19 41 245°7 
+7554 101°5 

115 - +58 31 123°9 
115 +39 16 1513 
45 fl —39 16 208-9 
25 —§831 237°5 
5 —75 54 255°8 

5 § t +19 41 102°6 
25 © 00 140°9 
45 —19 41 181°0 


155 +19 41 200°1 
185 0 00 232°2 

165 +46 03 119°2 
25 45 —50 —28 o1 186°0 
5 25 =o So —46 03 259°0 

5 25 +10 —30 — 9 24 105°2 


* Position angles corrected for convergency of meridians. 
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TABLE II 
Analysis of the Drift Curves—Drift I 
F8-Gs5 Ko-M 

AV, 4, kV, 4, N, 
270 = 130 
310-1454 
2824 116 
120 «112 
140 87 
160 =:165 


1'05 
270 «=. 285 
1°6 
1°3 
I'l 160 


a= 00 
uw 


o8 265 

1'7 

19 95 

270 = «185 
180 


oa 
wn 


275 
270 

95 
180 
215 
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‘TABLE ITI (cont.) 
Analysis of the Drift Curves—Drift II 
F8-G5 Ko-M 
AV, 9% hV, 9, Nz 
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06 = 160 06 160 
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Taste III 
The Drift Constants 
Drift I 
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TABLE III (cont.) 
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Tase [VB 
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TABLE IVE 
Galactic Zones 
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TABLE V 
The Vertex of Star-Streaming 
As-F5 F8-G5 Ko—-M A5-M 
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TABLE V (cont.) 
The Vertex of Star-Streaming 
Zone I Zone II Zone III Zone IV 
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A TWO-PLATE CATADIOPTRIC SYSTEM FOR ASTRONOMICAL 
PHOTOGRAPHY 


P. A. Wayman 
(Received 1931 September 28) 


Summary 


It is shown that a two-plate aplanatic optical system with a single spherical 
mirror can be designed to give good images over a considerable field and with a 
smaller tube-length than is used in the classical Schmidt camera. Alter- 
natively, an anastigmatic design with a non-spherical mirror is possible. 





Wright (1) and Vaisala (2) suggested a method for overcoming the disadvantage, 
in the classical Schmidt camera, of great overall tube-length. The condition of 
zero Seidel astigmatism is abandoned and the corrector-plate at the centre of 
curvature of the mirror is replaced by a corrector-plate nearer to the focal surface. 
The curvature of the focal surface remains unchanged but the mirror has to be 
figured with a non-spherical profile. Primary astigmatism is then the dominant 
aberration and is large enough to prohibit application of this principle to systems 
of appreciable power. With a total tube-length of about six feet, the circle of 
least confusion 1}° off-axis in a 12-in. {/6 system would be 50 microns in diameter. 
If the total tube-length is increased to eight feet the maximum image-diameter is 
reduced to the tolerable value of 33 microns. For systems witha focal length much 
less than six feet, the disadvantage of long tube-length in the classical Schmidt 
camera is not serious. ‘Thus it would appear that the useful applications of the 
Wright—Viaisala system are definitely limited. Furthermore, there is an additional 
disadvantage in that the main mirror has the figure of an oblate spheroid, presenting 
difficulties with regard to optical testing. 

With the use of two corrector-plates there will be more freedom in the optical 
design. The “ plate-diagram” of such a system is shown in Fig. 1. (See Burch 
(3) for the original discussion of plate-diagram analysis.) 


= y, % 
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We take as unit of “strength” for the corrector-plate that of the anastig- 
matising plate for a spherical mirror of radius R, which we suppose to be the 
paraxial radius of curvature of the main mirror. , and y, are the strengths of two 
corrector-plates, at distances d, and d, from the point which represents the 
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paraxial centre of curvature of the mirror. « is the strength of the figuring on the 
mirror in the same units. (A value «=1 would apply to a paraboloidal mirror 
and « =0 to a spherical mirror.) 
For zero spherical aberration, ‘in Seidel approximation, we must have 

at+y, + V¥e—1=0. (1) 
For zero primary coma 

714) + Yodg+aR=0. (2) 
‘The astigmatism will then be proportional to the quantity 

A =y,d,? + yods? + aR. (3) 
The Wright-Vaisala system with a total tube length 4 times the focal length is 
given by 
¥2=0 é d, =4R. 

‘The eqns. (1) and (2) require that 


c= } and V1 = :. 


Substituting these values in (3), we have A = — $R?. 

Suppose now that we consider a two-plate system with a spherical mirror also 
having an overall length 4 times the focal length. ‘The position of the first plate is 
given by d, = }R, and we will choose to put the second plate immediately in front 
of the focal surface and so will write d,=}R. Also, « =o for a spherical mirror. 

From (1) and (2) we have, inserting these values, 


Nt+y2=1, 
and hy, + 42 =0, 
giving m=3 j.and y= 
From (3) A= —R?/6. 


‘Thus with a spherical mirror and two corrector-plates we have less astigmatism 
(by a factor of two) and no greater overall length than in the one-plate Wright— 
Vaisala system. With a total field of 3°-75 and an aperture of 15 in. a system with 
a focal length of six feet could be built in which the diameter of the largest astig- 
matic images would be 33 microns. With an overall length of eight feet and with 
a spherical mirror such a system would constitute a convenient and useful 
instrument for objective-prism photography. 

Alternatively, we can investigate the solutions of the eqns. (1), (2) and (3) with 
the astigmatism A=o. If, as before, we take d,=4R, d,=}R we easily obtain 
vm1=45, Y2= — 4, a =0°5, 

Thus an anastigmatic design is possible in which the mirror has a figuring-layer 
half that of a normal paraboloid, i.e. that of a prolate spheroid. ‘This does not 
have the same disadvantage with regard to testing as an oblate spheroid. 

Full consideration of the image-quality of these systems will be treated in a 
further paper, but a brief examination of the chromatic aberrations to be expected 
seems appropriate here. 

In the case already discussed of the 15-in. f 4-8 system with a spherical mirror, 
the on-axis chromatism need be little more than that of the corresponding classical 
Schmidt camera, since the “strengths”’ of the two plates partly compensate one 
another. Off the axis there would be a slight “shearing” of the plates which 

2* 
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would lead to a small asymmetrical colour-error. Now the maximum plate-slope of 
a corrector-plate determines the greatest amount of lateral chromatic aberration 
which that plate introduces. It can easily be shown that the maximum plate-slope 
occurring on the first and stronger plate of the system here considered is less than 
that on the plate of an f/3 classical Schmidt camera with the same focal length, in 
which the colour-error is known to be satisfactorily small. ‘I‘hus, even ignoring 
the fact that in the case of the off-axis pencils the chromatic aberrations of the two 
plates still compensate one another to a certain extent, it can be seen that the total 
chromatic aberrations need not be excessive. 

In the case of the anastigmatic design with a non-spherical mirror, it is, in fact, 
the chromatic aberrations which will set a lower limit to the focal ratio which can be 
used. A detailed investigation is necessary in order to show exactly what this 
limit is likely to be, but the considerations of the preceding paragraph, together 
with the fact that in this case the strengths of the plates, though larger, are more 
nearly equal, indicate that it will not be unduly large. Furthermore, the system is 
intended to be relatively inexpensive and if larger optical elements are contemplated 
other designs would probably be preferable. 


Yerkes Observatory, 
Williams Bay, Wisconsin : 
1951 September 15. 
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THE VELOCITY DISTRIBUTION OF SPORADIC METEORS. II 
Mary Almond, . G. Davies and A. C. B. Lovell 


(Communicated by A. C. B. Lovell) 


(Received 1951 September 28) 


Summary 


The work on the velocity distribution of sporadic meteors described in 
Part I has been extended to fainter meteors. ‘The results are similar and 
give no evidence for a significant hyperbolic velocity component. A fuller 
discussion of the magnitude range is given and it is concluded that the 
majority of meteors in these new experiments were probably in the magnitude 
range 6:0 to 7°5. 





i. Introduction.—Part | of these papers* gave an account of the application 
of the radio echo technique to the study of the velocity distribution of sporadic 
meteors. In this paper the extension of the work to fainter meteors is described. 
An alternative method of estimating the magnitude range leads to more definite 
information regarding the meteor magnitudes covered by the series of experiments. 

2. Apparatus. —The apparatus was basically similar to the 8-13 metre apparatus 


described in PartI. ‘The aerial system was fixed to produce a beam at an elevation 
of 10° with its main axis in an easterly direction (azimuth 90°) as before, but its 
gain was increased by the addition of more elements. 

The parameters of the apparatus were as follows :— 


Wave-length A=8-13 metres 
Peak transmitter power P=240kW 
Aerial beam width Azimuth = +12 
to half power i Elevation = + 5 
Aerial power gain over 
isotropic source G = 36°6 
Receiver noise level €) = 1-3 x 10°! watts 


The additional sensitivity over the apparatus described in Part I was achieved 
by increases in transmitter power P and aerial gain G. ‘The apparatus and 
technique for velocity measurement remained the same. 

With this equipment two series of experiments were carried out: (i) during 
the autumn mornings of 1950 (apex experiment), (ii) during the spring evenings 
of 1951 (antapex experiment). 

3. The 1950 apex experiment.._Due to the specular reflection properties of 
meteor trailst an infinitely narrow beam directed at azimuth go° at 06 would 
record only those meteors lying in a great circle plane through the apex of the 
Earth’s way, cutting the horizon in the north and south. The meteors recorded 

* Almond, Mary, Davies, J. G. and Lovell, A. C. B., M.N., 211, 585, 1951. 


t Hey, J.S. and Stewart, G.S., Proc. Phys. Soc., §9, 858, 1947. Lovell, A. C. B., Banwell, C. J. 
and Clegg. J. A., M.N., 107, 164, 1947. 
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by a beam of finite width, such as that used in the present experiments, will be 
those whose radiants lie between two small circle planes enclosing the apex. 
It is desirable to carry out such an experiment when the apex is in the neighbour- 
hood of its greatest altitude so that the sporadic meteor rates are a maximum. 
The first experiment was therefore made between the hours of approximately 
05" and 07" in the autumn mornings. ‘The details of the experiment which was 
carried out between 1950 November 8 and 1950 December 14 are given in 
‘Table I. In 10-7 hours of observation 335 velocity measurements were made. 
The times and details of these measurements are given in Table IJ and the 
number distribution of velocities in Fig. 1. 


TABLE | 
Details of apex experiment ; autumn 1950 
(Wave-length 8°13 m) 


Limits of Total No. of | Hourly No. of 
watch 


UT time echoes rate velocities 


~ 


hm hm 
25 00-06 48 
©5 00-07 00 
©5 00-05 20 
O5 00-07 00 
05 03-07 08 
O5 07-07 10 
05 47-06 14 


600 943 8o 
344 39 
29 6 
564 
57 
653 
64 
Totals 10 43 2824 335 
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The velocity distribution to be expected under such conditions, on the 
assumption that the meteors are moving at the parabolic velocity limit and in 
random directions, has been calculated by Clegg.* ‘This calculation takes into 
account the shape of the aerial beam and the altitude of the apex, and the result 
for the specific conditions of this experiment is shown in Fig. 1 with the measured 
distribution. ‘The experimental errors of measurement have been combined 
with the theoretical distribution as described in Part I. 

There is no sign of a hyperbolic velocity component, but in view of criticismst 
that limitations of technique were responsible for the high-velocity cut-off, a 
further experiment on the antapex was carried out as in Part I. 

4. The 1951 antapex experiment.—'The experiment described above was 
repeated between the hours of approximately 17" and 19" from 1951 March 29 
to May 18. At these tires the antapex occupies the position in the sky previously 
taken by the apex. In 56-5 hours of observation 57 velocity measurements were 
made. The details of the experiment are given in ‘lable II] and of the measured 
velocities in Table IV. ‘The velocity distribution is shown in Fig. 2, together 
with the theoretical distribution calculated as before by the method described by 
Clegg.* There is again no sign of any hyperbolic velocity component, and it is 
evident from a comparison of Figs. 1 and 2 that instrumental cut-off, even if it 
existed, could not be responsible for the low velocity distribution observed in the 
antapex experiment. 


* Clegg, J. A., M.N., 112, No. 4, 1952. 
t Opik, E. J., Irish Astr. 7., 1, 80, 1950. 
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TABLE II 
Data for velocity measurements during apex experiment; autumn 1950 
Amplitude _ AR*? Electron line 
A oo ae density (a) 
(Signal to (Rin units * 107° electrons/ 
noise ratio) of g5 km) cm path 


3 4 5 6 7 


Velocity Range 
V+dV R 
km/sec km 


Time 


ss 


49°8+2'1 11 
50°5+2°4 Il 
69°6+1°7 830 13 
43°0+2°0 10 
53°4+4°6 6 
60°3+2°3 15 
49°8+2°5 5 10 
4r°6+1°9 13 
25°74 1°4 12 
41'2+2°9 12 
53°9+4°3 14 
51'2+3°9 12 
57°8+2°1 14 
53°243°4 14 
58'8+2°0 15 
550+ 1°9 14 
5470+ 1°3 14 
54°5+3°2 10 
55°'2t1'4 15 
60°3+3°6 14 
62°5 + 3°0 il 
53°8+2°6 i2 
37 8t1°9 9 
630+ 3°7 13 
61°'2+3°7 12° 
48°344°1 13 
57°3+4'8 12 
63°2+2°6 14 
64°6+6°6 9 
66-1 + 3°7 12 
53°2+2°0 
51°3+1°8 8 
72°64+2°'8 12 
61°'2+2°3 11 
54°7+2°1 14 
39° 713 13 
63°5+2°7 11 
13 
15 
9 
6 
13 
12 
14 
14 
13 
14 
12 
14 
13 
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TABLE II (cont.) 
Velocity Renee — AR* Electron line 


cue density («) 
— had (Signal to (Rin units * 10" electrons/ 


noise ratio) of 95 km) cm path 
3 4 5 7 


63°1+3°7 580 iI 

70°8+ 38 635 14 ’ 40°0 
69°0+ 4°7 390 13 
601+ 4°0 410 14 
60°1+2°7 13 
43°2+1°9 500 12 
59°4+3°4 390 13 
63°3+2°2 15 
51°2+3°5 13 
43°5+1°8 II 
58°7+5°8 13 
59°2+2°8 13 
66°8+ 2°5 14 
29°7+1°0 13 
63°1+3°9 13 
28:5+1°2 8 
59°24+2°4 14 
702+ 3°5 12 
569+ 3°7 15 
62°14 3°5 14 
60°5 + 2°3 13 
53°4+24 15 
64°0+ 4°3 14 
53°7+2°7 14 
686+ 4°7 10 
55°2+ 31 13 
68:8+7°9 15 
59°2+3°9 6 
64°64 2°7 13 
586+ 2°53 14 
55°4+2°0 14 
47°6+2°8 14 
25°6+2°0 8 
68:9+4°6 12 
71°6+3°2 14 
47°0+2°2 II 
54°74 2-2 14 
55°4+3°4 15 
59°14 4°7 II 
58:2+2°0 14 
49°O+ 2°4 13 
62°5+2°6 13 
645+ 2°5 14 
57°24 3°4 14 
40°9t 2°9 12 
57°74 63 9 
51°5t2°7 14 
53°444°4 Io 
47°5+4°9 8 
49°5+2°3 14 
60°7+ 3°1 12 
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TABLE II (cont.) 
‘ Amplitude AR??? Electron line 
Velocity Range y= : 
VidV R A G density (a) 


onteie wd (Signal to (Rin units X 10° electrons/ 
noise ratio) of g5 km) cm path 


3 4 5 6 7 


61°8+ 3-4 9 
57°2+4°7 540 8 
52°44 2°5 445 II 
63°0+ 3°7 II 
49°3+3°1 14 
52°8+ 2-0 360 13 
62°1+3°5 600 12 
61'1+2°6 13 
481+ 4°2 620 II 
54°94 2°7 365 12 
61°4+4°0 430 12 
62°9+2°'1 435 

55°8+2°'8 620 

62°5+1°8 850 

59°84 3°3 440 

551t31 

568+ 3:2 500 

470+ 2°6 345 

62°0+ 3°7 390 

561+ 3°9 570 

54°7+2°8 

57°5t2°2 

57°8+1°8 560 

30°0+ 1°8 

788+ 571 430 

67°342°2 

39°4+1°7 

52°5+1°9 

58'1+23 

40°7+1°4 

561+ 2°9 

51°3+2°3 

30°5t1°7 

69°8+4"1 

54°8+1°4 

62°5+ 3°6 

54°7+3°1 

80°4+ 81 

62°8+ 6°0 

33°1+3°2 

31°4t1'5 

68-2+4°2 

63°8+2°7 

68°4+ 3:2 

57°8t21 

57°6+ 3:2 

582+ 3°4 

60°3+ 3°0 

47°3£3°7 

51'4t5'1 

62°5+2°7 


DHA OMIM De ow 
- Om Oe OO OW DO On 


6 


PUWAIW CNWONWHNHE 
n O QR eRe Ne eM Ae 





Mary Almond, J. G. Davies and A. C. B. Lovell Vol. 112 


TABLE II (cont.) 


; Amplitude AR*!? Electron line 
ones Velocity Range y= — e 
rime V+dV R A G density (c) 
U.T. fe (Signal to (Rin units 102 electrons/ 
km/sec km . : 
noise ratio) of os km) cm path 


3 4 6 7 


63°8+ 4:0 520 4°0 
61°5+3°0 39 
§2°8+2°4 375 3°5 
7o'1+2°8 640 
266+ 1°9 440 2°5 
35°241°7 
28'0+0'9 
59°8+4°1 
410+ 2°0 
59°9+ 671 
50°9+2°5 
48:2+2°1 
58°8+6'1 
55°0+5°4 
53°1+3°0 
55°5+5°9 
47°9+2°5 
75°344°2 
48°4+ 2°5 
57°446°2 
64°6+ 6:0 
§7°342°7 
5o°2+2°9 
61°7+6°7 
35°94 1°5 
64°1+3°7 
61°7+4'5 
69°6+9°4 
68°4+ 5°5 
60°1+5°6 
56°5+3°5 
43°0+8°7 
59°9+5°9 
7338442 
63°3+7°0 
40°342°4 
*3+1°8 
+ 2°0 
o+3°9 
"3429 
3°5+3°3 
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TABLE II (cont.) 

Amplitude AR®*"3 Electron line 
A ie density («) 

(Signal to (Rin units * 10"° electrons/ 

noise ratio) of g5 km) cm path 


3 4 5 7 


Velocity Range 
V+tdVv R 
km/sec km 


61°7+6°6 500 II 15°8 
66°2+3°2 680 12 . 44°7 
57°8+3°3 13 . 20°0 
49°0+2°3 12 . .15°8 
60°9 + 3°2 6 ° 24°0 
63'2+6°6 8 . 10°0 
550+ 1°6 13 . 16°5 
63°3+2°6 11 ; 13°6 
41°'7+4°2 II ' 19°7 
59°5+3°0 II , 32°6 
58°5+5°5 7 125°0 
42°5+1°5 12 . 17°5 
49°8+2°8 12 ‘ 17°5 
71°5+4°6 12 , 17°2 
50°8+1°9 13 5° 18-6 
64°5+3°3 12 , 17°5 
62°1+3°1 14 ’ 17°5 
39°8+3°8 12 40°8 
640+ 3'0 12 ’ 14°7 
52°342°2 14 . 25°0 
630+ 3°6 14 : 20°4 
43°9t2°4 10 5° 17°9 
62°5+1°5 14 ‘ 17°5 
60°8+ 1°0 14 ; 17'2 
40°2+1°5 13 ’ 20°4 
249+ 13 6 2° 9°3 
60°6+ 2°7 7 35°38 
56°8+1°9 13 , 19°7 
56°8+ 3°4 8 : 19°7 
57°6+3°5 13 : 23°3 
55°7+2°9 14 , 

55:1+2°6 9 
49°3+3'5 8 
34°8+1°9 14 
56°4+5°6 8 
599+ 2°0 13 
52°3+2°6 9 
83°6+3°9 13 
589+ 2°7 13 
50°0+ 56 9 
585+ 6'0 13 
60°5 + 2°6 13 
52°4+5°5 10 
64°5+2°0 14 
63°6+ 3°2 13 
60°1+4°0 10 
37°9+2°2 9 
53°345°1 10 
580+ 2°2 10 
54°6+2°4 8 
41'1t+23 10 
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TABLE II (cont.) 
: Rang Amplitude AR*'* Electron line 
y re R A roe density (c) 


km/sec we (Signal to (R in units * 10" electrons/ 
noise ratio) of gs km) cm path 


3 5 


10 
6 
10 
10 
11 
9 
7 
10 
13 
10 
10 
6 
II 
9 
9 
II 
10 
10 
9 
I! 
8 
10 
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TABLE II (cont.) 


: , ‘ 
Velocity Renge Amplitude puss . Electron line 
V+dV R A G density (a) 
hinine rane (Signal to (Rinunits * 10" electrons/ 
noise ratio) of gs km) cm path 


3 4 5 7 


49°9+5°'0 370 9 § 
55°5+4'8 445 10 ; 12°5 
50°7+2°4 390 10 . 12°9 
55°9+4°9 620 10 4 26°4 
62°34 2°5 500 ; 15°7 
53°5+2°3 330 ' 17°9 
46°8+2°9 490 14°0 
49°4+1°9 360 , 14°7 
61°5+5°7 660 ‘ 34°0 
41°4+3°1 120 ove ove 
53°44 2°0 390 : 14°0 
67°64+2°7 410 ° 13°6 
45°6+2'4 330 : 17°5 
34°2+3°3 410 ; 12°2 
50°9+1°4 442 ; 12°5 
57°1+4°9 600 : 24'0 
48:24+4°3 430 : 10°7 
59°8+2°0 400 ‘ 12°2 
60°7+4°9 370 : 12°5 
565+ 5°1 330 iI : 19°7 
450+ 24 390 10 12°9 
63°1+3°6 450 10 ‘ 12°9 
32°3+2°2 330 10 ‘ 17°9 
69°2+2'1 610 15 36°9 
54°3+3°0 330 16 ; 28-6 
40°4+5°7 310 15 33°6 
06 00 34°8+1°4 140 14 soe eee 
58-9+2°8 560 10 , 19°0 
06 07 36°8+3°9 130 14 ene eee 
06 o9 30°0+ 0'7 120 16 


5. The discrepancy between the theoretical and experimental distributions.— 
Although Figs. 1 and 2 show satisfactorily that there can be no appreciable number 
of meteors with velocities in excess of the parabolic limit, there is, nevertheless, 
a marked discrepancy between the theoretical and experimental curves; the 
experimental distributions lying at much lower velocities than those predicted 
by the theory. A similar discrepancy was remarked upon in Part I in connection 
with the antapex distribution. Since the publication of Part i an error has been 
found in the method of summation used to obtain the theoretical distributions. 
The correct theoretical distributions are similar to those used in the present 
paper. Although the theoretical distribution for the antapex is not changed 
appreciably, that for the apex shows a much greater concentration near the 
parabolic limit. Hence, in the case of the apex results of Part I the experimental 
distributions also lie at much lower velocities than those predicted by theory. 

Thus in all the experiments described in Parts I and II the results not only 
show the absence of hyperbolic velocities, but also indicate that a large 
proportion of sporadic meteors must have relatively short-period orbits. 
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6. The influence of meteors from the helion and antihelion points on the measured 
distribution._A number of meteors included in Tables II and IV were observed 
at ranges less than about 250km. ‘These represent meteors which are not 
observed in the main beam of the aerial system and are not therefore included in 








0 10 20 


Fic. 1.—The measured velocity distribution (histogram) for the 1950 apex experiment. The 
non-hatched areas show the influence of the antihelion component. The theoretical distribution, 
calculated on the assumption of uniformly distributed parabolic meteors, is shown in the smooth 
curve. The theoretical curve is normalized to the equivalent number of observations. 

Ordinates : number of meteors. 
Abscissae : velocity in km|sec. 











10 §=620 30 tC 60 70 80 
Fic, 2.—-The measured velocity distribution (histogram) for the 1951 antapex experiment. 
The non-hatched areas show the influence of the helion component. The theoretical distribution, 
calculated on the assumption of uniformly distributed parabolic meteors, is shown in the smooth curve. 
The theoretical curve is normalized to the equivalent number of (hatched) observations. 
Ordinates : number of meteors. 
Abscissae : velocity in km/sec. 


the theoretical distribution. Since these meteors are observed in a lobe of the 
aerial beam at high elevation in an easterly direction their most probable radiant 
is low in the west; that is, from the antihelion point in the apex experiment, and 
from the helion point in the antapex experiment. Marked concentrations of 
activity from these regions have been observed in the radio echo survey of the 
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sporadic activity.* lor a true comparison with the theoretical velocity distri- 
bution this short-range component must be subtracted from the main distribution. 
Table II (apex) contains 15 of these short-range meteors out of 335 total. The 
influence of this antihelion component in the apex distribution, illustrated in 
Fig. 1, is therefore very small. In the case of the antapex experiment Table IV 
contains 22 short-range meteors out of 57 total, and the influence of this helion 
component on the antapex distribution is marked, as shown in Fig. 2. Similar 
adjustments apply to the distributions measured in Part I. It is evident that 
the removal of these components from the velocity distribution in no way affects 
the conclusions reached about the non-existence of hyperbolic velocities. As 
regards the helion and antihelion components shown in Figs. 1 and 2, calculations 
show that the most probable heliocentric velocity is also non-hyperbolic. 

7. The magnitude range.—In Part | the magnitude range of the meteors was 
estimated by computation of the electron densities in the trails, followed by a 
conversion of these densities to visual meteor magnitudes. ‘This process is 
repeated for the present data obtained with more sensitive equipment. In 
addition, alternative estimates from a comparison of the visual and radio echo 
rates are made both for the previous results in Part | and for the measurements 
described in this paper. 

(a) Magnitude estimates from electron densities in the trails.-For a received 
signal of strength «, and for an apparatus with parameters: wave-length (A), 
transmitter power (P) and aerial gain (G), the line density of electrons in a meteor 
trail at range R is given by t 


siesta 2 RS 
=Vi20 (") a Je electrons/cm path. (1) 


This only applies to a trail where the volume density is less than the critical 
density for the radio wave-length employed. Discussions of the situation which 
arises tor higher electron densities have been given by Herlofson{, Feinstein§, 
and Kaiser and Closs.|| From their results it is evident that for the electron 
densities with which we are concerned in this work, no appreciable error will be 
introduced by using equation (1) to compute the line densities from the measured 
signal strengths. 

As in Part | a quantity y = AR**/G is calculated for each meteor and is given 
in column 6 of ‘Tables I] and IV. Here, A is the amplitude of the radio echo 
from the first Fresnel zone compared with the receiver noise level, R the range 
of the meteor and G the effective gain of the aerial system. G depends on the 
range R and on the azimuth of the meteor with respect to the main beam of the 
aerial system. The method of computing y is described in the Appendix to 
Part I. The only additional information required in this paper is the appropriate 
sensitivity-range curve for the modified aerial system, which is given in Fig. 5. 

The line density « for each meteor is then calculated from 


— 2 ae 


* Hawkins, G. S. and Aspinall, A. (unpublished). 

t+ Lovell, A. C. B. and Clegg, J. A., Proc. Phys. Soc., 60, 491, 1948. 

} Herlofson, N., The Observatory, 68, 230, 1948; Arkiv for Fysik, 3, 247, 1951. 
§ Feinstein, J., 7. Geophys. Res., §6, 37, 1951. 

|| Kaiser, T. and Closs, R. L., Phil. Mag., 43, 1, 1952. 
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TABLE III 
Details of antapex experiment; spring 1951 
(Wave-length 8-13 m) 
Date Limitsofwatch ‘Total No. of | Hourly No. of 
1951 1 i i echoes rate velocities 
hm hm 
Mar. 29 17 16-17 50 : 6 10°6 
30 17 00-18 58 27 
31 17 03-18 58 15 
Apr. 2 17 10-18 54 3 
4 17 15-18 58 32 
5 17 08-19 02 21 
6 17 00-19 00 13 
9 17 00-19 03 21 
10 17 00-18 57 
II 17 21-19 00 
12 17 04-19 00 
13 17 19-19 03 
14 17 10-19 02 
16 17 01-18 55 
17 17 00-19 00 
18 17 00-19 00 
19 17 00-19 00 
20 17 00-19 00 
17 28-19 00 
17 02-18 52 
17 OI-19 O1 
17 03-19 00 
17 09-19 04 
18 23-19 00 
17 00-19 00 
17 17-19 00 
17 13-19 00 
17 00-19 05 
17 00-19 00 
17 37-19 00 
17 04-19 00 
17 32-19 00 
17 10-19 05 
17 45-19 00 
Totals 56 33 1586 
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TaBLe IV 
velocity measurements during antapex experiment ; spring 1951 

: 3/2 : 

Velocity Renge Amplitude ee. Electron line 
V+dV R - A G density («) 
km/sec “val (Signal to (Rin units * 10"° electrons/ 
noise ratio) of 95 km) cm path 
3 4 7 


37°4£2°0 : 3°6 
35°7t0°9 350 ‘ 10°7 
35 ‘ot 1°6 395 ° 10°00 
35°8+2°4 45° 7°5 
23°7+1'2 140 ose eee 
27°5+2°6 ° 5:0 
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TABLE IV (cont.) 


Veloci Amplitude AR* Electron line 
elocity Range = ‘ 
V+dV R A G density (a) 
- (Signal to (Rin units X10" electrons/ 
noise ratio) of g5 km) cm path 

2 3 6 7 

hm 
18 40 26-5+2°1 6°5 
17 00 24°6+1°9 2°1 75 
18 29 34°64 3°2 10°O 35'8 
27°9+0°8 3 I'l 3°9 
36°6+2°8 4°2 
18:2+1°6 
551° 
30°8+1 
31641 
39°O + 3° 
37°64 1° 
20°7+0 
24°0+0° 
21-°1+o0 
28:8+1 
31°32 
35°52: 
18°8+0° 
19°6+1° 
33°9+2° 
28-9+1 
26°3+1 
24°52 
315 t2 
26°6+ 1° 
26°5+0° 
389+ 1° 
25°5+0° 
2QrO Tt 1° 
32°94 2° 
19°9+0° 
35°8+3° 
35°9+2° 
27°25 
30°4+2° 
22° + 3 
360+ 1° 
49°8+ 4° 
30°6+1° 
22°7+0° 
27°7+1° 
32°0+2° 
44°4$ 2° 
37°O+ 14 
26°2+2°4 
40°7+3°4 
35°st29 
36°5+1°6 
25°4+0°8 
31°0+0'8 
40°5+1°0 
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using the values of A, P and «g given in Section 2; and is given in column 7 
of Tables I] and IV. ‘The corresponding number distribution against electron 
density for the apex and antapex experiments is shown in Figs. 3 and 4. 

The difficulties in the conversion of these line densities to visual magnitudes 
have been discussed in Part I. The most recent published data on which the 
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Fic, 3.— Number distribution of electron line density («) for the 1950 apex velocity measurements. 
(Fifteen meteors for which there is no information about «, and one meteor for which the only 
information is «> 18-6 x 10° are not included. Eight other meteors with values of x= 104, 109, 
125, 165, 193, 238, 247, 286 X 10! lie outside the limits of the abscissa.) 

Ordinates : number of meteors. 
Abscissae : line density « electrons/cm path ( x 1°57 X 10"). 
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Fic. 4.—Number distribution of electron line density (a) for the 1951 antapex velocity 


measurements. (Twenty-two meteors for which there is no information about « are not included.) 
Ordinates : number of meteors. 


Abscissae : line density = electrons/cm path ( X 1°57 < 10"). 
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conversion may be based are given by Herlofson*, according to whom the peak 
of the distribution curves in Figs. 3 and 4 lies at zenithal magnitude 3-5. However, 
recent experimental and theoretical workt has shown that a considerable 
revision of these estimates is necessary, and that the ionization produced by a 
meteor of a given magnitude is about 40 times greater than that given by 
Herlofson. Thus the peak of the distribution curves in Figs. 3 and 4 must lie 
in the region of magnitude 7. 


6°5 

















100 200 300 400 500 600 700 800 900 


Fic. 5.—Relation between most probable sensitivity S and range R, for the modified 8-metre 
aerial system. From these curves the values of y are calculated for each meteor. 
Ordinates : most probable sensitivity S. 

Abscissae : range in kilometres. 


In order to examine the velocity distribution for changes related to electron 
density the 335 velocities comprising Fig. I (with the 15 short-range meteors 
rejected) have been divided into three groups covering different ranges of a. 
The three distributions are shown in Fig. 6. ‘There is no indication of any 
change in velocity distribution with magnitude. 

(b) Magnitude estimates by comparison of visual and radio echo rates.—This 
method depends on a comparison of the relative collecting areas for meteors of 
different magnitudes for a visual observer and for the radio equipment. The 


* Herlofson, N., Phys. Soc. Rep. Prog. Phys., 11, 444, 1948. 
+ Unpublished work at Jodrell Bank. 


3* 
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estimates of the observer’s collecting area for various magnitudes have been based 
on the limiting brightness of a meteor visible at a given angle from the line of 
vision. ‘The available data on this, taken from the work of Backhouse*, Opikt, 
and Hoffmeister{, are plotted in Fig. 7.§ From these data values have been 
adopted as given in Table V. 








0 10 80 


Fic. 6.—The oolecity distribution for different values of line density (x) for the 1950 apex experiment. 
a<14°4X 10!° electrons/cm path (99 velocities). All 
a 14°4<4< 22°2 X 10" electrons/cm path (124 velocities). normalized 
a” >22'2X 101° electrons/cm path (96 velocities). to 100 total. 
distribution for all as. 
Ordinates ; number of meteors. 
Abscissae : velocity in km|sec. 
(One meteor for which the only information is «> 18-6 x 10! is excluded, together with the 15 
short-range meteors for which no values of « are available.) 


A 
ies 


If the observer’s line of sight is directed upwards at 45°, the intersections of 
these cones of limiting visibility on the surface 95km above the Earth can be 
plotted. ‘These apparent magnitudes are corrected to zenithal magnitude by 
using data given by Porter|| and Opikt+, which are in agreement to within a quarter 
of a magnitude. From the above data the areas in square kilometres over which 
one observer can see meteors of a given magnitude have been measured graphically. 
The results are given in column 2 of Table VI. Column 3 gives the relative 
number of meteors in the various magnitude groups falling on a square kilometre 

* Backhouse, T. W., The Observatory, 7, 299, 1884. 

+ Opik, E. J., Publ. Tartu Obs., 25, No. 4, 1923. 

t Hoffmeister, C. (private communication). 

§ In plotting Hoffmeister’s data it has been assumed that a meteor of zero magnitude can just 


be seen at 60° from the line of vision. 
|| Porter, J. G., Mem. B.A.A., 34, Pt. 4, 1942. 
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of the sky. ‘These are based on the measurements of Watson* which show that 
the average increase in numbers of meteors is 2-5 times for each fainter magnitude. 
Thus the relative numbers (N,,N,) between two magnitude groups (m,, m,) 
will be given by 

N, / Ne an (2-5) ma). (2) 








0 10 20 30 40 50 60 70 80 
Fic. 7.—Relation between limiting brightness of a meteor visible at a given angle from an 
observer's line of vision according to Hoffmeister ( ); Opik (——-) and Backhouse (- « - * -). 
Ordinates : meteor magnitude, 
Abscissae : angle from observer's line of vision. 


TABLE V 
Limiting meteor magmitude visible at a given angle from an observer's line of vision 
Angle from line of vision 60° 58 53° 44 30 ae’ °° 
Limiting magnitude ° I 2 3 4 5 5°5 
TABLE VI 
Numbers of meteors of different magnitudes seen by a visual observer 
‘Prue number 
per sq. km 
per hour for 
hourly rate =1 


Zenithal Collecting Relative Relative 
area number number 


magnitude (sq. km) per sq. km seen 


2 3 4 
( x 1078) 
161 158 25 000 _ grs 
5 160 100 516 000 20°2 
9 845 63°1 620 000 12°7 
16 125 39°5 641 000 8-02 
24 200 25°! 607 500 5°05 
36 150 15°8 571 000 3°18 
50 000 10°0 500 000 2°02 
66 000 6°31 420 000 1°27 
87 500 3°98 348 000 0-80 
115 700 2°51 291 000 0°50 
135 800 1°58 214 000 0°32 
200 000 1°00 200 000 0°20 


Total 4 953 500 


* Watson, F., Harv. Ann., 105, 624, 1937: Proc. Amer. Philos. Soc., 81, 493, 1939- 
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If one meteor of magnitude —1 falls per sq. km per hour, the number of 
meteors seen in each magnitude per hour can be calculated, and is given in column 4. 
Summation of this column gives the total rate seen by the observer. If this is 
taken to be one meteor per hour, the relative numbers in column 3 can be scaled 
down to give the actual numbers of meteors falling on a square kilometre per hour, 
as in column 5. 

Assuming that the apex is at an elevation of 45° (its mean elevation during the 
apex experiments) and that the meteors fall on the Earth from random heliocentric 
directions with parabolic velocities, it is possible to deduce the numbers of meteors 
passing per hour through a square kilometre of the 95 km surface from different 
directions. ‘These true number$ may now be applied in place of the relative 
numbers given by Clegg* in calculating the theoretical velocity distribution and 
hence, making the arbitrary assumption that meteors of a certain magnitude are 
seen to a given sensitivity contour, the number distribution with magnitude can 
be plotted. ‘This magnitude distribution has a peak at + 0-7. 

With the arbitrarily assumed visual rate of one meteor per hour, and for this 
magnitude distribution, the equivalent radio echo rate derived in this way came 
to 0-021 per hour for the 1950 apex experiment described in Section 3. If a 
visual rate of 10 per hour is adopted for the early mornings in autumn when the 
experiment was performed, the equivalent radio rate becomes 0-21 for this magni- 
tude range. ‘The actual observed magnitude range can now be estimated by 
comparing this rate with the hourly radio echo rate during the experiment. 
The total rate observed during the apex experiment was 264 per hour. The 
shift in the magnitude scale required to bring this rate into line with the above 
results can be derived from equation (2) which transforms to 


mM, — Mz = 2°5 logyg(N,) No). 


The shift required is therefore given by Am =2-5 log,)(204,0°21)=7:7. ‘The 
assumption of a visual rate of 10 per hour is not critical. For example, a doubling 
of the assumed rate would change the magnitude shift by only 0-75 magnitude. 
‘lhe assumed distribution peaked at magnitude 0-7, so the observed distribution 
peaked at magnitude &-4. 

Of these echoes only 11 per cent gave velocity measurements; although an 
additional 29 per cent were of sufficient amplitude and duration they did not 
yield velocities owing to distortion of the trail or lack of range measurement. 
Thus the rate of velocity type echoes is 40 per cent of the total and the 
corresponding magnitude distribution peaks at magnitude 7-4. 

The same analysis can be applied to the case of 1949 apex experiment described 
in Part I. In this case the peak of the magnitude distribution for the velocity 
type echoes is at magnitude 6-1 and for all echoes is at magnitude 7-1. In the 
case of the antapex experiments, the calculation has not been attempted because 
of the difficulty of obtaining a reliable visual rate, and also because a large pro- 
portion of echoes occurs outside the main lobe of the aerial. 

In Table VII the peaks of the magnitude distributions in the various 
experiments derived from visual rates are compared with those calculated 
from the peak electron densities as described in (a) above. 

In Part I it was concluded that the 1949-1950 measurements probably 
referred to meteors of about magnitude 4:5. ‘The estimates now made indicate 


* Clegg, J. A., M.N., 112, No. 4, 1952. 
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that these meteors were of fifth and sixth magnitudes. In the 1950-1951 
experiments described in this paper the meteors are fainter on both estimates 
by 1-5 magnitudes. The uncertainty in estimating the magnitudes arises on 
the one hand because of the incomplete state of the meteor evaporation theory, 
and on the other because of the lack of data on the collecting area of a visual 
observer for meteors of a given magnitude. It is therefore only possible to 
conclude that the measurements described in Part I probably referred to meteor 
magnitudes between 4:5 and 6:0, and in Part II to meteor magnitudes between 
6-0 and 7°5. 
TaBLe VII 
Comparison of magnitudes from electron densities and from visual rates 
From visual From electron density at peak of 
rate distribution curve 
(a) (b) 
1949 apex 6:0 5°5 2°0 
1950 apex a 7° 3°5 
(a) Assuming the new estimates for electron density (see page 35). 
(6) Using Herlofson’s estimates for the ionization produced by a meteor 
of given magnitude. 


8. Conclusion.—The experiments described in this paper and in Part I give 
no evidence that any appreciable numbers of meteors are moving with velocities 
in excess of the parabolic limit. ‘This conclusion is in conflict with the work of 
Opik*t and Hoffmeister.{ Although it appears possible to reconcile Hoffmeister’s 
results with this conclusion by abandoning the assumption of uniform distribution 
(see Part I), there does not appear to be any ground on which the results of 
Opik can be aligned with those presented in this paper. ‘lhe majority of the 
measurements probably refer to meteors in the magnitude range 6-0 to 7:5, with 
limits of 7-5 tog-0. ‘The results therefore cover not only the visual range available 
to Opik, but also the telescopic results of Boothroyd.§ Although in his earlier 
analysis* Opik concluded that the velocity distribution changed with meteor 
magnitude, his final decisiont, from the analysis of both the 'T'artu and Arizona 
work, was that no significant change in velocity distribution with magnitude 
existed. In this respect only are the present radio echo results in agreement 
with the work of Opik. 
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ABSORPTION OF RADIATION BY AN ATMOSPHERE OF 
H, H+ AND H{—SEMI-CLASSICAL TREATMENT 


D. R. Bates 
(Received 1951 September 238) 


Summary 


Calculations are carried out on the absorption of an atmosphere of H, H* 
and H} using the accurate quantal dipole moment but treating the nuclei as 
classical particles. Absorption coefficients are tabulated for wave numbers 
between soocm~' and 26000cm~', and for temperatures between 
2 500 deg. K and 12000deg. KK. The relative importance of bound-free 
and free-free transitions is also determined. 





1. /ntroduction.—Some time ago Wildt (1) suggested that a significant 
contribution to the opacity of stellar atmospheres might arise from bound-free 
absorption by Hj ions, 

H;(1so)+hv-+-(H + H* \(2po) (1) 
and from free-free absorption by colliding H atoms and H~* ions 
(H + H*)(1s0) + hv-+(H + H*\(2po). (2) 
The present paper is devoted to the investigation of these two processes on the 
assumption that the nuclei behave as classical particles. 

2. Basis of calculations.—The absorption coefficient «(v, 7’) associated with 
(1) and (2) together is related to the photon emissivity, j(v, 7), due to the inverse 
processes, by the well-known equation 


K(v, T) =(c?/8rv*){exp (Av/RT)—1hj(v, T), (3) 

v being the frequency of the radiation, T the temperature, 4 Planck’s constant, 
and k Boltzmann’s constant. Formulae from which j(v, 7) may be obtained 
have been developed elsewhere (2). On introducing »(v, T') dv, the rate coefficient 
giving the photon emission in the frequency range v to v + dv, we see that 

(», T) =n(H)n(H*)y(x, T), (4) 
the n’s being the concentrations of the substances indicated. Further, if R is 
the internuclear distance corresponding to the frequency v, and y(R, T)dR is 
the rate coefficient giving the photon emission from H atoms and H* ions whose 
separation lies between R and R+dR, then 

vv, T)=(R, T)(dv/dR). (5) 
Now from the earlier paper (2) just mentioned 


yV(R, T) =47gA(R)R?* exp { — U(2zpo| R)/RT}, 


where g is the ratio of the statistical weight of the repulsive potential energy 
curve U(2po|R) to the sum of the statistical weights of this curve and of the 
attractive potential energy curve U(1tso|R), and where A(R) is the Einstein 
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spontaneous transition probability. Denoting the dipole moment involved by 
D(R), writing 
v(Iso| R, T)? = — U(1s0| R)/kT, (7) 
x(2pa| R, T)? = — U(2po| R)/kT, (8) 
and using the relations 


“ 


A(R) = D(R), (9)* 


hv=L ame — U(1so| R), (10) 
we tind 


u(v, 1) = 2 ree ‘lar +R D(R)lexp {x( 180] R, 7)?| — exp {x(2pa] R, T2}]n(H)n(H*). 
(11) 


This formula can of course be derived directly using classical statistics to obtain 
the distribution of internuclear distances. On taking g to equal 4, and substi- 
tuting for the fundamental physical constants, we have 

«(v, T) = 2°51 x 1074 {y/(dv/dR)}R*D(R)lexp {x(1so| R, 7)*} 

— exp {x(2po| R, T)?}|n(H)n(H*) cm ', (12) 
where R and D(R) are in the appropriate atomic units (a) and a,"e? respectively), 
n(H) and n(H*) are particle numbers per cm*, and the other quantities are in 
any self-consistent set of units. ‘The stimulated transitions are of course 
included in «(v, 7’) as negative absorptions. 

An asymptotic expression suitable for use at low frequencies can be derived. 
From elementary perturbation theory (3) we have that as R tends to infinity 


vie~(8 3)(R—3/2R)exp( — R) Rydbergsf, (13) 


v/(dv dR) ~ R(2R? — 3) (2R® — 2R? — 3R-—3). (14) 
Noting also (4) that 
R) ~4R%, (15) 
we see on expanding the exponentials in (12) that 
hv R°(2R? — 3) 
) ~6: x it oo + —1 
k(v, T) ~ 6°28 x 10 (ar) apse aa fia jem", (16) 


the units being as specified above. From (13) R is approximately 





— In(3v/8c) + InkR, (17) 

so that the value appropriate to any chosen v can be found rapidly by iteration. 

For certain purposes it is desirable to be able to separate the contributions 

from bound-free and free-free absorption. ‘This can readily be done, for, as 

can be seen from the earlier paper (2), the fraction of the total due to the former 

is given by 

rp(v, T)=@{x(1s0| R, T)} —27 !x(1so| R, T)exp {—x(1so|R, T)*}, (18) 

where ® is the error integral. At low frequencies r,(v, 7’) is clearly very small. 

An asymptotic expression can be obtained, but it is of little interest and need 
not be given. 


* The orbital degeneracy factor is omitted as it equals unity for the transition under discussion 
+ 1 Rydberg= 109 677°58 cm™. 
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All the functions appearing in (12) and (18) are known with some precision ; 
thus the potential energy curves have been evaluated in the course of a recent 
quantal study of the two-centre problem (§), and the dipole moment has been 
calculated (6) from electronic wave-functions made available by this investigation. 
A number of convenient values of the wave number v/c covering the range 
500cm'! to 26000cm™! were chosen, and the corresponding values of R, 
U(1sc|R), U(2pe|R), v(dvdR) and D(R) were computed by inverse inter- 
polation and numerical differentiation. ‘The figures obtained (which are 
displayed in Table 1) were then used in (12) and (18) to determine «(v, T) and 


TABLE I 
Data on 1sa—2po transition of Hy 
Internuclear 
Wave number __ distance v/(dv/dR) Potentials Dipole moment 
vic R (atomic U(iso | R) + U(a2po | R) D(R) 

(cm~') (atomic units) (Rydbergs) (Rydbergs) (atomic units) 
units) 

500 8-633 

I 000 7°839 

I 500 7°367 
2 000 7°026 
2 500 6°760 
3 000 6°5416 
3 500 6°3548 
4 000 "1QIg 
5 000 OI79 
6 000 O91, 
7 000 
8 000 
9 000 
10 000 
12 000 
14 000 
16 000 
18 000 
20 000 
22 000 
24 000 
26 000 


15 0°00315 O°OoI4! 17°4 

16 ‘00583 "00329 14°2 

176 00843 "00525 12°5 

“18, ‘01097 "00727 I1°3 

‘197 "01346 "00933 10°4 

"206 ‘O1592 ‘O1143 9°77 
‘215 "01834 "01357 Q's 
"229 02074 ‘01573 8°76 
‘236 "02545 "02014 7°92 
‘249 ‘03004 "02467 7°31 
"328, 27) ‘03891 "03403 6°37 
285 "04322 "03884 6°09 
"2B, "04744 ‘04374 
"30% "05560 "05381 5‘la 
"315 ‘06347 06418 4°71 
"325 ‘07100 -07488 4°36 
‘33s ‘07828 ‘08584 4°05 
"340 ‘08527 ‘09708 3°79 
345 ‘09200 "10859 3°56 
351 ‘09847 "12035 3°36 
"356 0*10470 0°13236 3°Is 


uuu D> 


5°64 


I 
I 
I 
I 
1 
I 
I 
I 
I 
I 

49798 1°26, 03451 02931 6°85 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 


ry(v, 7) at various temperatures between 2500deg.K and 12000 deg. K. 
‘Tables II and III contain the final results. ‘Those in the former table may be 
a few per cent in error (particularly at the lower wave numbers) because of 
possible slight inaccuracies in the original data; those in the latter should be 
correct to the number of decimal places given. _ It is difficult to assess the extent 
of the uncertainties introduced by the assumption that the movement of the 
nuclei is strictly classical. The de Broglie wave-length of a proton near the 
maximum of the Maxwellian distribution at 2 500 deg. K (the lowest temperature 
considered) is about one atomic unit, which is certainly not negligible compared 
with the internuclear separations of importance; but in spite of this the classical 
treatment of the free states should actually be quite successful since the use of 
the Winans~Stueckelberg 5-function (7) gives a fair representation of the photo- 
dissociation of even light molecules (8). The quantization of the bound states 
is unlikely to be very important owing to the large number involved. We consider 
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TABLE II * 
Absorption coefficient in units 10-** em—"/(H atom/em*)/(H*+ ion/cm'*) 
‘Temperature in 10° deg. K 
2°5 3°5 4 5 6 7 10 12 
1°14 ’ 080 40°69) «60°55 046 039 0'270 0°226 
1°61 I'I2 097 0°77 63 "54 . "38 “gt 
1°35 117 0°92 ‘76 64 “44 ‘37 
1°53 "32 1'03 “85 "72 . "50 “41 
1°69 1°45 "13 ‘92 ‘78 ‘54 44 
184 1°56 ‘21 0°99 83 ' ‘57 “47 
I "O04 88 ‘ ‘60 “49 
*Q2 ; 62 “51 
98 66 54 
4 69 ‘57 
08 . ws 58 
"12 : ‘73 59 
“16 : "74 ‘60 
‘19 , "75 ‘61 
25 , "77 62 
-30 ‘ 78 -62 
36 . 78 62 
“41 ; "719 “61 
46 "14 "79 ‘61 
"52 “16 ‘79 ‘bo 
‘57 18 ‘79 59 
63 ‘21 0'79 0°58 


3 500 
4 000 . A 2°10 ‘77 
5 000 : . 2°36 ‘96 
6 000 : : 2°63 15 
7 900 es ; 2°91 = 2°33 
8 000 ; 3°2 2°53 
9 900 ; 3°5 2°74 


‘97 ‘67 I 
I 
I 
I 
I 
I 
I 
10 000 : 5° 3°9 295 ‘ - 
I 
I 
I 
I 
2 


— me 


12 000 ‘ 4°7 3°4 

14 000 , 57 8640 

16 000 . 70 4°7 2°64 

18000 35 . 84 54 2-91 

20000 47 7 10°2 6°3 3°2 03 
22000 64 : 12‘2 7°3 3°5 2°16 
24000 86 14°6 8-4 3°8 2°29 
26 000 114 17°3 9°6 4°2 2°42 


emt 


* At wave numbers below 500 cm~! the asymptotic formula (16) may be used. 


Tabet III 
Fraction of total absorption due to bound-free transitions 
‘Temperature in 10° deg. K 
3 3°5 4 5 6 7 

0°044 0°037 0°031 0°022 O'017 O'OI4 
"107 ‘087 "072 053 ‘O41 033 
‘171 “141 "118 ‘o88 ‘069 056 
+2360 *196 *166 "125 “098 ‘o80 
298 ‘250 "214 "162 "129 “105 
°357 "303 ‘260 "200 ‘159 ‘131 
“413 ‘353 *395 ‘237 "190 157 
464 *400 349 ‘273 ‘221 "183 
556 486 "429 "342 ‘280 234 
632 561 ‘501 406 330 *284 
696 625 "564 464 388 331 
‘748 «680 619,517) 437 °375 
‘792, -°727,- 667, *564 “482 “417 
827, 767, *709--_— 607. 524 “456 
“881 829 "777 ‘680 ‘597 "526 
‘O17 874 829 “739 658 586 
"942 ‘906 867 -786 ‘708 638 
‘959 930 “897, *8240—Ci*751 “683 
‘970 «947-0 ‘919 *8540— *786 “721 
978 ‘960 936 879 816 "754 
‘984 “969 ‘949 *898 “841 -782 
0°976 + 0°959 0'914 0°862 0806 


‘ ry 

wwNnNN 
wn wu OM 3 
8888s 
0oOwW 
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therefore that the results presented are probably accurate at the higher tempera- 
tures and lower frequencies and that they are probably correct to well within 
one part in ten even at the lower temperature and higher frequencies. Quantal 
calculations in the former region would be very laborious. ‘Though still very 
formidable in the latter region they have, nevertheless, been undertaken by 
Dr R. A. Buckingham* and Mr S. Reed and should show clearly the range of 
validity of the classical approximation. ‘The two investigations are essentially 
supplementary and together should provide reliable data at all temperatures 
and frequencies of interest. 

The application of the absorption cross-section data to stellar atmospheres 
is being carried out by Professor B. Strémgren. 


Dept. of Applied Mathematics, 
Queen’s University, 
Belfast, N. Ireland : 
1951 September 26. 
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ON EMISSION LINES IN STELLAR SPECTRA 
Ida W. Busbridge 
(Received 1951 October 27) 


Summary 


Conditions are investigated for the formation of absorption lines when 
the fraction € of the line absorption coefficient /, is due to absorption and the 
fraction 1—€ is due to scattering, continuous absorption and continuous 
isotropic scattering being taken into account when estimating the flux of 
continuous radiation. It is shown that if << A,=«,/(x,~+-a) the line must be 
dark, but that if €> 4A, the line may be bright. The case €=1, in which the 
line is formed by pure absorption, is considered in detail. It is shown that 
the line will be bright if f(A,, By) >f(A,, By), where f(A, f) is a certain function, 
Ay=(K»t+l,)/(k»+o+/,), and f, is a constant depending on the constitution 
of the star and the frequency v considered. Some numerical examples are 
discussed. 





1. Introduction.—In a recent paper (3) Miss Underhill has considered the 
formation of absorption lines when continuous absorption and continuous 
isotropic scattering are both taken into account in estimating the flux of 
continuous radiation, and when the lines are formed (a) by absorption and (6) 
by scattering. Assuming that o, the coefficient of continuous scattering, is 
independent of both the frequency v and the depth, and that the coefficient of 
absorption «, is independent of the depth, she has shown that lines formed by 
scattering must be dark. She has also shown by an extreme case that lines formed 
by absorption may be bright, but she merely says that in other cases ‘the 
appearance of emission becomes more stringent ”’. 

The main object of this paper is the investigation of these more stringent 
conditions, but I have begun by considering what happens when the line is 
formed partly by absorption and partly by scattering. 

2. The residual intensity.—Let €1,(0 <€ <1) be the part of the line absorption 
coefficient /, which is due to absorption and let the remainder (1 —€)/, be due 
to scattering. ‘Then the equation of transfer is (with the usual notation) 


dI ra 
some ie — — — 
H de Moto +h, Slo + (1 El] | Ted (5+ EL) By (2.1) 
where x is the depth below the surface of the star of the layer considered. 
On writing 
AO Kk, +él, 
: k,totl, 
and 
t,= | p(x, +o +1,) dx, 
“0 


(2.1) takes the standard form 
dl, 


“1 
i a ee ee) 
Th. b( )\ ye 





46 Ida W. Busbridge Vol. 112 


Following previous writers, Underhill assumed that, for the surface layers, 
a linear approximation for the Planck function B, is sufficiently good. Accepting 
this, though remembering that the exact solution of (2.4) depends on the values 
taken by B, for all values of t, from 0 to 0, we write 


B,=b +6,7,, (2.5) 
where — J ” («,-+0)p de. 
0 


If /,, «, and o are all independent of the depth, then 
T,=N,t,, (2.6) 
where 
K,+o 
n,= Soret’ (2.7) 
Subject to these conditions, (2.4) can be solved exactly (2) and the emergent 
flux is 
F,= 26, {n, R(A,) + (6,0/b,)R(A,)}, (2.8) 
where 


Ry(A) = Mag + HT Aon’, (2.9) 
R,{A) =A"*a,, (2.10) 
and a, % are the first and second moments of H(), the solution of the equation 


Hu) = 1+ 4A) [AOD ay (2.11) 


By a theorem of Chandrasekhar*, a, and a, are related by the equation 
N2a5 + 4(1 —A)a,? = 4, (2.12) 
so that R,(A) can also be written 
R,(A) = 4(1 —A)a,? + 4. 2.13) 
Tables for R,(A) and R,(A) are given by Underhill (3) and tables for «, and «, 
by Chandrasekhar (2). 
The emergent flux of v-radiation in the continuum is obtained by putting 
1,=0 in the preceding equations. Writing 
A. =k,/(k, +e), -I4) 
this is 
Front. = 26, {Ri (Ae) + (0,0 /b, RA.) }- 
‘The residual intensity R, in the absorption line is 
a F, RA.) + (b,0/6, RA) =e 
Fron. RA) + Os] 0RA) es 
This is identical with the expression obtained by Underhill (though her A, and 
A,’ are obtained by taking §=1 and €=0 in (2.2)). The form (2.15) is, however, 
not the most useful form for R, and we proceed to transform it. 
By (2.2), (2.7) and (2.14), 





n, =(§—A,)/(E-A,) (2.16) 


* (x), theorem 2. Here ¥(p)=$(1—A). 








No. I, 1952 On emission lines in stellar spectra 47 


except when €=A,. ‘Then =A, also. Hence R, can be written, when & AA,, 

in the form 

(€-A,)R, (A) + BORA) — 
(€-A)R(A) + BORA)” 


a a (€ —A,)b™, 6.™. (2.18) 
We shall tirst show that if <A,, then R,<1 and the line must be dark, but 
that if >A, the line may be bright. 
Since 0<n,<1, it follows from (2.16) that if €<A,, then <A, and 
E<A, <X.. (2.19) 





R,= 


where 


If €>A,, then >A, and 
E>A,>DA,. (2.20) 
When é<A,, 8,“* is negative and R, can be written 
— AM -ORA) + |B, RAM) 
"Ae €)Ri(Ae) + 1B | RalAc) 
From the tables given by Underhill for R,(A) and R,(A) it is at once seen that 
R,(A) is a decreasing function in (0,1) and that R,(A) is an increasing function. 
Moreover a simple calculation shows that AR,(A) is increasing in (0,1). It 
follows that 





(A—€)Ry(A) + 1B, RA) 
is increasing in (0, 1) and, since <A,“ <A,, that R,<1. ‘Thus when £<A, the 
absorption line is dark. In particular, lines formed by pure scattering (€ =0) 
are necessarily dark. ‘This is the result established by Underhill. 
When >A, we return to the form (2.17) for R,, since 8,“ is now positive. 
Since (€—A)R,(A) is a decreasing function of A while R,(A) is increasing, 
(E-N)RQ) + BRA) 
will be increasing or decreasing in (0,1) according to the values of & and £,“*. 
Thus it is possible to have R,>1, so that the line is bright. 
We shall discuss only the case of pure absorption (€=1) in detail. Let 
J(A, B) = (1 —A)R,(A) + BRA). (2.21) 
By (2.10) and (2.13), this can be written 
fs B)= (tA)? + JAA) + BM. (2.22) 
From (2.17) with €=1 and (2.21) we see that 
R, =f0.", BL) flo BO). 
Since, by (2.20), 1>A,‘”>A,, the condition for the formation of an emission 
line can be written 
Sv Be) >flrer Br) (Av>Ac)s (2.23) 
where, by (2.2), (2.14) and (2.18), 


Kk, +l, as 
K,t+o+l,’ ° x,+0 


A, =A, = - 


and 


B,=B, =(1—A,)b,/6,. (2. 


25 
Note that «, in (2.22) is a function of A, since H(z) is a function of p and A 


) 





48 Ida W. Busbridge Vol. 112 


3. The constant B,.—The usual definition of 6,6, is 


b5,© my hig! f4T\— 
5,0 =|I —exp(- bv TIE (Fe) 


7T=0 


where & is a mean absorption coefficient and 


a 
T= | Kp dx. 
0 
Assuming the approximate formula 
T4 = T,*(1 +7), 


where, for a grey atmosphere, a = 3/2, we have 


dT 
(=)... = jaT,. 


Hence, if x=hv/kT,, 


Bp 


:= 


x=. 2 
kT, (3.2) 
In this form f, is completely determined by the frequency v considered and by 
the constitution of the star, subject to the assumptions made above. 
It is easy to see that £8, lies between o and 4, assuming that a, which is 3 for 
a grey atmosphere, cannot be less than unity. Since 
es. # 


Pils oe Gat oe =, 4° aa 
x"(1-e°-7)=1 ats? as x>O 


d f 

dx ' 
therefore (1—e~*)/x decreases steadily from I to 0 as x increases from 0 to %. 
‘The ratio o/(* +o) of the continuous scattering to the total continuous absorption 
can take any value between 0 and 1, and therefore £, lies in the stated range. 

We have to discover the range of values of 8, for which (2.23) can be true. 

Since x (i.e. hv/kT)) decreases with increasing wave-length, for a given star 
8, will increase with increasing wave-length, and the largest value of 8, for visible 
light will be that corresponding to the red end of the spectrum. From this it 
will subsequently appear that emission lines in high-temperature stars are more 
likely to appear in the red part of the spectrum, a result which accords with 
observation. 

4. The function f(A, 8).—The values of f(A, 8) for values of A lying between 
© and 1 and for certain values of f are given in Table I and are shown by the full 
lines in Fig. 1. 

When § <1, the behaviour near A =0 needs careful investigation, and for this 
we have to go back to the behaviour of H(,) for small A. 

Let H,() be the solution of (2.11) whenA=0. ‘The exact solution of Milne’s 
problem in the grey case is /(0, 4) =} 3FH)() and for this reason the second 
moment of H,(u) is known with great precision, while the first is known exactly. 
‘These are 


x (1-e°*)}=—x 2e *(e’—1-x)<o for x>0, 


a;%=2/4/3, 2 =0-820 352. (4.1) 
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TABLE I 


Values of f(A, B) 





0°710446 ° : ° “ . 3258889 





0666... 0666... | 0°666... | 0°666... | 0°666... 
05442 ‘69857 77482 ‘85107 "92732 
64343 "70182 80265 "90348 "00431 
63314 "70124 81883 93641 "05400 
62336 "69893 "82943 95992 "09042 
*58707 68235 "84688 ‘O1I4I "17594 
"55385 66148 "84734 "03320 ‘21907 
52254 63913 "84045 ‘04177 24310 
“49261 61617 ‘82954 ‘04291 *25628 
*46371 "59295 ‘$1611 "03928 ‘26244 
"43566 *56966 ‘SO105 "03244 *26382 
*40830 "54639 ‘78481 "02325 °26168 
“35151 *52315 °76772 ‘012209 *25687 

| 0°35522 0*50000 | 0°75000 ‘00000 ‘25000 


‘666... 
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‘I7159 
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34947 
40493 
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“50144 
*50000 


666... 
04305 
"15735 
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28848 
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"50116 
54866 
‘58012 
60116 
‘61502 
"62357 
“62808 
‘62944 
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Fic. 1.—The function f(A, B). The full curves represent the function for 
different values of B. The maxima lie in the dotted curve. 
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The equation (2.11) can be expressed in the alternative form (1) 


cl ‘ ‘ 
=AN2 4 4(1 —d) | oe) dy’, (4.2) 
0 


ne 
H(u) 
and putting A=0, we have 
I ‘s b(n’) , 
=$| —— dy’. 4-3) 
Alu) ~ J, B+p 
An iteration starting from H,(y) as a first approximation suggests that, when A 
is small, H(u) can be expressed in the form 
H(p) = Ho(u) — A"? Ay(u) +AH,(u) —.--, 
but it is doubtful whether the sequence of iterations converges. If, however, 
we substitute for H(j) in (4.2) 
H(p) = Ho(p) — A?) + Ah, {), (4-4) 
where H,(y) is the solution of the integral equation 
€ . ’H : , 
Hy) Hou) {x1 [ Dg, (4. 
o R+pP 
it can be shown that /,() is bounded* as A-~ 0, so that the last term on the right- 
hand side of (4.4) is O(A). 
Let «,™, a be the first and second moments of H,(u). Then multiplying 
(4.4) by » and »? in turn and integrating over (0, I), we get 
Oy = 4 — AN2g 4 OA) =2/ 4/3 —AN®a, + O(A), (4.6) 
dy = ay — ANFa9 + OA), (4-7) 
and substituting these into (2.12) we have the identity 
Ml — dar” + (82) + (1 —A){H — A824) 4/3 + O(A)} =F. 
In this the coefficient of A"* must vanish and sot 
a) = 4/30. = 1-420 892 (4.8) 
by (4.1). 
Differentiating (4.6), we have 
da, /dX = — ha,Pr 424+ O(1). (4.9) 
From (2.22) 
d > da . da 
SMA, B) = — M1 Any? + (1 —A)%0y St — + LBA, + MASI (4.10) 
and by (4.6) and (4.9) 


d I 


57 Bm" + O(1). 





* From (4.4) it is only known that A/j() is bounded as A > 0 foro <j <1. From this and from 
the integral equation which /:;(4) is found to satisfy, it can be shown that /4() is itself bounded as 
Ao, 

+ An alternative way of establishing (4.8) which provides an independent check of (4.4) is as 
follows. By means of (4.3), (4.5) can be written 

er 
Hp) =3 [ —— (Hy (wH(n’) Hp’) H(p)} dp’. 
Joh re 
Multiplying by »? and integrating over (0, 1), we have 
Pd ty’ 
as 4 | | HE Hu )Hi(n!) + Hole’ Hu) de dp’. 
JoJonTrp 
In this » and ,»’ can be interchanged and the two equations added. This gives 
22,0 ${2, 0% a(x (0) a1 (er (0) 20,0), V3 


by (4.1). 
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It follows that the graph of f(A, 8) has a finite gradient at A=o only if 

B = $a, = 0-710 446, 
If 8 >0-710 446, f(A, 8) has an infinite positive gradient at A=o; if 8B <o-710 446 
it has an infinite negative gradient. Thus, provided that 8,>o-710 446, it is 
possible to satisfy (2.23) for sufficiently small A, and A... 

From the graphs in Fig. 1 it is seen that, when 8 >0-710 446, f(A, 8) has a 
single maximum for 0<A<r. As 8 increases, the abscissa of this maximum 
also increases and it is important to determine the value of 8 for which it is exactly 
atA=1. For all greater values of £, f(A, 8) is steadily increasing in (0, 1) and the 
condition (2.23) for emission lines (with 8,= ) will always be satisfied since 

A. 

From (4.10), with A=1, we have 


| 6] > 1+A| 4 A Neg + AN? 2 dey Z|... 


When A= 1, «, =} (see (2)), a so 


[a | =-3+f E + (3)... |- (4.11) 


Differentiating (2.11) with respect to A and putting A=1 we have, since H({j) =1 
when A=I, 


[ dy 5 , : 1 
E Hw)| —4 bro —tuln(1+p-!). 


Hence 


(FR)... -| : Hu) | pedy 


-3f pe? ln (1 + po!) de 
J0 


ae 1 
=5— $In2. 


Substituting into (4.11), we have 


d 
[qe] = —4+/(4—4In2), 


and this is zero when 
B=1/(1—-In2)=3: 
Thus for 8,>3°258889 the condition (2 
satisfied ; if 


258 889. 
23 


3) for emission will always be 


0°710 446 <B, <3°258 889, (4.12) 
the condition will be satisfied if A, and A, both lie to the left of the maximum, 
and it is possible for it to be satisfied if they are separated by the maximum. ‘The 
maxima of f(A, f) as f varies lie on the dotted curve in Fig. 1. 

5. Some numerical examples.—(i) If Ty=5 00odeg., then for visible light 
x=hv/kT, varies between about 3-6 and 7:2; hence (I —e~")/x varies between 
0-271 and 0-139. Taking a=% in (3.2) we see that f, varies between 
0°7230 |(K +a) (corresponding to the red end of the spectrum) and 0-3670/(K +o) 
(corresponding to the violet end). Even the largest value of £, will not fall within 
the range given by (4.12) unless o/(K +o) >0-9 826 making the ratio o : x greater 
than 70: I1—an impossible value because there are very few free electrons at this 


4* 
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temperature to produce the continuous scattering. ‘Thus emission lines will 
not be formed by pure absorption. 

Without investigating (2.17) in detail when €<1, it is possible to see that 
emission lines cannot be formed in such a star partly by absorption and partly 
by scattering. Sinceg is small, it follows from (2.14) that A, is very nearly unity. 
But we have shown that the formation of emission lines is impossible unless 
>A, and so <1. This brings us back to the case which we have already 
dismissed as impossible. ‘There is strong evidence, however, that scattering 
plays an important part in line formation in such a star, so that € is, in fact, small, 
and certainly €<A,. ‘Thus lines formed in this way must be dark, and an explana- 
tion of any bright lines which do occur must be based on some other hypothesis 
(see Section 6). 

(ii) If 7,= 30000 deg., then for visible light x varies between 0-6 and I:-2, 
and (1—e~*)/x varies between 0°752 and 0-582. Again taking a=, we see 
that 8, varies between 2-0050,(K +a) and 1°5520,(kK+a). ‘These values of /, 
all lie within the range (4.12) if 0 /(K +0) >0°457 8, Le. if o/K >0°844 3. 

(iti) Underhill has recently published details of a model for an early O-type 
star (4), which she has classified as Og:5. For this 7,=30000deg. The 
coefficients* o and k vary with depth, but the average value of a /(K +o) over the 
range 0 <7 <5 is 0°327. With this value, together with the results of (ii), we 
find that the greatest value of f, for visible light is 0-655 6; hence emission lines 
are not likely to occur. 

(iv) It is unlikely that the upper value 3-258 88g will be exceeded by the 
B, of any star. If, for example, 7, = 40000 deg., the largest value of (I —e~*)/x 
for visible light is about 0-806 (x =0-448). With this value, £, > 3-258 889 only if 


: 2 ,, 3:258 889 
ak+a™ 4x0806 


<T1-0!, 


an inequality which is impossible unless a<1 and o is very large compared 
with &. 

6. Discussion._From the examples just considered we can see that conditions 
favourable for the formation of emission lines are to be expected in a star of 
high surface temperature with an atmosphere in which electron scattering is of 
greater importance than continuous absorption. ‘The conditions will, moreover, 
be more favourable in the longer wave-lengths, so that emission lines are likely 
to appear first in the red. 

Unsdld, in his investigations of the B-type star t Scorpii (§), concluded from 
the observed central intensities of the absorpzion lines that the lines are formed 
by pure absorption. ‘The fact that emission lines in a high-temperature star can 
be explained theoretically on the hypothesis of pure absorption but not on that 
of pure scattering gives independent support for Unsdld’s inference. If this 
hypothesis is indeed correct, that of detached shells, which has been evoked as 
an explanation of emission lines in a hot star, becomes unnecessary. 

In low-temperature stars the formation of emission lines is impossible 
according to the theory developed here. ‘The theory depends upon assumptions, 
such as the independence of depth of the coefficients o and x,, which are known 
to be untrue, but it should provide a first approximation to the stellar reality. 
‘The figures given in Section 5(i) are sufficient to show that no refinements of 


* Underhill uses the Chandrasekhar—mean absorption coefficient for x. 
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the theory are likely to account for bright lines in the spectrum of a star whose 
surface temperature is 5 000 deg. or less. It follows that we must look for an 
entirely different hypothesis to account for the emission lines in the spectra of 
some long-period variables of low surface temperature. This must take into 
account the peculiar conditions obtaining near the surface of the variable at 
times of maximum and minimum brightness. 

7. Mixed absorption and scattering.—In Section 4 we have only considered 
lines formed by pure absorption. We can, however, show that, if a line is found 
to be dark on a hypothesis of pure absorption, it will be dark if it is formed partly 
by absorption and partly by scattering. In other words, no hypothesis inter- 
mediate between pure absorption and pure scattering will account for the forma- 
tion of bright lines, if they cannot be explained on the hypothesis of pure absorption. 

We have seen that a line is dark if € <A, and we must therefore consider values 
of € for which A, <€ <1, or values of A,“ for which 

A. A. As (7.1) 
where A,“ and A, are given by (2.2) and (2.24) respectively. 

From (2.7), (2.14) and (2.24) we have 

n,=(1—A,)/(I—A,) 
and (2.15) can therefore be written 
p= EARL?) + BRA) 
" (L=A_)Ri(Ae) + BeRa(Ac) 
where 8, is given by (2.25) or by (3.2). From (2.21), 
JO, By) — Qu ARAL) _ fs, Be) aa 
; fre Br) Fre Be) 

We have now to distinguish two cases: £,<0-710446 and 8, >0-710 446. 
In Section 4 we have shown that /(A,f,) is steadily decreasing in (0,1) when 
f£,<0-710446. In that case 

AAS”, B,) < fs Ba), 
since A,“ >A,, and it follows from (7.3) that R,<1. ‘Thus when £, <0-710 446 
lines are necessarily dark whether formed by pure absorption or by part absorption 
and part scattering. 

Now suppose that 8, >0-710446 and consider the numerator N(A,) of 
(7.2), where 


(7.2) 








N(A) = (1 —A,)Rj(A) + BRa(A). (7-4) 
We shall show that (A) is an increasing function of Ain (0,1). From this it wil! 
follow that, when A,“ <A,, 
N(A,) <N(A,) =f(Aw By) 
and hence that 
flr, B,) 
R,< * 
Src Br) 
If a line is found to be dark on the hypothesis of pure absorption, then 
SA, B,) ; flrcr B,) 
and it will follow that R,<1 and that the line is dark on a hypothesis of part 
absorption and part scattering. 
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In order to establish that N(A) is an increasing function of A in (0,1), we 
write (7.4) in the form 
N(A) =A(A) —A,R,(A), (7-5) 


where 


h(A) = R,(A) + B,R3(A). (7.6) 
Since R,(A) is steadily decreasing in (0,1), the result will follow if h(A) is steadily 
increasing. Now from (2.12) and (2.13) we have 
R,(A) = 3 —AM? ag, 
and so, by (2.10), 
(A) = § + AY Ba — a9) 
= § + RAA)(B, — 9/1). (7-7) 
From the tables for «, and a, given in (2) it is found that «,;«, decreases steadily 
from the value 0-710 446 when A=o to the value § when A=1. It follows that 
6, —&q/a, is positive and steadily increasing when 8, >0-710446. Since the same 


is true of R,(A), it follows from (7.7) that A(A) is also steadily increasing in (0, 1). 
This completes the proof. 


In conclusion I should like to thank Professor H. H. Plaskett for drawing 
my attention to Unsdld’s paper on + Scorpii and for some valuable criticisms. 


St. Hugh’s College, 
Oxford : 
1951 October 20. 
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Summary 

An account is given of observations made on quiescent prominences in the 
Ha line with the combined spectrograph and spectrohelioscope at the Royal 
Observatory, Edinburgh. Prominences were photographed in emission on 
the limb and in absorption on the disk. Corrections were determined and 
applied for the effect of grating ghosts, scattered light and limited resolving 
power. The average corrected intensity for prominences observed in 
emission on the limb was 10-4 per cent of the continuous spectrum of the 
centre of the Sun’s disk 13 A outside Ha. The half-width of the Ha 
emission line ranges between 0°64 and 1-30 A. 





1. Introduction.—Line profiles of Ha in quiescent prominences have been 
measured previously by Woolley and Newton (1) at Greenwich and by Unsdéld 
(2) at Mount Wilson. The measurements of Woolley and Newton were made 
visually with a spectrohelioscope. Intensities were found by means of a wedge 
photometer, prominence intensities being compared with the brightness of a 
lamp, which in turn was compared with the Sun’s continuous spectrum 15 A 
outside Ha at the centre of the disk. Observed profiles were corrected for the 
finite width of the spectrohelioscope slits, and for the effect of heterochromatic 
light in the instrument. 

Unséld measured profiles photographically but applied no corrections. 
Only three plates were obtained, and two of these show a considerable background 
of scattered light. 

The present work was undertaken with a view to determining line profiles 
to a greater degree of accuracy than has been attained by either of the above 
methods. ‘The improvements may be summed up as follows. 

(a) By using the photographic method, only one slit was used. In Woolley 
and Newton’s visual method, corrections had to be made for the slits of the 
spectrohelioscope. Further, their correction was for finite slit-width only; 
thus they assumed that the optical performance of the instrument was limited 
only by the slit. In this present work the instrumental profile was measured 
directly, so that the errors due to the imperfect performance of the instrument, 
as well as those due to the finite width of the slit, have been accounted for. 

(6) Scattered light was measured directly. 

(c) A correction was made for grating ghosts. 

(d) In the photographic method, the spectrograph slit could be set to gather 
light from any desired part of the prominence, therefore intensity and profile 
measurements made on such a plate represent the conditions at a particular region 
in the prominence. ‘The intensity of a prominence is by no means uniform over 
its whole area, though it may appear uniformly bright when observed in the 
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spectrohelioscope. By Woolley and Newton’s visual method the whole pro- 
minence was observed at one time. 

(e) When making exposures on limb prominences, an image rotator was used 
to orient the Sun’s image until the limb near the prominence was brought 
parallel to the slit. In this way all photospheric light was excluded from the 
spectrograph during prominence exposures, thus avoiding a background of 
scattered light which may be as great as 60 per cent of the maximum brightness 
of the prominence, as was the case in one of Unséld’s prominences. 

2. Procedure.—The observations described here were carried out with the 
combined spectrograph and spectrohelioscope at the Royal Observatory, Edin- 
burgh. ‘This combined instrument was constructed by Dr M. A. Ellison and 
is designed for the photography of the spectra of flares and other transient 
phenomena. It is shown diagrammatically in Fig. 1. ‘The Sun’s light is directed 
by a siderostat through a 5-inch achromatic lens of 18 feet focal length, which 
forms an image of the Sun 2 inches in diameter on the primary slit 5,._ ‘The 
instrument, which is autocollimating, uses a 54-inch Rowland plane grating G of 
80 000 lines. ‘The achromatic lens L, which has a focal length of 16 feet, fills 
the grating with parallel light and also serves to focus the spectrum either on the 
viewing slit S, or on the plate in the camera at P. The Sun’s image may be 


| 


| 
nen be 


Fic. 1.—Combined spectrohelioscope and spectrograph. 


viewed in Ha light in the spectrohelioscope. For photography, the grating 
may be rotated until the required part of the spectrum is brought on to the plate; 
this is easily done by means of a handle at the eye end. ‘The primary slit is 0-1 mm 
wide, a suitable width for spectrohelioscope work, which is the main function of 
this instrument. ‘This width is, however, greater than the optimum value for 
a spectrograph of this type, which from simple theory (3) is0-o2mm. The slit 
jaws were adjusted to parallelism by closing them against a steel shim, previously 
tested for uniformity of thickness with a micrometer. ‘The jaws were clamped 
in this position and have not been, altered throughout the course of the work 
presented in this paper. 

Each prominence to be photographed was first located visually in the spectro- 
helioscope, and its dimensions and characteristics noted and entered on a drawing 
of the Sun’s disk. During exposure an image rotator was mounted directly in 
front of the slit on an optical bench (B in Fig. 1). ‘This image rotator (constructed 
by Dr Ellison) uses a right-angled prism with sides of 2 inches mounted with its 
hypotenuse face parallel to the axis of the incoming beam. A weak positive 
lens is fixed in front of the prism te compensate for the change in focus. During 
exposures the image was accurately guided on the slit by tilting a plane parallel 
plate of glass in front of the slit. For each exposure the height above the limb 
could be estimated to the nearest ten thousand kilometres by means of two pairs 
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of lines ruled on the faces of the slit jaws at distances corresponding to distances 
of thirty and sixty thousand kilometres on the Sun. For small prominences 
exposures at one height only were made; for large prominences a series of 
exposures was taken at different heights above the limb, ranging from ten to 
seventy thousand kilometres. 

In the case of prominences photographed on the disk, it was not possible 
to eliminate scattered light, since the slit was fully illuminated. ‘There was, 
therefore, no point in using the image rotator for disk prominences. 

Exposures were made in both first and second orders. Exposure times varied 
from 5 to 30 seconds, and were in the majority of cases about 20 seconds. ‘The 
camera shutter was operated by hand, timing being done with an electric clock 
with centre seconds hand. 

All prominence spectra were standardized by photographing, on another 
part of the same plate and through a Hilger step wedge, the spectrum of the 
centre of the Sun’s disk. ‘The wedge has six steps and a clear glass portion, 
the transmissions of successive steps being in the approximate ratio of 2:1. 
The transmissions were accurately measured with a rectifier photocell in the red 
end of the spectrum. ‘The intensity at any wave-length in the profile of a pro- 
minence could thus be found in terms of the intensity of the Ha continuum at 
the centre of the Sun’s disk. ‘This intensity was considerably higher than that 
of the prominence, and therefore it was necessary to reduce the light from the disk 
during standardization by means of neutral screens, in order that the exposure 
times on prominence and disk be the same. ‘The transmissions of the screens 
were measured (in the red end of the spectrum) with the rectifier photocell. 
Plates were standardized as far as possible immediately after exposure, to avoid 
any error due to changing sky conditions. Where the image rotator was used 
for an exposure, it was left in position during standardization. ‘The plates used 
throughout were Kodak type III I; they were always developed in M.Q. 
developer at 18 deg. C for five minutes. 

‘The majority of plates were measured with the Kipp and Zonen recording 
microphotometer which was brought into use in the Observatory at the beginning 
of 1949. A few of the plates obtained before this date were measured with the 
non-recording microphotometer of the Observatory. ‘The recording microphoto- 
meter gives a magnification of 50 between plate and bromide tracing. Distance 
lines impressed on the tracings every 5 mm correspond to distances of 0-I mm on 
the plate. ‘The full deflection between clear glass and infinite density was kept 
as near as possible to 10cm. ‘It was found that, with the particular dispersion 
and resolution of the spectrograph, a width of up to 20 was permissible for the 
analysing slit. It was thus possible to adjust the full-scale deflection to locm 
by varying the slit width between 16 and 20. ‘I'he width was never allowed to 
exceed 20;:, so in the case of very short prominence images a full-scale deflection 
was not possible. ‘the illuminating slit was always set so that the width of the 
exploring beam was double that of the analysing slit. 

The standardization spectrum was measured in two regions free from 
Fraunhofer lines and equidistant from Ha. ‘The wave-lengths chosen were 
A0550 and A6576. A short run on each step of the standardization spectrum at 
both these wave-lengths was made. ‘Two clear-glass records were also made, 
one near the prominence spectrum and one near the standardization spectrum 
and near the centre of the plate to avoid edge effects. 
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Fig. 2 (Plate 1) shows a specimen of microphotometer tracings for a 
prominence Hz profile and its standardization. 

Deflections on the tracings were converted to A-values, where the function 
A, known as the blackening, is log,, {(1—7)/7}, T being the transparency, or 
fraction of light transmitted. For the standardization spectra, the A-values 
Corresponding to the wave-length of Hx were taken to be the mean of the 
A-values measured at A6550 and A6576. A-values were plotted against wedge 
“‘ magnitudes’’ (magnitude m= 2-5 log,,(1/t), where ¢ is the transmission of a 
step of the wedge), to give calibration curves which are approximately linear 
and of about unit slope. ‘This procedure was adopted on the suggestion of 
Dr E. A. Baker, who first introduced the A-function into spectrophotometry (4). 

3. Correction for grating ghosts.—The grating used in the Edinburgh spectro- 
graph had been ruled by Schneider on Professor Rowland’s engine and therefore 
it was expected to be comparatively free from ruling errors. The ghosts were 
photographed in both orders for the green line of mercury A5461. As the first 
ghosts on either side of the line seemed to have an intensity of less than I per cent 
of the main line, a very narrow strip of neutral glass of measured transmission about 
I per cent was mounted directly in front of the plate so as to intercept and transmit 
the light from the main line. ‘This had the effect of reducing the intensity of 
the main line so as to be comparable with that of its ghosts. The ghosts were 
found to be symmetrically placed about the main line at distances from it of 
nd, where n( = 1, 2, 3, etc.) is the order of the ghost, d= 7A for the first order and 
d= 3°5A for the second order. ‘They were somewhat brighter on the red side 
of the line than on the violet side. Exposures of 15 to 30 minutes were given on 
the ghosts and the plates were standardized for the same length of time, using 
light from the centre of the Sun’s disk. It was necessary to reduce the intensity 
of the Sun’s light by stopping down the aperture of the objective to I or 2 inches. 
Four plates were obtained in each order. One of the first-order plates was 
exposed without the I per cent screen; this plate gave the relative intensities of the 
ghosts but not their absolute values. ‘The plates show a background of scattered 
light, which is termed linear scattering and which will be discussed later. ‘This 
background of scattered light had to be subtracted from the apparent ghost 
intensities to give true values. It is fortunate that the ghosts are weak, so that 
the errors in the estimation of their intensity due to the presence of scattered light 
are without importance. 

‘Though seven or eight ghosts are visible on either side of the line, only four 
were measurable on the plates. ‘The intensities of these four are given in 
Table I. The totals given include contributions of the fainter ghosts which 
were found by extrapolation. 

TABLE I 
Intensities of Grating Ghosts in Percentages of the Main Line 
Red Violet 

Ghost 1 2 3 4 I 2 3 4 Total 

First Order 0°20 0°07 0°04 0°03 ©'10 0°05 0°03 0°02 0°60 

Second Order 0°62 0°23 O'14 0°09 0°52 0°26 0°14 0°09 2°34 

The ratio of the total ghost intensities in the first and second orders is in good 
agreement with Minnaert’s theory (§) which predicts a ratio of 1: 4. 

‘The intensities of the ghosts in the Edinburgh spectrograph may be compared 
with those of the Mount Wilson 150-foot tower telescope which is equipped with 
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a Michelson plane grating. The total ghost intensities for this spectrograph, as 
measured by Minnaert, Mulders and Houtgast (6), are 0-61 per cent in the first and 
2-34 per cent in the second order, though the relative intensities of the individual 
ghosts are different from ours. For the Utrecht spectrograph, which is similar 
to the Edinburgh instrument and which has an 8-cm plane Rowland grating of 
45 000 lines, the total ghost intensities found by Houtgast (7) are 1-25 per cent 
in the first and 5-1 per cent in the second order. 

If the total ghost intensity is g, that is, if each part of the spectrum loses a 
fraction g of its light to its ghosts, then the true intensity 7 of an emission line is 
found from the observed intensity 7’ from the relation 


i=1'/(1—g). (1) 
The value for the ghost intensities derived in Table I are the intensities of the 
ghosts compared with the residual line. The last column of the table, therefore, 
does not represent g but g/(I—g). The difference between the two is of the 
order g, and in practice it is sufficient to take g to be 0-006 in the first order, and 
0023 in the second order. 

The correction of absorption lines for ghosts depends on how the intensities 
have been measured in the first instance. Normally, the intensity of an absorption 
line is measured as a fraction of the continuous spectrum outside its own wings, 
and when this is the case the correction is the usual one (6,7), 


r=(r' —g)/(I—8), (2) 
where ¢ is the true intensity of the line and r’ is the observed intensity, the intensity 
of the continuous spectrum on this scale being 1. Intensities measured in terms 
of the local continuum in this way will be referred to in this paper as “ relative”’ 
intensities. 

In measuring the intensity of a prominence which is not located on the centre 
of the Sun’s disk, the intensity may also be given in terms of the continuous 
spectrum for the centre of the Sun’s disk as unit. ‘The value of this intensity 
is found directly from the standardization spectrum, and will be referred to as 
the ‘‘absolute”’ intensity. In this case the intensity of the continuous spectrum 
at the centre of the Sun’s disk is put equal to 1, and the intensity of the continuous 
spectrum at the position of the prominence is now not I but o, the limb-darkening 
coefficient. ‘The quantity o was known for each prominence, as will be described 
in Section g. ‘The centre of the absorption line whose true intensity is 7 loses 
therefore an amount of light gi and gains an amount go from the local continuous 
spectrum. Hence the observed intensity 7’ is given by 7’ =(¢—gi) + og or 


t= (#' — og)/(1 —g). (3) 
4. Correction for heterochromatic light.—\t is well known that one of the main 
sources of error in spectrophotometry arises from heterochromatic light being 
scattered from the walls and surfaces of the instrument on to the plate. General 
precautions were taken to reduce this scattering to a minimum. 
(a) ‘The walls of the spectrograph were covered throughout with black velvet. 
(6) A stop O, (Fig. 1) was placed at the etfective focus of the outer concave 
surface of the spectrograph lens to eliminate the real image of the first slit formed 
by reflection from this surface. This stop (a black rectangle 10mm = 5 mm) is 
placed 30 inches in front of the lens on the axis of the spectrograph. Another 
stop O, in the form of a horizontal strip 6mm wide, mounted at the level of the 
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centre of the lens, cuts out the light from the virtual images of the first slit formed 
behind the lens by reflection from the other lens surfaces. 

(c) A red filter (Schott and Gen. R.G. 2) which transmits only wave-lengths 
greater than A6200 was mounted in front of the slit. Since the amount of scattered 
light depends on the total light entering the instrument, this had the effect of 
reducing scattered light as well as preventing overlapping of orders. 

The intensity of the scattered light reaching the plate was determined by 
experiment. A metal strip a few millimetres wide was fixed horizontally across 
the centre of the (vertical) slit of the spectrograph. When an exposure was made 
on the centre of the Sun’s disk, it was found that a certain amount of light had been 
scattered into the unexposed band across the centre of the spectrum corresponding 
to the part of the slit which had been occulted. Exposures of ten times the 
normal exposure time had to be given in order to show up the scattered light, and 
the plates were standardized in the usual way, using the centre of the Sun’s disk. 

On examining the occulted region on the plates it was found that Fraunhofer 
lines could be traced to distances of 1-2mm from either edge. This effect is 
astigmatism, which causes the spectrum to be spread in a direction perpendicular 
to the dispersion. Beyond 12mm the scattered light was uniform both with 
distance from the edges and along the dispersion. ‘This indicates that the 
scattered light which we want to measure is spread uniformly over the spectrum. 
In order to measure the intensity of this uniformly scattered light, it is obvious 
that the minimum width of occulting strip must exceed 2-4mm. In these 
experiments two strips were used of widths 3mm and 65mm. ‘The exact 
width of the occulted region was measured on each plate with a travelling 
microscope. A series of plates was taken in both orders, the total length of slit 
being varied by means of a pair of horizontal jaws mounted in front of it. 

The results found are plotted in Fig. 3. Each point on the curves represents 
the mean of two plates taken under identical conditions. ‘The length of slit 
illuminated is the total length minus the width of the occulting strip, and the 
intensities are measured in terms of the continuous spectrum outside Hx at the 
centre of the disk. Scattered light was assumed by Shane (8) for the Mount 
Wilson spectrograph to be proportional to the length of slit illuminated, but the 
curves in Fig. 3 show that this is only approximately true for our spectrograph. 

A somewhat similar method for the determination of scattered light was used 
by Shane (8) in his study of absorption lines made with the Mount Wilson spectro- 
graph. In Shane’s experiments over-exposed plates were taken of the Sun’s 
spectrum and the light scattered beyond the edges of the spectrum was measured. 
This light is made up of two components—general scattered light which is the 
quantity required, and light scattered by astigmatism in a direction perpendicular 
to the spectrum. Perpendicular scattering was measured separately from an 
over-exposed plate of a single emission line, which was, of course, free from 
general scattering. ‘lhe difference between the total scattering and the per- 
pendicular scattering gave the general scattering. ‘Ihe advantage of the occulting 
strip method used here is that the general scattering could be measured directly, 
independently of astigmatism. Moreover, the scattering was measured at the 
centre of the tield, just where we require to know its value, and not beyond the 
edge of the spectrum, where its etlect may be appreciably different. 

A certain amount of inaccuracy arises in estimating the scattered light correction 
to be applied when the full slit is illuminated. ‘The total slit length is 26mm, 
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and, as has been explained before, the minimum occulting strip which could be 
used in our experiments had to exceed 2-4mm in width. Hence the last point 
on the curves in Fig. 3 corresponds to only 23 mm length of illuminated slit, so 
that the scattered light for the total slit must be found by extrapolation. ‘lhe 
obvious way to avoid this is to use a shorter slit for exposures on the Sun’s disk. 
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Since the image rotator was used in photographing prominences on the limb, 
and no direct photospheric light was allowed to enter the instrument, there is for 
such emission prominences no appreciable heterochromatic light scattered on to 
the plate. ‘The standardization spectrum, however, is enhanced by hetero- 
chromatic light, and the intensit‘es measured from it must be corrected accord- 
ingly. ‘The total effective length of slit illuminated during standardization 
is 94mm. _ If the intensity of the scattered light corresponding to a slit length 
of 9'4mm is « times the intensity of the Hx continuum, then the apparent 
intensities found from the standardization must be increased by « times the 
intensity of the H« continuum for the light incident on the slit. ‘The intensities of 
scattered light plotted in Fig. 3 were found, using the uncorrected standardization 
(as given by the photoelectrically determined transmissions of the wedge steps’, 
and so they need correction, with the result that the value of « is not that corre- 
sponding to 9-4 mm on the appropriate graph in Fig. 3. In the first order, exposures 
for scattered light were standardized through a neutral screen of transmission 
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0°25, so that all intensities on the first-order graph are too small by an amount 
0:25a. From the curve we find the apparent scattering at 9-4 mm to be 0:56 per 
cent, giving « =0°56+0-25«, so that «=0-75 percent. Similarly, in the second 
order the screen used for standardization was of transmission 0-364. The apparent 
intensity of the scattered light at 9-4 mm is I-20 per cent, giving « = I-20 + 0-364a, 
or « =1-89 per cent. 

In general, then, where a screen of transmission 7' is used to reduce the 
intensity of the standardization spectrum, the corrected intensity / is found from 
the apparent intensity J’ by using the formula J/=/'+«T7, the unit of intensity 
in this formula being the intensity of the Ha continuum at the centre of the disk. 
This correction affects every point on the line profile and is proportionally 
greatest at the low-intensity values. ‘The best method of procedure would have 
been to adopt new effective transmission values for the steps of the calibrating 
wedge and thus avoid completely the necessity of correcting each intensity 
measurement made, but as most of the reductions of plates were done before the 
scattered light was accurately measured, this was not practicable. ‘The intensity 
measurements set down as “ observed ” intensities in ‘l'ables III and V incorporate 
this standardization correction. 

When measured intensities in an absorption line have been corrected for the 
standardization effect described above, a further correction is necessary for the 
light scattered over the spectrum itself. In photographing prominences on the 
disk, the whole slit was always fully illuminated. If the scattered light corre- 
sponding to the total slit is 87 (where J denotes the intensity of the continuum), 
then each part of the spectrum is enhanced by an amount of light 8/._ Therefore, 
if r is the true relative (see Section 3) intensity, the observed relative intensity 
r’ is given by 7’ =(r+)/(1+8) or 


r=r'(1+B)—f. (4) 


The absolute (see Section 3) intensity is affected in a different way. If 7 is the 
true absolute intensity, 7’ the observed intensity and o the intensity of the local 
continuum, then 

i=i' — fo. (: 


The value of f is 2-8 per cent in the first and 7-7 per cent in the second order. 
These values were found by reading off the graph in Fig. 3 the percentage 
scattered light corresponding to 26mm, which is the total length of the slit, 
and adding the correction «7’. 

5. Discussion of other forms of scattered light 

(a) Astigmatism.—Light scattered in a direction perpendicular to the disper- 
sion by astigmatism does not affect profile measurements at all and affects absolute 
intensity measurements only if the object photographed is very small. It has 
been stated above that the effect of astigmatism can be traced out to 1-2 mm from 
the edge of the spectrum. _ Its actual intensity has been measured from the plates 
used for measuring scattered light, and is found to fall off rapidly with distance 
from the edge of the spectrum. At a distance of 0-3 mm from the spectrum the 
effect of perpendicular scattering (found by subtracting the known general 
scattering from the total scattering) is only 0-5 per cent. As the length of the Hx 
line of a prominence is usually a few millimetres on the plate, there will be no 
change of the density of the central part of the image due to this effect. 
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(6) Linear scattering —The light of a particular wave-length scattered by 
the grating along the dispersion into other wave-lengths is termed linear scattering. 
It arises from random irregularities in the ruling of the grating in the same way 
as periodic ruling errors give rise to ghosts. Its magnitude was estimated by 
Shane (8) from highly exposed plates of the mercury green line used for measuring 
ghosts. The effect is apparent on our ghost plates as well, but, as ‘Thackeray (9) 
has pointed out, it is impossible to distinguish true linear scattering from the 
extended wings of a broad line such as the mercury line. No attempt was made 
to measure linear scattering in this way, because since it falls off with the distance 
from the centre of the line, it is included in the wings of the instrumental profile. 
For this reason it was not necessary to investigate linear scattering separately. 

(c) Scattered light due to indirect photospheric light.—\t has been pointed out 
already that no direct photospheric light was allowed to enter the spectrograph 
during exposures on limb prominences. Photospheric light enters the slit 
indirectly, however, as a result of diffusion in the atmosphere and of scattering by 
the telescope lens and the siderostat mirror. Atmospheric scattering has been 
measured by Fracastoro (10) and found to be 10° * times the intensity of the Sun’s 
photosphere at 2’ beyond the limb under normal seeing conditions at a level of 
2700 metres. ‘The same figure has been obtained by Kienle (11) for the minimum 
sky brightness at a distance of 0-1 R fromthe Sun’slimb. Under normal conditions 
he estimated the intensity to be at least ten times this value, or equal to 10 ? times 
the photospheric brightness. ‘Though this effect is too small to show up on the 
photographic plate in the exposure time given to the prominences, it means, 
nevertheless, that the intensity values measured for prominences are all 
systematically too great by an amount which may be up to I percent. ‘The amount 
varies with the sky conditions and with distance from the limb, so that ideally 
it should be measured for each individual exposure. 

Scattering by the lens and mirror can be estimated from measurements made 
by Lyot (12) on an 8-inch lens of focal length 12 feet and of medium quality. 
He found the scattering in this lens to be of the order of 10 * at 2’ beyond the 
Sun’s limb, a distance which corresponds approximately to 100000 kilometres 
on the Sun. The object glass which forms the Sun’s image at Edinburgh has 
a very perfect polish and is free from bubbles, so that the scattered light intensity 
probably does not exceed Lyot’s figure. A similar source of scattering is the 
siderostat mirror, whose aluminized surface was kept in good condition so that 
it probably scatters light to the same small degree as the lens. ‘The total effect, 
therefore, appears to be small, and as it decreases rapidly with height above the 
limb, it will be negligible at distances of, say, 20000 kilometres upwards. 

6. Correction for limited resolving power.—The instrumental profile was 
measured for the krypton green line 45571, which has a natural half-width of 
o-o1 A (13). ‘Three plates were obtained in each order; exposure times were 
30 minutes in the first order and 2 to 2} hours in the second order. ‘These 
exposures were too long to allow the plates to be standardized in the usual way, 
using the centre of the Sun’s disk; they were therefore standardized on the 
multiple-slit calibration spectrograph of the Observatory. ‘The wings of the 
profile (intensities below 7 per cent of the centre of the line) were found from over- 
exposed plates, one in each order, which gave the intensity down to I percent. ‘The 
exposures which had to be given were 4 hours in the first and 7 hours in the second 
order. Exposures over a longer period, to show up the wings farther out, were 
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not made as it was considered that such long exposures would result in loss of 
accuracy due to temperature changes in the spectrograph, or shifts in the plate. 
The line profiles reproduced in Fig. 4 show that most of the broadening is due to 
the finite width of the slit (the slit width is marked on the figure). The profiles 
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are not perfectly symmetrical at the wings, but it was sufficient for purposes of 
reduction to take the mean of the two wings for the standard profile. This also 
simplified the working out of corrections. 

The observed profiles were corrected by Burger and van Cittert’s method 
of successive approximations as outlined by van Albada (14). If the true profile 
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is W, the instrumental profile A, and W(A) the profile W broadened by the 
instrument, the problem is to solve the equation W(A)=TJ, J being the observed 


ripe 
profile. A is normalized so that | Adx=1, where x represents the distance 


~ « 
from the centre of the line in arbitrary units. In this case the units chosen were 
o-10A in the first order and 0-05 A in the second order. ‘Table II shows the 
correction to a typical prominence in the second order and will serve as an 
illustration of the method. Column 1 gives distances from the centre of the line 
in angstroms. Since both the prominence and the instrumental profiles are 
taken to be symmetrical, only one-half of the prominence profile has been tabulated. 
The intensities at the various points along the observed prominence profile, 
corresponding to the displacements in column 1, are set down in column 2. 
All intensities have been multiplied by a constant factor which reduces the 
intensity of the central maximum to 1-00. Intensities which were too low to be 
measured on the plate were extrapolated and are given in brackets. 


TABLE II 
Correction of an Emission Line for Instrumental Profile 


I 2 3 4 5 6 

AA I I(A) D,=I1(A)—I Ws. W,(A) 

° 1°00 0'98 0°02 1°02 1°00 

o'oSs 0°99 0°97 o°o2 I‘o!l o’99 

o'10 098 0:96 002 1'00 0°99 

Or15 0°96 0°94 0°02 0°98 0°98 

0°20 0°94 °'90 0°04 0°98 0°94 

0°25 0°89 0°85 0°04 0°93 0°89 

0°30 0°84 0°77 0'07 o'"9! o'$2 

0°35 0°73 0°68 0°05 0°78 0°73 

0°40 0°58 0°57 oro! 0°59 0°58 

0°45 0°43 0°46 + 0°03 0°40 0°44 

0°50 0°30 0°35 o'oOs5 o'25 0°29 

0°55 o'21 0°25 0°04 O'17 o'21 

o'60 o'13 O'l7 0°04 0°09 O'13 

0°65 (0°08) o'll 0°03 0°05 0'07 

0°70 (0°04) 0°07 + 0°03 o'ol 0°04 

The observed profile 7 was taken as the first approximation. Under the 
influence of instrumental broadening this becomes /(A). The difference 
D, =1(A)—J (column 4) is subtracted from J to give the second approximation 
W, (column 5). W, is now broadened to W (A) and the difference D, = W,(.4)—/ 
is subtracted from W, to give the third approximation W,=W,—D,. This 
process is continued until W,, is reached, where W,(A)=J. The calculation of 
the (A) values involves integrating the contributions to intensity made by every 
part of the profile W. ‘To avoid plotting several series of points and measuring 
areas in every case, the integration was carried out by Simpson’s rule, using 
intervals of I unit (o-r A in the first order and 0-05 A in the second order). Asa 
tinal check, areas could be plotted and measured with an Amsler planimeter. 
It was important to verify occasionally during the approximations that the 
areas under the curves / and IV were the same, as the area represents the total 
amount of light emitted and cannot be changed. 
The number of successive approximations which had to be made depended 

on the shape of the profile in question. In the first order, in the case of a profile 


=] 
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0-98 A wide at half-maximum, nine approximations were required. For second- 
order spectra, the instrumental profile correction is much smaller, and two or 
three approximations were usually sufficient. Fig. 5 shows observed and 
corrected profiles for one typical prominence in each order. 

On account of the labour involved in making the instrumental profile correction, 
the observed profiles of limb prominences were not individually corrected. 
Instead, three profiles in each order, which we shall call standard profiles, of 
half-widths (i.e. widths at half maximum) 0-8, 0-9 and 1-0 A, approximately, 
were chosen and corrected. ‘Their intensities were reduced to give a central 
maximum of 1-00, as described in connection with ‘lable II], and their observed 
profiles were plotted on thin paper. All other profiles, similarly reduced to a 
maximum of 1-00, were plotted on the same scale as the standard profiles. ‘T'o 
correct a given profile, the standard profiles in the same order spectrum were 
superimposed on it in turn, and the best fit, as regards the shape of the curves 
near the centre, was picked out. Since the increase in central intensity as a 
result of instrumental correction was known for the standard profiles, the correction 
for the profile in question could be estimated satisfactorily. Correction to the 
half-width could not be estimated in this way, but, as the curves in Fig. 5 show, 
the correction is small, so that the compiete instrumental profile correction was 
carried out only in particular cases. 

7. Results: prominences on the limb.—'Table III gives the results obtained 
from all quiescent prominences photographed on the limb. Prominences are 
detined here as being quiescent when they show no Doppler motions in the line 
of sight when observed in the spectrohelioscope with the line-shifter, and when 
their Hz line profile is symmetrical. All of these belong to Class V of the Pettit 
classification (1§): that is, they are long-lived stable prominences which often 
form part of the large dark filaments observed on the disk, and their appearance 
remains unchanged for at least several hours. ‘The table is divided into two 
parts—observations made from first- and from second-order plates. 

Columns I and 2 of the table give the date of observation and the position 
of the prominence (its latitude and the limb on which it was observed) respectively. 
The prominence heights entered in column 3 represent the height of the slit 
abovethe Sun’slimb. ‘There are apparent heights only ; the true heights measured 
from the base of the chromosphere may be greater than this if the base of the 
prominence is situated either behind the limb or in front of it, so that the heights 
given are minimum values. In some cases two or more measures are entered 
opposite a particular height. ‘These have been made on the same plate, but each 
corresponds to a separate uniform part of the same prominence. ‘The observed 
central intensities are given in column 4. Column 5 gives these figures corrected 
for ghosts, according to equation (1) in Section 3, and column 6 gives the same 
figures corrected for both ghosts and instrumental profile. Column 6 therefore 
gives fully corrected central intensities, all intensities being expressed in 
percentages of the continuous spectrum of the centre of the Sun’s disk, 13A 
outside Ha. The half-widths, uncorrected for the reason already described, 
are entered in column 7. 

8. Discussion of results.—We see from column 6 of ‘Table III that the final 
corrected intensities range from 3 to 28 per cent, the average intensity from all the 
measurements being 10-4 per cent. All quiescent prominences observed under 
favourable conditions were photographed, so that these results are considered to be 
representative of quiescent prominences for the years 1948 and 1949. It is well 
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Tape III 
Central Intensities and Half-widths of Prominences on the Limb 


2 3 4 5 6 
Central Intensity 
Position Height Corrected 
(Limb and in —— for Ghosts Half-width 
latitude) 10° km Observed : and Observed 
for Ghosts 
Instrumental 


Profile 





(a) Observations from first-order plates 


1948 Sept. W. —30° 15 5°3 5°3 5°6 
1948 Oct. E. + 20° 12 8-6 8-7 8-7 
1948 Dec. W. —29° 20 10°9 I1‘o 11°23 
1949 Feb. W. +28° 30 8-7 8-8 9'2 
7°6 7°6 8:0 

1949 Feb. W. + 26° 30 26°6 26°8 (27°0) 
60 9°7 9°8 (9°9) 

1949 Feb. W. +30° 10 17'I 17°2 (17°4) 
20 10°3 10°4 10°6 

30 8-8 8-9 9°2 

40 g'1 9°2 9°6 

50 3°2 3°2 3°2 

60 4°0 40 4°4 


(6) Observations from second-order plates 


1948 Apr. - 54° 5 8:8 9'0 g'1 
1948 May t. 41° 10 15°7 16°1 
30 15°3 15°7 158 
1948 May .— 8° 30 10°6 10°8 
1948 July _> 2 15 4°4 4°5 4°6 
4°3 4°4 4°6 
2°9 3°0 3°0 
5°3 5°4 5°5 
1948 July S.* 14°0 14°5 
17°3 ; 17°9 
1948 Aug. S. +16° 10°8 . 11‘! 
1948 Sept. ‘Son 3°4 3°5 
1948 Oct. .- 7°4 , 7°6 
1949 Apr. ] 7°6 . 8-0 
1949 May I. + 5 5*2 , 5°4 
1949 May 1. + 6° 9°6 , 10°0 
77 “ 8-0 
1949 May 3. 142° 15°6 16°1 
6-1 . 6"4 
1949 May be 2‘8 (2°9) 
81 8°3 (8-4) 
1949 June he 16°8 17‘2 (17°4) 
(19"1) (19°5) (19°7) 
19°4 19°9 19°9 
1949 June & E 81 8-3 8°5 
8:2 8-4 8°5 
10°4 10°6 10°7 
1949 June : 3°7 3°8 3°9 
1949 July = 8-9 g'I 9'2 
10'8 Ii‘! I1‘2 
10°6 10°8 10°9 
1949 July 1, +14° 12°2 12°5 12°6 


re 
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Studies of Hx line profiles in prominences 
TABLE III (cont.) 


4 5 
Central Intensity 
—— Se eee -—-- 
Position Height Corrected 
(Limb and in —! for Ghosts Half-width 
latitude) 10° km Observed . : and Observed 
for Ghosts 
Instrumental 
Profile 
sl 0°76 
8-9 o'88 
12°3 o'82 
16°3 o'89 
15° 0°89 
9°: 0°72 
RS" 096 
nS 0°86 
0°73 
9° o'So 
0°79 
0'70 
0°73 
0°84 
‘SI 
‘96 
12 
‘Ol 
“89 
‘88 
‘00 
‘97 
‘96 
‘93 


) 


1949 Aug. 10 W. —35 
1949 Aug. W. —4! 
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1949 Aug. ‘. 15 
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©0000 #00 
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1950 Mar. 7 E 28 23 23°6 


Note to Table I11,—-Quantities which were found by extrapolation of the standardization 
curve are given in brackets, Where half-widths were not measurable, the instrumental 
profile correction was only approximate; in these cases the intensities in column 6 are also 
given in brackets. 


known that the frequency of low-latitude (less than 40°) prominences follows 
the sunspot cycle; high-latitude prominences, on the other hand, are rare at 
times of sunspot maximum, reaching maximum frequency about three years 
after the sunspot minimum. Although no evidence has yet been brought 
torward on the subjeci, it is possible that prominence intensities are also connected 
with the sunspot cycle. ‘The years 1948 and 1949 were years of very high solar 
activity, their relative sunspot numbers averaging 135 (16), which is higher than 
the relative numbers for the six sunspot maxima preceding 1947. If such a 
correlation exists between prominence intensity and sunspot number, the results 
in ‘Table Il] must not be taken as a general representation of prominence intensity. 
‘They are, however, a true representation for years of high solar activity. 

The frequency distribution of quiescent prominences, classified according to 
their intensity, is given in ‘lable [V. ‘The first column gives the intensity range, 
intensities being given as usual in percentages of the Sun’s continuous spectrum 
outside Ha at the centre of the disk. ‘lhe frequency of occurrence of intensities 
within these ranges is recorded in column 2, and the last column gives similar 
figures found by Woolley and Newton (1). 
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Woolley and Newton’s observations were made in 1935 and 1936, about two 
years before the sunspot maximum. (Relative sunspot numbers for 1935 and 
1936 were 36 and 8o respectively.) 

With regard to the half-width of prominence profiles (given in column 7 of 
Table III), it will be seen that it varies between 0-64 and 1-30 angstroms. In 
another paper it is proposed to make a theoretical analysis of the observed line 
profiles with a view to gaining some information on the temperature conditions 
within a prominence, and on the effect of self-absorption on the H« emission line. 
Further data are to be obtained, and observations are to be extended to the D, 
line of helium and the H and K lines of calcium. ‘The problem of making 
corrections for atmospheric scattering is to be pursued. 

TABLE IV 
Frequency Distribution of Prominences 
Number 


A 





c ~~, 
Intensity Range Present Observations Woolley and Newton 


I-4 10 10 
5-10 31 51 
11-15 16 25 
16-21 9 10 
>21 3 ° 
g. Results: prominences on the disk._'Yhe intensities of prominences photo- 
graphed as dark markings on the disk are recorded in ‘Table V. Column I gives 
the date of observation and column 2 the heliocentric angle (i.e. the angular 
distance from the centre of the apparent disk) of the part of the prominence 
whose spectrum was obtained. ‘The coefficient of limb-darkening, o, is given 
in column 3, and was found from the centre-to-limb variation of the continuous 
spectrum at wave-lengths near Ha as measured by previous workers (17). ‘The 
observed central intensity of the prominences is given in column 4. On each 
plate the intensity of the centre of H« on an undisturbed part of the Sun’s disk 
near the prominence was measured ; this is called the normal Ha and its observed 
value is recorded in column 5. In the majority of plates when the prominence 
was exposed to give a suitable density on the plate for photometry, the continuum 
was too dense for measurement. For this reason most of the intensities in 
columns 4 and 5 are absolute intensities, as defined in Section 3, that is, the 
intensities have been measured with respect to the continuous spectrum of the 
centre of the Sun’s disk, by using the standardization spectrum. In a few cases, 
however, it was possible to make reliable observations of both the prominence and 
the continuum on the same exposure, and thus to obtain relative intensities. 
Such measurements involve no errors due to differences in timing or sky conditions 
between exposure and standardization, and so are more accurate than the absolute 
intensity measurements. Observed relative intensities are underlined in the 
table. , 
By combining the corrections worked out in Sections 3 and 4 for ghosts and 
scattered light, we find the combined correction for the relative intensity (see 
equations (2) and (4)) to be 


=ir'(I+P)—(g+B)j/(t-8), (0) 
and for the absolute intensity (see equations (3) and (5)) 


i= {t’—o(g + B)}/(I—g). 
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These corrections have been carried out for all observed intensities and are 
entered in columns 6 and 7 of Table V for the prominences, and in columns 8 
and 9 forthe normal H«. Except in cases where relative intensities were observed, 
the corrected relative intensities were found from the corrected absolute intensities 
(as given by equation (7)) by using the relation i/r =o, which follows from the 
definition of the limb-darkening coefficient. When relative intensities were 
observed, they were corrected by equation (6) and then, from the relative 
intensities so corrected, the corrected absolute intensities were deduced by using 
the relation i/r=o. ‘The instrumental profile correction was not carried out, 
but was estimated to be equivalent to a reduction in central intensity of about 
24 per cent in the first, and 4 per cent in the second order. 


TABLE V 
Central Intensities of Prominences on the Disk 
2 3 4 5 6 7 8 9 
Helio- Corrected Intensity 


centric ¢ Sens Shamil ti Prominence Normal Ha 
angle Abs. Rel. Abs. Rel. 
(a) Observations from first-order plates 
1949 Mar. 56° 0°78 8-2 16°0 5°6 71 13°4 17°2 
1949 Mar. 42° 0°88 11'S 18°5 8-6 9°7 15°6 17°7 
1949 Mar. 42° 0°88 13°6 22°9 10°7 12°1 20°0 22°8 
1949 Mar. 38° 0°90 160 25°4 13°0 14°5 22°5 25°0 
1949 Mar. 44° 0°86 15°6 22°1 12°8 14°8 19°3 22°4 


Central Intensity 


(b) Observations from second-order plates 

1948 Mar. 13 0°99 19°5 27°3 9°8 9°9 17°8 18-0 
1948 Mar. 34 o'92 189 28°3 9°9 10°8 19°5 212 
1948 Sept. 46° 0°85 19°4 25°4 11°2 13°2 17°3 20°4 
1948 Sept. 40 0°89 20°7 26°8 12'1 13°6 18°3 20°5 
23°3 30°2 13°7, 15°4 20°5 23°0 
1948 Sept. 46° 0°85 20°5 29°9 12°3 1474 21°99 258 
1948 Sept. 53° 080 22°5 30°2 14°8 18°5 22°7 28°4 
1949 Mar. 42° 0°88 12°9 19‘2 4°2 4°8 10°6 12° 
1949 Mar. 42° 0°88 17°97 27°0 g'l 10°3 186 8 2i'1 
1949 Mar. 51° 0°82 14°6 19'I 6°5 8-0 =—sar1'2 13°7 
1949 Mar. 53° o'80 13°7 24°0 58 7°3 164 20°5 
1949 Mar. 44° 086 16°8 22°7 8°4 9°7 14°4 16°7 
1949 Apr. 48° 0°84 20°7 26°5 12°6 15‘O 18°5 22°0 
1949 May 65° o*70 27° 30°5 13°77 19°5 163 23°3 
1949 May 090 19°8 24°2 Ir‘O 12°3 15's 17'2 
1949 May 39° 089 18°6 26°3 9°9 II‘! 17°8 20°0 
1949 Aug. 55° 14°5 20°9 6°8 8-6 13°3 16°8 
1949 Sept. 50° 0°83 16°9 22°7 8-8 10°6 14°7 17°7 
17°7 23°6 7°7 93 131 15°8 
1949 Sept. 50° 0°83 16-2 21°! 8-2 9°9 13°1 15°8 
19°4 25°1 93 +12 #144 «17% 

Columns 8 and g of ‘lable V show a scatter in the central intensity of the 

normal H« which is too large and bears no relation to the heliocentric angle, and 

therefore cannot be accounted for by the centre-to-limb variation. Furthermore, 
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a comparison of columns 6 and 8 (or columns 7 and 9) shows a correlation between 
the prominence and normal intensities; for instance, where a particularly low 
value is recorded for prominence intensity a correspondingly low value is found 
for the normal Hz intensity on the same plate. This is probably the result of 
inaccurate timing of the exposures; prominence and standardization spectra 
were taken separately, and the exposures were rather too short (5 seconds) to 
allow for accurate timing with a hand-operated camera shutter. Part of the 
error may also be due in some cases to changed sky conditions between the 
prominence exposure and the standardization. ‘The intensities recorded are 
therefore to be taken only as preliminary; the table has been included mainly 
as an illustration of the order of magnitude of the ghosts and scattered light 
corrections, which are seen to be considerable. 
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SOME EXPERIMENTS ON THE SCINTILLATION OF 
STARS AND PLANETS 
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(Communicated by the Astronomer Royal for Scotland) 


(Received 1951 November 6) 


Summary 


Cathode-ray oscillograms of stellar and planetary scintillation have been 
recorded over a period of a year. For this investigation selected multiplier 
photocells (type RCA 931—A) of special sensitivity and high signal/noise 
ratio have been used at the Cassegrain focus of the 36-inch reflector of the 
Royal Observatory. ‘The output voltages, after amplification, are applied 
to the beam of an oscillograph, the ‘‘ Y-deflections”” of the spot being 
recorded on moving film. ‘The sensitivity and dark current of the cells and 
the time constant of the circuit are discussed 

Scintillation in integrated light is exhibited as a function of altitude and 
of seeing conditions. Scintillation amplitudes of 5 to 10 per cent are always 
present in the zenith, when recording with the full aperture of 36 inches. 

The amplitude of scintillation is found to vary with telescope aperture A 
approximately in the ratio 1/4, but appears to reach a maximum at an aperture 
near 3 inches. Under normal conditions at this aperture, high-altitude stars 
show sudden flashes of light, lasting about one-hundredth of a second, which 
cause momentary increases in brightness of 1 to 2 magnitudes. 

The average size (3 inches) and distance (5 km) of the atmospheric 
irregularities causing scintillation have been deduced from a comparison of 
the scintillation patterns of stars, planets and satellites. 

Some experiments on colour scintillation, using ‘‘ monochromatic ”’ 
light, are described. An attempt has also been made to imitate the main 
features of scintillation by generating refractional irregularities in a model 
atmosphere of water. 


“ 





1. Introduction.—It is convenient to classify the effects of atmospheric 
turbulence upon star images under three headings: 

(a) Changes of apparent brightness (brightness scintillation) ; 

(6) Changes of colour (colour scintillation) ; 

(c) Changes of position, including “ dancing”’ and “ pulsation’. 

Aristotle, in De Caelo (1), refers to “‘ the apparent twinkling of the fixed stars 
and the absence of twinkling in the planets’’.* As we shall see (Section 2c), 
this statement is true of the major planets, not merely of observations made with 
the unaided eye but also for instruments of much larger aperture: it also applies 
over the whole recorded range of scintillation frequencies. 

The use of a 1P21 multiplier photocell, whose output is amplitied and applied 
to a cathode-ray oscillograph, has made possible an exact study of the intensity 
variations and their frequency in scintillating starlight. By reason of the exceed- 
ingly small time constant (~10 microseconds) which is attainable in such a circuit, 

* We are indebted to Professor W. B. Stanford, F.T.C.D., for finding this reference. Aristotle’s 
notion of the cause of scintillation, though attractive, is scarcely consistent with modern ideas in 


optics. He says ‘‘ The planets are near, so that the visual ray reaches them in its full vigour, but 
when it comes to the fixed stars it is quivering because of the distance and its excessive extension ”’ 
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the equipment is admirably suited to cover frequencies of scintillation ranging 
from those which are so low as to be easily perceptible to the eye, say, I cycle per 
second, to the highest values recorded, which are in the region of 1000c/s. The 
first recordings of this type, using an electromagnetic oscillograph, were made by 
Whitford and Stebbins (2) in 1936, but the long time constant of their equipment 
would not permit investigation of frequencies higher than 20¢;s. 

Changes in colour are easily discernible * in bright stars seen at low altitude 
with the unaided eye. From this we may infer that the fluctuations of brightness 
are not synchronous at all wave-lengths throughout the visible spectrum, and that 
at any instant some colours fail to reach the eye, either through interference or by 
refractional de-focusing effects. Colour scintillation can be more readily observed 
by a simple method due to Montigny (3), who mounted behind the eyepiece of a 
small telescope a thick plane-parallel plate whose normal was set at a small angle 
to the optical axis of the telescope. When this plate is caused to spin about the 
telescope axis at five to ten revolutions per second the star image is spread out 
along the circumference of a circle, and its light is thereby separated into its 
time elements. Colours of many hues are observed at different points on the time 
base. ‘The same phenomenon may be demonstrated with even simpler apparatus, 
by observing the star with a pair of field glasses held in the hands, while giving to 
the objective end a rapid rotatory motion (4). 

More precise information about colour scintillation may be obtained by means 
of the objective prism experiment described in Section 4, by Zwicky’s photo- 
graphic method (5), or by means of the multiplier photocell equipment used with 
suitable narrow-band filters. 

Lateral movements of star images have been measured photographically by 
Schlesinger (6). ‘These movements have a period which is roughly of the order 
of one minute of time and an amplitude of about one second of arc. ‘They have 
been recorded even near the zenith and under good conditions of seeing: they are 
also familiar to all astronomers who are accustomed to accurate guiding upon a 
star when the image is bisected by a micrometer wire. Such movements are not 
to be confused with the much more rapid dancing and pulsation which are observed 
under conditions of bad seeing. It seems probable that the former arise from the 
interposition in the beam of large wedge-shaped air masses of abnormal refractive 
index, having dimensions of the order of 100 metres. Dancing and pulsation, on 
the other hand, having frequencies comparable with those of scintillation and 
colour change, most probably have a common origin in the small-scale turbulence 
(dimensions ~ 10cm) which we believe to be the cause of scintillation. 


2. Brightness scintillation in integrated light 


(a) Instrumental arrangements.—-A multiplier photocell (type RCA 931~A) 
was built into a light-tight tube 15 inches long ( Fig. 1) which can be inserted in the 
converging beam of the 36-inch Cassegrain reflector of the Royal Observatory. 
This instrument has an equivalent focal length of 54 feet. ‘The rays are deflected 
by means of a g0° prism and brought to a focus in the plane of an iris diaphragm, 
which serves both for centring the image and for guiding. ‘The iris aperture is 
kept large enough (0-16inch, equivalent to 50") to ensure that its edge will not 

* In the well-known nursery rhyme, the twinkling star is described as “‘ like a diamond in the 


sky’. This is an apposite comparison; for a well-cut diamond held in a beam of sunlight in a 
darkened room exhibits just such flashes of sparkling colour. 
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occult any part of the star’s light through dancing and turbulence of the image. 
Guiding is effected by light reflected from the unsilvered face of a second go° prism ; 
the main beam traverses this prism and passes through a holder for a neutral screen 
or filter. It is finally focused as a spot 5 mm in diameter upon the sensitive surface 
of the photocathode. Precise focus upon the cathode is impossible because of the 
inferior glass in the cell walls; and it is in any case undesirable, in view of the 
possibility of local variations in the electron emission sensitivity of the cathode 
surface. 
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Fic. 1.— Arrangement of the photomultiplier and its circuits at the 
Cassegrain focus of the 36-inch reflector. 

The cell output is passed through a high resistance (10°ohms) to earth. 
The voltages across the ends of this resistance are amplified and registered as 
““Y-deflections”’ of the light spot on a cathode-ray screen. A mechanical contact 
is arranged to short-circuit the ends of the resistance at intervals of 0%-8.* This 
device continuously imprints on the record both the time-scale and the level of 
zero light intensity (cf. Figs. 3, 6, 9). 

The Y-deflections are recorded photographically by means of a camera which 
operates a continuously moving film: the film runs through horizontally at a speed 
which is normally about 6 inches per second. In later work a Cossor oscilloscope 
(Type 1049) with d.c. amplification has been used. ‘This instrument has the 
advantage over the a.c. amplifier of maintaining the zero level at a fixed point on 
the screen, despite changes in the amplifier gain. ‘The Y-scale may be varied by 


* In the earlier work a pendulum contact was used, giving time intervals of o*-5. In the present 
arrangement a wheel carrying a spring contact is driven by a reduction gear from the synchronous 
mains motor of the camera, and this contact makes the short-circuit (of duration 5 milliseconds) 
at o* 8 intervals. 
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control of the accelerating voltages applied to the phototube, by means of the 
amplifier gain, or by inserting neutral glass screens in front of the phototube. 

Power for the photomultiplier was, for the earlier experiments, provided by dry 
batteries, but owing to the high initial cost of these, their tendency to develop a 
variable internal resistance, and the danger of having a power supply of such lethal 
capabilities in the darkened dome, a mains-driven regulated high-voltage supply, 
based on the design of DeWitt and Seyfert (7) was constructed. This supply 
delivers any voltage required between 400 and 1600 volts with current up to 
about I milliamp. When adjusted to 800 volts, variations in output are less than 
tvolt. Switching and overload precautions are effected entirely by relays. The 
high voltage is conveyed to the photomultiplier by coaxial cable with polythene 
dielectric. 

(b) Observational procedure.—Over two hundred records of scintillation have 
been made with the above apparatus during a period of one year, under a wide 
variety of weather conditions and at many different zenith distances. Notes have 
been kept of the degree of scintillation observed visually at the time, of the appear- 
ances of the star images in the 7-inch guiding refractor attached to the 36-inch 
reflector, and of general weather phenomena, such as temperature, surface wind 
speed and direction. 

The photocells used in this work were specially selected for high sensitivity and 
high signal/noise ratio from a large number purchased from government disposal 
centres. ‘They are, therefore, equivalent to the 1P21 cells so extensively used in 
stellar photometric work. Our best cell, for example, gave a current of 
2:5 x 10°°amp from Deneb (A type, mag. 1-33, altitude 49°) with goo volts across 
the ten stages of the cell, using the full 36-inch aperture of the telescope. ‘The 
dark current was measured as 5x 10-Samp. ‘The time constant of the whole 
apparatus, due to the 21-foot coaxial cable and the input resistor of the oscilloscope, 
was 0°13 millisecond. ‘lhe timing short-circuit, of duration 5 milliseconds, was 
also found to give a perfectly square pulse on the film record. ‘The oscilloscope is 
stated by the makers to produce no distortion so long as the wave form is less than 
10°c/s. It seems probable that scintillation frequencies are normally well below 
this value. 

It was thought desirable to check the uniformity of the response from different 
parts of the photocathode in each cell, so as to eliminate the possibility of sharp 
irregularities in the response from adjacent small areas of the sensitive surface. 
This investigation was kindly carried out by Mr M. J. Smyth.* His work makes 
it clear that, so long as the starlight is not focused sharply upon the cathode 
surface, minute image shifts alone cannot contribute any appreciable irregularities 
to the cell output currents. 

‘The proportionality between spot deflection on the cathode-ray screen and 
light intensity incident upon the photocell has been established for all cells used in 
this work, (a) by the use of neutral screens, or () by a variable-aperture method. 
It was established that the response is linear, so long as the output voltage of the 
cell is not allowed to exceed 60 per cent of the voltage applied between each 


pair of dynodes. This condition obtains at the light levels experienced in the 
present work. 


Fig. 2 exhibits the spectral sensitivity characteristic curve for the 931-A 
photocells having an S-4 response. ‘This curve has been plotted from data 


* See p. go. 
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supplied by the makers, and is given for equal values of radiant flux per angstrom 
at all wave-lengths. ‘The prisms and lenses of crown glass, necessarily employed 
to focus the beam, will somewhat reduce the response towards the violet end of 
this curve. It is clear that the main response is to green, blue and violet light with 
a much lower response to yellow and red. 

The film records have been examined by projection through a photographic 
enlarging camera, the images being focused on to a white screen. ‘T'racings of 
each record were made upon sheets of squared paper, altering the magnification, 
when necessary, so as to exhibit all records with equal intensity ordinates 
(cf. Figs. 3, 6,9). From each tracing the amplitude of maximum light variation 
was read off and expressed as a percentage of the mean light level ( + above, and 
— below the mean level). An estimate was also made of the maximum frequency 
of the light fluctuations in cycles per second. 
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Fic. 2.—Spectral sensitivity curve for a 931-A photomultiplier, as given by the makers, for 
equal values of radiant flux per angstrom at all wave-lengths. Effective filter transmissions shown 
by vertical lines. 

Ordinates: Relative sensitivity. 
Abscissae: Angstroms. 


(c) Summary of results.—Using the full spectral response of the cell and the 
full aperture of the 36-inch reflector: 

(1) On all nights the amplitude of scintillation decreases with increasing 
altitude of the star (cf. Fig. 3, Nos. 1-6). In Fig. 4 the connected points show the 
trend of the amplitude—altitude relation for different nights and for different 
conditions of seeing. ‘These curves show that even on nights of steady seeing the 
amplitude of scintillation in the zenith is of the order of +5 per cent with the full 
36-inch aperture, and, as we shall see later, it is much greater than this for smaller 
apertures. At an altitude of 10° it lies within the range + 20to +50 percent. In 
Fig. 5 the amplitudes are plotted against sec Z (Z = zenith distance in degrees). 
These curves illustrate more clearly the convergence to a scintillation value of 
+5 to +10 per cent in the zenith on an average night, recording with the 36-inch 
aperture. It may be emphasized that this is not an instrumental effect, since the 
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cell dark current is quite negligible at these high light-levels (see, for example, 
Fig. 9, No. 1). 

(2) In low-altitude records (<10°) low-frequency scintillation (~5c;s) 
predominates, though higher frequencies (~100c/s) are usually present: with 
increasing altitude the low-frequency fluctuations die out while the high-frequency 
component remains, even up to the zenith. 

(3) The amplitude of scintillation appears to be closely correlated with the 
quality of the star images, as observed visually with the 7-inch auxiliary refractor. 
When the diffraction disk and rings are visible in this telescope, the scintillation 
amplitude does not exceed 10 per cent in the larger instrument: amplitude 
variations of 30 per cent or more are associated with the most violent movements 
and boiling of the visual image. 














Fic. 3.—Cathode-ray oscillograms of scintillation in integrated light using the 
full 36-inch aperture of the telescope. 
Time horizontally : ordinates proportional to light intensity. 
. Deneb at 68° altitude: zero marks 08-8 apart. (Film No. 138.) 
. Vega at 49° altitude: zero marks 08-8 apart. (Film No. 137.) 
. Altair at 32° altitude : zero marks 08:8 apart. (Film No. 139.) 
. Altair at 25° altitude: zero marks 08-8 apart. (Film No. 149.) 
. Betelgeuse at 13 altitude: total time 08-8. (Film No. 1§2.) 
6. Betelgeuse at 7 altitude: total time 15-0. (Film No. 150.) 
These records were obtained on the same night during a period of 14 hours. 
(1n Figs. 3, 6 and 9 the horizontal straight lines give the zero levels of intensity.) 
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(4) Under poor conditions of seeing, when light from a star is fluctuating with 
an amplitude of +20 per cent, that from Jupiter or Saturn remains constant to 
less than I per cent (see Fig. 9, No. 1). 

The following conclusions are drawn from records obtained with the full 
spectral response of the cell and using a variable telescope aperture. For these 
measurements a series of circular stops were used. These were placed across the 
mouth of the telescope tube, centred upon the circumference of a circle of 9 inches 
radius so as to avoid as far as possible the shadow of the 10-inch secondary mirror. 
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(5) Decreasing the aperture of the telescope increases the amplitude of scintil- 
lation for all altitudes (Fig. 6, Nos. 1-5). ‘The plotted points of Fig. 7 illustrate the 
amplitude—aperture relation for a number of stars at various altitudes and for the 
planet Venus (31” apparent diameter, crescent 12” wide). 

or 
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Fic. 4.—Amplitude of scintillation in integrated light as a function of altitude in degrees. 
Recordings made with the full 36-inch aperture of the telescope. Records for different nights are 
shown by characteristic symbols. 
Ordinates : Altitude in degrees. 
Abscissae : Amplitude of scintillation as percentage of mean light level. 





(6) The most remarkable features in the stellar records from small apertures 
(of 3 and 6 inches) are the sudden increases of light intensity, of 200 to 400 per cent 
above the mean level and having durations of the order of 0*-o1 (Fig. 6, No. 5). 
This phenomenon has also been reported by H. E. Butler (11). ‘These flashes are 
usually most pronounced at an aperture of about 3 inches, and there is a suggestion 
in the shapes of the low-altitude records of Sirius and Rigel that at this aperture 
we are approaching the maximum amplitude of scintillation. Unfortunately, we 
cannot work at still smaller apertures to make certain of this conclusion, for we 
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have here reached such a low light-level that the cell noise, as distinct from the 
‘seeing’ noise, is beginning to predominate. We are also at the limit of satis- 
factory amplification with the existing apparatus. 

However, it seems clear from these results that we have to deal with a pattern 
of light variations crossing the telescope mirror, the effective size of the disturbances 
being of the order of 3 inches rather than 36 inches. Shadowgrams of such waves 
crossing the mirror of the Cordoba 60-inch reflector have been photographed by 
FE. Gaviola (8), and more recently also by Mikesell, Hoag and Hall *, using the 
40-inch reflector of the Naval Observatory, Washington. Observations of the 
dark shadow-bands seen at total eclipses of the Sun also give widths of 2 to 6 inches, 
separated by distances.up to I or 2 feet.t 
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Fic. 5.—Amplitude of scintillation as a function of sec Z. From the same data as used in 
Fig. 4. Note the convergence to a value of § to 10 per cent in the zenith (sec Z=1). 
Ordinates: (on left) sec Z. 
(on right) Altitude in degrees. 
Abscissae: Amplitude of scintillation as percentage of mean light level. 


The general effect of using a large aperture in our experiments is to smooth out 
the light variations, as shown in Fig. 7. We may call this the “integration effect’’, 
and we can predict that with apertures of the order of 100 inches brightness 
scintillation should be negligible. 

3. The size and location of the atmospheric irregularities.—The evidence 
summarized in the previous section leads to an estimate of 3 inches as the order of 
the average true size of the atmospheric irregularities responsible for the scintillation 
phenomena we have been discussing. So far, we have visualized scintillation as 
brought about by the shadows of these irregularities cast upon the telescope mirror 
or objective by stellar point sources of light. When the apparent diameter of the 
source of light becomes appreciable, as for the planets, the moving shadows cast 


* Privately communicated. 
t See, for example, Pernter and Exner, Meteorologische Optik, p. 166, 1910. 
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upon the telescope aperture will lose their sharpness and the phenomena of scintil- 
lation will become less marked. Light from each point in the extended source will 
still be subject to scintillation, as from a star, but rays emanating from widely 
separated points upon a planet’s disk will traverse different atmospheric irregular- 
ities whose variations will be out of phase with one another. ‘Therefore, we should 
expect brightness scintillation to begin to die away at some stage as the apparent 
diameter of the source increases, just as, when using a point source (star), the 
scintillation decreases with increasing telescope aperture (integration effect). 





5 

















Fic. 6.—Scintillation patterns in integrated light, recorded with telescope apertures from 
36 inches down to 3 inches. All refer to Altair (mean altitude 28°) within a period of 45 minutes. 
Time horizontally : ordinates proportional to light inter-ity. 
1. 36-inch aperture. (Film No. 139.) Time marks 08-8 apart. 
2. 18-inch aperture. (Film No. 140.) Total time 08-7. 
3. 12-inch aperture. (Film No. 141.) Time marks 95:8 apart. 
4. 6-inch aperture. (Film No. 145.) Total time 08-33. 
5. 3-inch aperture. (Film No. 144.) Note the greatest peak four 
times the mean level. Total time 08-33. 
See also Fig. 7. 


It seems clear that this stage must begin to be reached once the apparent dia- 
meter of the source becomes equal to that subtended at the receiver by the irregular- 
ities concerned. ‘To investigate this problem we have performed numerous 
experiments using planets as sources. In each case the amplitude of scintillation 
has been measured and compared with that of a nearby star, both being taken at an 
altitude of about 30°. Jupiter, apparent diameter 49”, and Saturn, apparent 

6 
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diameter 20”, showed no measurable scintillation with 36 inches aperture, while 
the two comparison stars gave scintillation amplitudes of +20 per cent and 
+10 per cent respectively. On another night, scintillation records were obtained 
of Jupiter, Jupiter II] and 20 Tauri, all at the same altitude of 33°. These 
traces (Fig. 9, Nos. 1, 2, 3) indicate that while the planet’s light again is steady, 
the satellite scintillates precisely like a star of the same magnitude. ‘The apparent 
diameter of Jupiter II] was 1”-5 at the time. 
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—Amplitude of scintillation plotted against telescope aperture for various nights and altitudes. 
Ordinates : Aperture of telescope in inches. 
Abscissae : Amplitude of scintillation as a percentage of mean light level. 


We infer from this that scintillation begins to die out somewhere within the 
range of apparent diameter 1"-5 to 20”. But we have been able to narrow the 
range still further by observations of the planet Mars. At this point we must 
emphasize the fact that the apparent diameter at which scintillation of a planet 
begins to disappear is also dependent upon the telescope aperture in use. For 
example, in Fig. 7 we have plotted the amplitude of scintillation of Venus (altitude 
21°) for telescope apertures varying from 36 inches downtotlinch. ‘The apparent 
diameter of the planet on this occasion was 31”, but the width of the illuminated 
crescent was only 12”, so that we might regard the planet for scintillation purposes 
as equivalent to a circular disk of about 20". With the full aperture there was no 
trace of variation in brightness, but scintillation was recorded at 18 inches aperture 
and increased steadily to a value of +23 per cent at an aperture of I inch. 
(Because of the great brilliance of the planet it was possible to work at these small 
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apertures.) It is, consequently, important that all such comparisons between stars 
and planets should be made with the same telescope aperture. A simultaneous 
record of Mars, which was only 4° away from Venus and had at the time a disk 
4”-4 in diameter, gave scintillation of + 14 per cent with the full 36-inch aperture, 
as compared with +23 per cent for Sirius at a comparable altitude on the same 
night. We are thus led to conclude that when the apparent diameter of a planet 
is in the region of 3” its scintillation will begin to decrease relative to that of a 
star at the same altitude, both recorded with the same telescope aperture. 

We can now make a tentative estimate of the distance from the telescope of the 
atmospheric irregularities concerned. We assume that their average diameter is 
3 inches (7-6 cm) and that they subtend an angle of 3": from this we arrive at a mean 
distance of 5-2km. ‘This figure may well be in error by a factor of 2, but it is of 
interest to note that it is of the same order as the value of 4 km derived by Rayleigh 
(9) from his discussion of the refraction theory of scintillation which was based 
upon the early visual observations of Montigny and Respighi. ‘This theory is 
referred to more fully in Section 5. 
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Fic. 8.—Schematic arrangement for the investigation of the Respight spectral band 
phenomenon (destruction of colours by scintillation). 


4. Colour scintillation.—The fact that we can observe colour changes in low- 
altitude stars with the unaided eye or with a small telescope indicates that, at any 
rate for small apertures, little or no phase coherency exists between the scintillations 
in different wave-lengths. ‘That is to say, the various colours must fluctuate in 
brightness independently of one another in time; for if all wave-lengths were to 
scintillate with equal amplitude and in phase we should expect to observe changes 
of apparent brightness in the integrated image but no changes of colour. 

To investigate the matter further, a 6-inch objective prism was placed in front 
of a 5-inch aperture objective of 25 inches focal length (Fig. 8). ‘The spectrum of 
Sirius was focused upon a fine-ground glass screen and observed through a magni- 
fying lens. ‘The lens could be set in vibration in a plane at right angles to the 
dispersion with a frequency of about 20 oscillations per second. In this way the 
time sequence of the colour changes could be rendered visible to the eye (4). At any 
instant these elementary spectra are seen to have one or more broad bands of colour 
missing: some exhibit only red and violet, others yellow and green, and so on. 
‘These dark bands, where the colours are destroyed, or partly so, are of the order of 
1000A wide and chase one another along the length of the spectrum in an 

6* 
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apparently haphazard fashion. ‘lhe frequency of fluctuations at any one wave- 
length is variable, but is in the neighbourhood of 50c/s.* It seems, therefore, 
that with this aperture also no phase coherency exists between the scintillations in 
different wave-lengths. ‘This is not to say that small fluctuations may not also 
occur in the apparent brightness of the spectrum as a whole, affecting all colours at 
the same time and to the same extent. 

Since this experiment was first described (4), F. Zwicky has published (§) some 
very beautiful spectrograms demonstrating the details of the same phenomenon 
His photographs were taken by means of an objective prism (and, in one case, a 
large replica grating) mounted in front of the corrector plate of an 18-inch Schmidt 
camera. The time base was provided, either by the diurnal drift of a bright star 
across the field, or by swinging the camera by hand in a direction at right angles to 
the plane of dispersion. ‘These photographs also record the presence of dark 
bands moving rapidly along the spectrum, showing that individual colours are 
successively obliterated or greatly reduced in brightness. 

Let us consider a simplified model in which the amplitude of scintillation at each 
separate wave-length is the same fraction x of the mean level at that wave-length. 
Then if the oscillations in intensity at the various wave-lengths are in phase, the 
amplitude of scintillation in integrated light will be the same fraction x of the mean 
integrated level. But if the oscillations in intensity at various wave-lengths are 
not in phase, then, for integrated light, the amplitude of scintillation will be less 
than x times the intensity of the mean level. 

It was decided to put this matter to the test, using the photomultiplier apparatus 
attached to the 36-inch reflector, in conjunction with various ‘‘ monochromatic ”’ 
filters. ‘The procedure was to make first a scintillation record in “ white”’ light, 


using the full range of colour response of the cell (see Fig. 2), but with a neutral 
screen interposed in the beam. A ‘“‘ monochromatic” filtert was then substituted 
for the neutral screen and the monochromatic record was made immediately 
afterwards. No change of gain of the C.R.O. amplifier was allowed between the 
records. In this way, by choosing.a neutral screen having the same effective 
transmission as the filter, one could eliminate the possibility of any instrumental 
alteration between the two records. ‘The results of these tests are entered in 


“ce 


Table I, where the amplitudes of scintillation are compared for “ white” and 
‘monochromatic ” light at various altitudes. Amplitudes are given as a percentage 
of the mean level, taken as 100 per cent. 

There is thus little difference between the amplitude of scintillationin “ white”’ 
light and in either of the two “‘ monochromatic”’ bands which were used, from which 
we infer that the fluctuations in intensity are probably coherent in adjacent 
wave-lengths. This conclusion is also in agreement with direct measures 
obtained recently by Mikeseli, Hoag and Hall {, using a split-beam photometer 
attached to the 40-inch reflector at the Naval Observatory, Washington, and 
recording simultaneously in two wave-lengths 800A apart. For the reasons 
given above it seems that this conclusion can scarcely be true for smaller apertures, 

* An earlier investigation of this ‘“‘ band’? phenomenon was carried out by Respighi (Assoc. 
Frangaise pour l’ Avancement des Sciences, 1, 148, 1872). A summary of his results is to be found in 
Physics of the Air (W. J. Humphreys, 1920), and in Meteorologische Optik (Pernter and Exner, 1910). 

t The following filters were used: a “ mercury violet’, transmitting a band approximately 
500 A wide, centred at 4500 A; a “‘ mercury green”’, passing 500 A, centred at 5 300 A, and a 


Schott O.G.2 transmitting from 5500 A to 7000 A (see Fig. 2). 
} Privately communicated. 
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in the range 3 to 6 inches, and it is just here that cathode-ray oscillograms taken 
simultaneously in two colours would be of the greatest interest. Unfortunately, 
there is insufficient light available with these small apertures to enable us to record 
monochromatically with the existing apparatus. 


TaBLe | 


Amplitude of Amplitude of 
scintillation in scintillation in Filter Star Altitude 
“ white ” light “* monochromatic ”’ 
light 
+15 +14 Mercury green Sirius 
15 16 Orange green 2 Sirius 
7 10 Orange green 2 Betelgeuse 

15 23 Orange green 2 Sirius 

30 30 Mercury green Sirius 

20 20 Mercury green Spica 

20 23 Mercury green Spica 

6 Orange green 2 Arcturus 


Means + 16-0 per cent +17°8 per cent 


5. Theories of scintillation.—The refraction theory of scintillation was first 
proposed by Montigny. His “ Principle’? has been summarized as follows: 
“The individual colours, which go to make up the uniformly coloured or ‘ white’ 
image seen by the eye, travel through the atmosphere along paths which become 
more separated the higher one follows them upwards. ‘These path separations are 
large enough at a certain height to enable small atmospheric irregularities to cross 
the ray paths of the separate colours in succession and to deviate these colours from 
their original direction independently of one another, so that first one colour, then 
another, is removed (or partly so) from vision. ‘Thus the integrated star image 
shows to the eye not only colour changes but also brightness changes.” * 

This theory received the general approval of Rayleigh, who developed it in 
mathematical form (g), so as to provide a satisfactory explanation of the spectral 
band phenomenon, first described by Respighi. Rayleigh calculated that at the 
top of the atmosphere the separation of the red and violet rays contributing to 
the light reaching the point O (see Fig. 10) would be of the order of I to 2 
metres if the apparent altitude of the star were 20°. But he also gives reasons 
for supposing that, in fact, the irregularities producing these deviations are 
much lower down in the atmosphere—at a distance from the observer of about 4 km 
and in dimensions about 4cm across. Clearly, on this model, colour scintillation 
would be less marked at smaller zenith distances, which is in good accord with 
observations, since the whole spectrum rather than single colours would then be 
affected simultaneously by the passing irregularities. Moreover, the degree ot 
scintillation should be greatest at small apertures (~4cm), the effect of a larger 
aperture being to include the pencils of rays refracted by a large number of 
irregularities simultaneously and thereby to smooth out in the combined effect the 
larger variations of colour and brightness. 

One of the principal difficulties confronting this theory, as Rayleigh himself 
pointed out +, is to show how sufficiently refracting irregularities can occur in the 
atmosphere, without having to assume improbably large temperature differences 


* Translated from Meteorologische Optik, p. 175. 
t Rayleigh, loc. cit., p. 139. 
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between adjacent air masses. C. G. Little has recently reviewed this difficulty 
(10), which he considers to be insuperable, and he abandons the refraction theory 
in favour of a diffraction theory of scintillation. ‘This may be described as a 
revival, in an improved form, of Arago’s interference theory (12) which Rayleigh 
believed he had killed.* Little’s theory is based upon a consideration of the 
distortion of the Fresnel half-period zone pattern of a plane wave viewed through a 
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FIG. 9. 

1. Scintillation pattern of Jupiter (apparent diameter 46") at 34° altitude with 36-inch aperture. 
(Film No. 157.) 

2. Jupiter IL] (mag. 5-3, apparent diameter 1":5). Altitude 33°. (Film No. 156.) 

3. 20 Tauri. Altitude 33°. (Film No. 158.) 

These three records were taken within a few minutes of one another. 
time marks in 1, 2 and 3 are 8. 

4. Scintillation of an artificial star shining through water 10 cm deep in a tank. Waves were 


induced upon the surface of the water by means of a tuning fork. (Film No. 105.) Time intervals 
here are 085. 


The intervals between the 


5. The flashes of a neon lamp recorded with the scintillation apparatus. The successive peaks 
are o8-o1 apart. (Film No. 174.) Note that the intensity of each flash depends upon the 
polarity of the alternating current. 

medium of slightly non-homogeneous refractive index. It has the merit of 
explaining the observed degrees of brightness and colour scintillation, assuming 
temperature gradients only one-thousandth of those required on the basis of 
refraction. On the other hand, it is not obvious that it can account for the complex 
spectral band phenomena of Respighi so well as the earlier ideas. 

We have investigated experimentally the question whether it is possible to 
imitate stellar scintillation effects by the generation of refractional irregularities in a 
* Rayleigh, doc. cit., p. 129. 
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model atmosphere of water. For this purpose a tank was used containing water 
to a depth of tocm. Small waves were generated in the water surface by the 
action of an electrically operated tuning fork (17¢/s) and paddle. Aon artificial 
star, consisting of a d.c. ‘‘ Pointolite’’ arc was placed 6 feet above the tank and the 
light from this ‘star’? was focused upon the cathode of the photomultiplier 
through a glass window in the flat bottom of the tank. 

Waves of small amplitude travelling outwards from the paddle were found to 
produce brightness scintillation patterns of the type shown in Fig. 9, No. 4, which 
give a very fair imitation of the scintillation recorded from a real star at low 
altitude (cf. Fig. 3, Nos. 4-6). It is interesting also to note that in Fig. 9, No. 4, 
the effect of interference between two wave trains can be clearly seen, one series of 
waves moving outwards from the paddle and the other, of lesser amplitude, 
returning after reflection from the walls of the vessel. 
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Fic. 10.— The refraction theory of scintillation. 


We do not wish to maintain that these simple experiments necessarily prove 
that brightness scintillation is caused by refractional focusing effects in the 
atmosphere, but they do suggest that such refractional irregularities may contri- 
bute largely to the scintillation effects we have recorded, especially at low altitudes. 
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NOTES ON THE 931-A PHO'TOMULTIPLIER 
M. F. Smyth 
(Communicated by the Astronomer Royal for Scotland) 
(Received 1951 December 5) 


Summary 


A number of ex-services RCA 931-A photomultipliers have been 
examined, and their characteristics are compared with those of the 1P21. 
The reduction of dark-current is discussed. A method of examining local 
cathode sensitivity reveals that some tubes have markedly uneven sensitivity 
and may be unsuitable for astronomical applications. ‘Types of stabilized 
power supply are discussed, and the effect of geoomagnetism on photomultiplier 
sensitivity is verified. 





1. Introduction.—In the past eight years multiplier phototubes have found 
many applications in astronomical photometry (1, 2, 3, 4, 5, 6); the specially 
selected RCA 1P2t1 in particular has been widely used. Owing to the high cost 
and difficulty of obtaining this tube, the similar RCA 931-A may be used, if a 
good specimen is selected. ‘This tube is available at a low price from ex-services 
equipment suppliers. A number of 931—A photomultipliers have been examined, 
and performance figures of the best tubes are comparable with those of the 1 P21. 


The method of selection is first considered. Necessary precautions are discussed 
in later sections. 

2. Selection of tube.—'The 931~A and 1P21 are small blue-sensitive phototubes 
with nine stages of multiplication (dynodes), and are mounted on moulded bakelite 
I1-pin bases. The 931~A is not specifically intended for applications involving 
low light levels, and specimens show widely different characteristics in this 
respect. A criterion of usefulness is the signal-to-noise ratio at low light levels. 
Instead of measuring the root-mean-square noise current directly, it was simpler 
to obtain an indication of its variation from one tube to another by measuring the 
anode dark-current, which is due to a combination of factors (see Section 3). 

Sixteen specimens were tested, using approximately 72 volts per dynode 
stage, and relative signal intensities and signal-current/dark-current ratios were 
measured. Dark-currents at about go volts per stage were also measured, at 
an ambient temperature about 18deg.C. ‘The results are shown in Table I, 
which illustrates the large variations found. It is remarkable that the first 
five tubes, which were purchased in light-tight metal holders, show lower dark- 
current than the remainder, which were purchased in cardboard boxes, although 
all tubes had subsequently been handled and exposed to daylight before testing. 

The dark-currents of the best tubes (e.g. No. 4) are comparable with the 
figure of 10-*amp given by DeWitt and Seyfert (6) for the 1P21, and below the 
maximum dark-current of 10 * amp quoted by the manufacturers. 

As a 1P21 was not available for testing, no direct comparison of luminous 
sensitivity could be made. Indirect comparisons can be made, but will be only 
approximate, since the sensitivity is critically dependent upon applied voltage, 
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and conditions of testing may vary. DeWitt and Seyfert (6) find that, using a 
1P21, a current of 10 ‘amp is produced by the light of a 4-5 mag. A-type star 
with a 12-inch reflector; while Ellison and Seddon (8), using the 36-inch reflector 
of the Royal Observatory, Edinburgh, obtain with the best 931—A tube (No. 4) 
a current of 2-5 x 10 °amp on a I-33mag. A-type star. This indicates slightly 
higher sensitivity for the particular 931—A. 
It thus appears that a selected 931—A tube should be equivalent in sensitivity 
and signal-to-noise ratio to a 1 P21. 
TABLE I 
‘Tube Relative Sensitivity Signal/Dark-current Dark-current at go V/stage 
Number (72 V stage) (72 V/stage) (in 10~* amp) 
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170 
270 

170 
220 
310 
310 
330 

130 

100 


Ove HN 


150 
190 

25 
190 
260 
180 
280 
35° 
300 
450° 


oo ws 
we WN - 
ooo Swen 0 SI NW 


3. Reduction of dark-current.—If the dark-current is large enough to introduce 
errors, steps must be taken to reduce it. ‘I'he sources of dark-current have been 
discussed by Engstrom (7) and by DeWitt and Seyfert (6). At low voltages 
ohmic leakage, particularly in the base of the tube, may predominate; amplified 
thermionic emission from the cathode becomes important at higher voltages ; 
and regenerative ionization effects may occur above 110 volts/stage. The latter 
effects need not be considered if the tube is operated below 100 volts/stage, as is 
desirable. Yates (5), DeWitt and Seyfert (6) and others have found that the 
thermionic dark-current may be reduced by a factor between 20 and 200 by 
cooling with solid CO,. But in most of the 931~A tubes tested, ohmic leakage 
was much greater than the other effects. ‘his was probably due to the proximity 
of the anode pin to the cathode pin, which may be at 1 000 volts lower potential. 
A considerable improvement, for some of the tubes, was produced by the technique 
of removing the anode pin from the base and insulating the anode lead, as 
recommended by Kron (3). The results of this operat,on for some tubes are 
given in ‘l'able II. 

"TABLE I] 
‘Tube Dark-current at 100 V stage 
Number Before removal of anode pin After removal of anode pin 
(10~-* amp) (107° amp) 

19 I 

30 8 

390 160 

140 2} 
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Signal-to-noise ratio may be expected to decrease with increasing applied 
voltage (7). However, one of the 931—A tubes measured showed a maximum of 
this ratio at 90 volts/stage; this was presumably due to the fact that leakage 
current was predominant at low voltages and increased less rapidly than the 
signal current with increasing voltage. Using a 931-A in a pulse-counting 
circuit similar to that described by Yates (§), an optimum voltage of 55 volts, stage 
was found for that particular combination. The pulse-counting technique, 
with a discriminator stage which removes smaller pulses due to amplifier noise 
and emission from the dynode stages of the tube, enables an improvement of 
signal-to-noise ratio to be obtained; also the sensitivity to small changes in 
overall voltage is diminished. 

In measuring the darlturrents, the tubes must be ‘“dark-aged”’ for some 
time before testing, as exposure to bright light renders them temporarily very 
‘“‘noisy’’ (3). Exposure to bright light when voltage is applied may, of course, 
permanently damage a tube. Over long periods of use with exposure to low 
light levels only, a steady improvement is found. 

931-A photomultipliers are also subject to fatigue after exposure to light; 
Engstrom (7) has found an average temporary loss of sensitivity of 50 per cent after 
40 minutes’ exposure with an output current of 100 microamps. ‘The effect 
varies from one tube to another; the author, using output currents about 
5 microamps, has found changes of about 1 percent in the relative sensitivity of two 
tubes within a few minutes’ exposu':. ‘Thus it is desirable in most cases to 
keep the output current as low as I microamp or less. 

4. Uniformity of cathode sensitivity.—In connection with the work of Ellison 
and Seddon (8) on the photoelectric recording of stellar scintillation, using a 
931-A photomultiplier, it was desirable to investigate the uniformity of the 
light sensitivity of photomultiplier cathodes. Large local variations in sensitivity 
would cause irregularities in the output current, due to minute image shifts, 
if the image focused on the cathode were too small. 

For this investigation an oscilloscope with short-afterglow tube was set up 
with its time-base adjusted to about 50 cycles per second, and with no deflection 
on the Y-plates. The oscillating spot of light was focused on to the rectangular 
cathode area of the photomultiplier, so as to scan a narrow strip of the cathode, 
parallel to its long dimension, 50 times per second. ‘The resulting current from 
the photomultiplier, at 75 volts/stage, was passed to the Y-plate amplifier of a 
second oscilloscope, whose X-plates were coupled to the time-base of the first. 
Thus the sensitivity along a narrow strip of the photomultiplier cathode was 
displayed on the second oscilloscope, the second beam of which was used to 
provide a zero line. 

Direct tracings, made with an enlarger from photographs of these patterns, 
are shown in Fig. 1. In this figure, (a) indicates the usual position of the scanned 
strip on the projected area of the cathode; for each tube the scanned strip was 
shifted to and fro parallel to itself, and for some tubes considerable lateral variations 
were thus found. Since the actual cathode surface is set at 45° to the direction 
of incident light, shifts along the shorter dimension of the cathode involve a 
change of focus. Fig. 1 (+) shows the trace obtained by focusing the scanning 
spot on a rhodium step-wedge and thence, using a field lens for uniform cathode 
illumination, on the photomultiplier; the response is sufficiently rapid and linear. 
Fig. 1(¢) is the trace for the tube with most uniform cathode sensitivity; the 
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maxima and minima arise from focusing the scanning spot on the wire grill in 
front of the cathode. Fig. 1(d) is the trace for our most sensitive tube (No. 4), 
and shows sensitivity uniformly increasing towards the lower end; the spot 
was again focused on the wire grill. Fig. 1(e) shows a central spot of nearly 
double sensitivity on an otherwise fairly uniform cathode; a small star image 
focused near this spot would give rise to large apparent intensity fluctuations. 
Fig. 1(f) shows a central region of low sensitivity ; the apparent sensitivity of 
this tube would be lowered if this region were used for measurements. Fig. 1(g), 
(A) and (2) are traces for the tube used in the stellar scintillation experiments, and 
show the response of parallel strips of the cathode about 3 mm apart. 


25mm al 
(a) raf) (c) 
TOP 




















(g) f\ t\ 


Fic. 1.—Tracings of oscillograms showing distribution of sensitivity along strips of the 
cathodes of 931-A photomultipliers. For explanation see text. 








Photographs of the cathode patterns of all our photomultipliers have been 
made and will provide a useful guide in their application to specific problems. 
No visible indication of high or low-sensitivity areas, as revealed in (e) and (f), 
can be seen on inspection of the cathodes. 

5. Power supplies.—The accelerating potentials for the dynodes may most 
simply be obtained from ten dry batteries. If variable potentials are required, 
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to control the gain, the dynodes may be fed from tappings on a potentiometer 
chain, in the usual manner. 

It should be noted that if a potentiometer chain is used, the currents drawn 
by the dynodes will modify the voltages; consequently, to preserve the intrinsic 
linearity of the photomultiplier (7), the potentiometer current should be at 
least 1000 times greater than the maximum (mean) anode current. This may 
impair the constancy of the voltage supplied by dry batteries; for this reason, 
and because of the bulk and lethal capabilities of batteries, a mains-driven stabilized 
power supply is often preferable. 

A unit is required which will supply about 1 milliamp at 1000 volts; since 
the sensitivity of the photomultiplier varies approximately as the ninth power of 
the overall voltage, a voltage stability of 1 part in ro* is desirable. Numerous 
circuits have been published (6, 9, 10, 11). ‘Two basic types of circuit which 
have been found useful by the author are illustrated in Fig. 2; no originality is 
claimed for these. ‘The rectifying and smoothing arrangements are conventional 
and have not been included; similarly, component values are omitted, since they 
may be chosen to suit individual circumstances. 




















a 


Fic. 2.—-Basic stabilizing circuits. Rectifying and smoothing circuits are omitted. 
R represents the potentiometer chain supplying the photomultiplier. 


Circuit (a) uses dry batteries, trom which no current is drawn, to stabilize the 
current through, and hence the voltage across, the photomultiplier potentiometer 
chain represented by R. ‘I'he current can be varied by adjusting potentiometer P ; 
e.g. if the batteries give 90 volts, and P is set at go000 ohms, the current is 
approximately 1 milliamp. ‘The circuit reduces input voltage variations, due to 
mains fluctuations, by a large factor. 

Circuit (6) is more flexible (cf. 6). ‘he 85A2, a miniature rare-gas voltage 
reference tube, is used at V,._ ‘This tube appears to be the most stable available 
(12,13); random voltage changes should not exceed 0-1 percent over periods of 100 
hours (14), and the voltage-current characteristic is very smooth above 1 milliamp. 
Other types of voltage-reference tube show irregular glow-shifts at low currents, 
and much greater voltage drifts. In this circuit, potentiometer P, adjusts the 
voltage across the photomultiplier potentiometer chain R, and P, can be adjusted, 
at each setting of P,, to provide complete stabilization against mains input 
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Huctuations; the settings of P, at each voltage may be recorded for reference. 
An EF-50 (VR-91) strapped as a triode has been used at V,. 

6. Positional effect—The manufacturers of photomultipliers recommend 
magnetic screening in the presence of strong magnetic fields; but it was found by 
Eggen (1§) that changes in orientation with respect to the Earth’s magnetic 
field, due for example to changing positions of a photomultiplier mounted on a 
telescope, can cause appreciable variations of sensitivity. ‘This ‘positional 
effect” is presumably due to defocusing of the electron streams, and was eliminated 
by Eggen by shielding with demagnetized iron. 

The author has versiied this effect, using a 931-A mounted in a box with a 
standard lamp. Litt!e effect was produced by changing the azimuth of the 
tube, but maximum variation was obtained between the cap-up and cap-down 
positions. ‘lhe excess sensitivity in the cap-up position, at constant output 
current, varied in this case between I1-5 per cent at 50 volts/stage and 9-5 per cent 
at gO volts/stage. ‘Ihe constancy of the standard lamp was checked by a barrier- 
layer cell. 

7. Acknowledgments.—I am indebted to Dr M. A. Ellison and Mr H. Seddon 
for valuable discussion and measurements with the 36-inch telescope, and to 
Professor W. M. H. Greaves for his continued interest and helpful suggestions ; 
also to the Moray Endowment of the University of Edinburgh for an equipment 
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STELLAR PARALLAXES DETERMINED PHOTOGRAPHICALLY 
AT THE CAPE OBSERVATORY (EIGHTEENTH LIST) 


J. Jackson and F. B. G. Turner 


(Communicated by H.M. Astronomer) 


(Received 1951 October 22) 


The principal list is in the same form as the previous lists beginning in M.N., 
89, 497, 1929. 

Approximately half of the stars in the present list are faint stars with a proper 
motion exceeding 0:5 annually, given by Luyten in Vol. III, No. 1 of the Publi- 
cations of the Astronomical Observatory of the University of Minnesota. About 
one-sixth are stars selected from the Cape Astrographic Zone as having an annual 
proper motion in excess of 0"-2, while most of the remainder were selected from 
the General Catalogue for a similar reason. ‘Three stars, Nos. 1635, 1670 and 
1671 appear on account of their relative proximity to Nos. 1634, 1669 and 1672 
respectively, and of the consequent uncertainty of the identity of the proper 
motion star in each case. 

Where possible, the magnitudes given are taken from the Henry Draper 
Catalogue ; otherwise these are generally as given by the authority quoted for the 
proper motion. 

Eleven of the stars have parallaxes exceeding 0”-1. Of these, nine are 
stars given by Luyten, and the other two are BD —21° 1051 (No. 1632) and 
BD +5° 1668 (No. 1650) a star with very large proper motion. 

The supplementary list contains twelve stars for which, for various reasons, 
second determinations of parallax were made. In each case a completely new 
series of plates was taken; but the reference stars and the ruled reference plates 
were the same as those used for the first determination made here. ‘The original 
determination is given under (1) and the second under (2). ‘The means are 
based upon the relative probable errors. 

NOTES 

1607. Dr Hoffleit (Harvard special classification for Cape Zone Photographic Catalogues, 
unpublished) classifies this star as Go. Its colour index, as determined at the Cape, 
Is + 0°25. 

1608. LPM "6 roa white dwarf. 

1668. This is CPD — 43° 7032, which is possibly CoD — 43° 9842. 

Mag. 

Vis. Type 

Phot. 


R.A., Dec. Proper Motion aaa ae ee 
(2900.0) ‘gpickacienMesitiaiate nit Para lax P.E. 
R.A. Dec. Weight 

hm 5 
LPM 32 ae © 35 oy 437 263 -29 + 46 
—44° 170 , —44° 47'°8 469 


Name and 


No. DM. No. 


LPM 49 ‘is < © 52 59 + 1050 170 tj 48 
62° 75 q — 62° 47°°0 + 1034 


LPM 50 ais © 53 28 1270 
-28 ° —28° 23:6 | 1299 


LPM 56 ins © 59 19 + 620 
—35° 360 —35°12"2 «+ 585 
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Mag. 2 Proper Motion P.E. 
4 ‘a " Vis. Type ae - tem, Unit Parallax 
en aes Phot. : R.A. Dec. Weight 
hm °s 
LPM 67 waa I 13 19 + 170 —690 +37 
—13° 24'°4 + 166 


LPM 73 ror 27 43 — 560 
—22° 526 . —22° 23'°7 - S97 


LPM 75 =e I 32 37 - 70 
—61° 122 . —61° 34'°4 — 99 


LPM 76 12° x 92°53 + §80 
; 5° 29°38 =+ «570 


Ci 20 126 ~ iia I 40 26 80 
Ross 552 . ~13° 59°°6 - $3 


GC 2286 7° ; 1 48 02 
— 2° 38 rh - 1° 48'°6 


GC 2303 . I 49 38 
ZC 890 . —46° 47'°5 


LPM 97 wi G4 2 00 50 
—28° 657 . — 28° 336 


LPM 100 ror k-m—s 2:«O4 21 
—17° 400 * —16° 48'-9 


LPM 104 ae 2 08 04 
—18° 09'-7 


LPM 107 _ 2 09 26 
—32° 828 "3 —32° 29'°6 


LPM 115 ‘ 2 14 00 
CoD —54 487 ° —54° 27'°6 


LPM 121 ‘ 2 24 19 
—20° 465 —20° 26’-0 
GC 3141 ° < 238 7 
— 3° 410 ve 3° 35°°6 


LPM 135 ie 2 50 O1 
—64° 06'°4 


GC 3734 j 3 02 35 
—14° 604 — 14° 08'-4 


LPM 144 aa 3 08 30 
— 38° 1058 . — 38° 290 


— 26° 1207 . ‘ 3 10 26 
~26° 49'°2 


LPM 159 a : 3 30 36 
—31° 1454 a: —31° 24°°8 


GC 4313 5 3 31 46 
10 Tauri 0° 051 
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nein anil Mag. et, R.A., Dec. Proper Motion Se : ™ : 
D.M. No Vis. l'ype (1900.0) _ Unit Parallax P.E. 
he thi Phot. ’ R.A. Dec, Weight 


hm °~.8 
1625 GC 4488 ‘96 ; 3 40 28 + 204 294 39 61 
38° 1264 | 38° 26": 225 


1626 LPM 173 ie : 3: 370 1080 
37° 1501 ’ : 409 


1627 FS 1228 hee ‘ 112 
50° 1302 a ' P 135 


1628 GC 5407 ° ; ao eae 


+3 


64° 334 66 


1629 LPM 198 wees 450 


391 


1630 ZC 2534 : F 30 
43° 1590 


1631 LPM 201 

1632 

1633 GC 6487 
PGC 1263 

1634 LPM 211 

1635 S.f. No. 1634 

1636 GC 6824 


10° 1204 


1637, LPM 218 
— 36° 2458 


1638 GC 7387 
25° 2734 


1639 6° 1339 

1640 9° 1261 

1641 GC 7479 
52° 805 


1642 GC 7738 
35° 2685 


1643p GC 7976 
— 65° 565 


1643f GC 7979 
— 65° 566 
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Seen ond Mag. Proper Motion P.E. 


Vis. 4 pansiicnti, Ctam Cie 
D.M. Ne. Phot. : R.A. Dec. Weight 


Mean 


LPM 234 as 6 17 45 
—22° 3005 ° —22° 40'°3 


L597-31 see 6 28 06 
—26° 57'1 


HH HH KH 


~26° 3112 iy see 6 29 00 
L. 276 , —26° 41°°6 


He Ht 


GC 8712 5° c 6 35 53 
PGC 1708 *s -30° 22'°4 


GC 9256 ‘ 6 57 08 
~25° 3913 ; -25° 48'°4 


rr ofr ee 


LPM 259 pa 7 16 38 
—15° 09'°4 


he 


5° 1668 . Ms 7 22 00 
+ §°32°3 


ZC 4651 7 ] 7 31 32 
—50° 2841 ; — 350" 37°*4 


GC 10425 5° 7 40 18 
PGC 2039 , —40° 41'°4 


HH HH HE 


GC 10524 . 7 44°55 
61° 888 . —61° 11'*4 


He H- 


LPM 282 sé 8 o8 18 
—21° 14"°0 

LPM 296 ee 8 36 39 
-23° 07'°7 


GC 13048 . ; 9 22 49 
PGC 2535 ° 05° 380 


LPM 321 se . 9 31 24 
—21° 13'°% 


LPM 334 jm 9 56 52 
29° 8019 ‘ 29° 56’-0 


GC 17025 8- 35 12 24 37 
LPM 432 - 02° 45°°9 


FS 10116 - . 12 29 30 
46° 8007 : 46 4471 


GC 17435 
°’ 2989 


GC 18144 
a Virginis 
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; Mag. Proper Motion P.E. 
<< Vis. Type yr = ate, Unit Poaiix PE. 
ane ‘ Phot. ‘ R.A. Dec. Weight 
hm °.8 
— 3° 3527 10°8 M213 39 55 11 — 64 +31 
oka —0o4° 069 149 


ZC 12041 ° j 14 02 29 331 +31 
—51° 7987 ; —§1° 27°°5 322 


LPM 532 ° 14 28 54 290 
Ci 20 870 a3 —12° 05':8 315 


LPM 544 oe 14 43 52 
—25° 40°°7 


LPM 543 vee 14 43 53 
~25° 10553 ‘ —25° 39''6 


Mean 


GC 20160 . 14 53 11 
ZC 12838 ° —43° 24'°5 


FS 13043 pve j 15 15 46 
—43° 303 

LPM 575 ee 1§ 32 13 
—20° 58'°5 


N.f. No. 1669 on os «6. 
—20° 56'°7 


S.p. No. 1672 ve on eam 
—13° 171 
LPM 585 x 15 41 35 
—13° 4246 : —13° 16’-9 
LPM 590 . 15 59 48 


—20° 4399 ° —20° 101 


LPM 597 shee : 16 13 30 
—37° 10765 , —37° 18'9 


LPM 598 wee 16 14 57 
Ci 20 980 . —17° 23':8 


LPM 617 fei 16 36 24 
—72° 48’ 

LPM 619 ove : 16 42 19 
—64° 155 


LPM 630 ane 17 07 31 
— 4 5971 


FS 16429 _ ; 17 15 18 
3 


43° 03°°3 


GC 23597 , . 17 20 16 
PGC 4420 —24° 05"°0 
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No. 


1681 


1682 


1683 


1684 


1685 


1686 


1687 


1688 


1689 


1690 


1691 


1692 


1693 


1694 


1695 


1696 


1697 


1698 


1699 


1700 


Mag. 
Vis. 
Phot. 


Name and 
D.M. No. 


LPM 661 


LPM 663 
—30° 15026 


LPM 671 


LPM 685 


GC 25941 
o Sagittarii 


FS 18390 
CAZ —47 


176 


GC 26161 
¢ Sagittari 


UOC 53 
Ci 20 1126 


Wolf 1062 
Ci 20 1136 


LPM 701 


FS 18826 


LPM 709 
59° 7535 


LPM 727 
20° 5833 


LPM 741 
—28° 16676 


LPM 746 
— 22° 5504 


LPM 748 


GC 29787 
0° 4195 


LPM 814 
Ci 20 1340 


GC 31165 
ZC 19833 


GC 31437 
ZC 19958 


Type 


G 


f-g 


R.A., Dec. 
(1900.0) 
im ®S 
17 50 12 
— 16° 22'°6 


17 52 27 
— 30° 08’-6 


18 13 58 
—05° 48-2 


18 44 24 
—ay i 9 


18 49 04 
26° 25'°3 


18 51 15 
— 46° 23'°9 


18 56 15 
— 30° O14 


18 59 03 
— 65° 244 


19 07 02 
- 2° 44°°2 


19 13 25 
-45° 39°°5 


19 38 04 
—40° 28'-6 


19 42 21 
—59° 25°°2 


20 04 33 
—20° 46'-6 


20 21 39 
— 28° 03:0 


20 35 47 
—22° 397 


20 37 18 
—79 40° 


21 12 52 
—00° 15'"1 


22 04 21 
—0§° 07''9 


22 It 22 
—49 08-8 


22 24 53 
— 49° 56'-4 


Proper Motion 
R.A. Dec. 


° —- 600 


os 


/ 


—650 


96 


- 420 
453 


220 
180 


620 
660 


670 
676 


450 
+ 476 


+ 1030 
+1139 


+ 190 
217 


- 
— 152 


Stellar parallaxes determined at the Cape Observatory 


P.E, 
Unit 
Weight 


£27 


Parallax 
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Second Determinations 


Name and rong Type RA» Dec. P.M. inp, Relative 
: , . —— Jaralle 
D.M. No. Phot. (1900.0) R.A, Parallax 
h m 8 
GC 1920 7°19 Gs I 30 23 140 i i (1) 


— 30° 529 8-o ~30° 25'°5 132 ‘ { (2) 
132 t b- Mean 


FS 1254 ies I 4 13 36 170 + + (1) 
CAZ — 41° 66 ° 41° 02''8 219 ~ (2) 
216 : 4 { Mean 


GC 9824 ‘ 7 17 09 358 t t (1) 
1° 1677 o1 40°°8 368 j _ 
367 t - § Mean 


Ci 18 918 7°25 < 7 34 59 : t 4 (1) 
— 3° 2001 . —0o3° 216 + 79 £ (2) 
78 t { Mean 


FS 10138 a Mo 12 31 21 88 : d (1) 
41° 7290 41° 098 ~ 86 + : { (2) 
86 t t+ 5 Mean 


FS 13662 ne ] 15 34 44 t t + (1) 
- 50° 9730 ~50° 34'°8 . t =o (2) 
t + 5 Mean 


FS 17279 ben ) 17 39 14 { t i (1) 
CAZ — 49° 362 . — 48° 03-0 t t } (2) 


Mean 


FS 18216 lar. 18 35 53 t t + (1) 
45° 12699 ee 45° 52°°5 i C { (2) 
t ‘ Mean 


FS 18799 wa J 19 35 44 $6 +6 i (1) 
45° 13383 -45° 18'5 + (2) 
t + 6 Mean 


FS 18800 be ; 19 35 50 L t ! (1) 
—45° 13385 45° 16'°7 t ; (2) 
t : Mean 


FS 19718 na ) 21 51 07 t t t (1) 
CAZ —51° 83 , 51° 28'*3 5 t t { (2) 
t } i Mean 


Royal Observatory, 
Cape of Good Hope : 
1951 October 15. 








THE PHOTO-CHEMISTRY OF SOME MINOR CONSTITUENTS 
OF THE EARTH’S ATMOSPHERE 
(CO,, CO, CH,, N,O) 


D. R. Bates and Agnes E. Witherspoon 
(Received 1951 October 3) 


Summary 


Carbon dioxide and carbon monoxide: The rate of photo-dissociation of 
CO, is appreciable only in the region above about 100 km. If local equi- 
librium prevails the carbon would occur mainly as the monoxide in this region, 
and as the dioxide below. However, oxidation may proceed so slowly that 
the life of a CO molecule is long compared with the characteristic time 
associated with atmospheric mixing effects. In this case the CO resulting 
from photo-dissociation would not be confined to above the 100 km level, 
but would extend much lower; its total abundance could scarcely be great 
enough to produce the observed absorption lines. It is estimated that the 
various forms of combustion taking place on the Earth would provide the CO 
content of the atmosphere within perhaps four years or even less. 

Methane : Dissociation of CH, in the upper atmosphere is brought about 
mainly by collision processes. ‘These prevent the existence of appreciable CH, 
above the 100 km level and probably keep the concentration low even down 
to70km. ‘The yield to date from oil wells and coal mines appears to be less 
than the amount of the gas now in the atmosphere. Seepage from fuel beds, 
and the anaerobic decay of vegetable matter, are the only obvious naturally 
occurring sources of significance. As far as can be judged at least ten years 
(and probably much longer) is required for the production of the atmospheric 
abundance. 

Nitrous oxide : Information on the necessary rate of formation of N,O is 
obtained by calculating the photo-dissociation rate. One hypothesis is that 
soil micro-organisms produce the gas, but it appears that the yield would only 
be sufficient if N,O were a major end-product of denitrification. Many of the 
homogeneous gas reactions suggested by earlier workers must be rejected, 
since there is an inadequate supply of the requisite parent particles. ‘The only 
acceptable parent particles seem to be those arising directly or indirectly from 
the photo-dissociation of O, in the Herzberg continuum. Various reactions 
involving O, O, and N, are discussed. ‘Though they are very slow they might 
nevertheless give rise to sufficient N,O. 





1. Introduction.._The photo-chemistry of atmospheric oxygen has been 
studied by many scientists (1), and though much detailed work remains to be 
done, the researches of Chapman and others have led to a general understanding 
of the factors involved in the formation of the ozone layer and of the atomic 
region. While the desirability of a similar study of nitrogen has long been 
realized, comparatively little has yet been achieved because of some special 
difficulties that arise and because of lack of certain basic data; and as nitrogen is 
obviously the next most important gas to investigate there has been a neglect of 
the other constituents—-though Meinel’s discovery (2) of the hydroxyl bands in 
the spectrum of the airglow recently stimulated Bates and Nicolet (3) to carry 
out an exploratory survey of the photo-chemistry of water vapour. The object 
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of the present paper is to remedy this omission partially by investigating the 
position regarding carbon dioxide, carbon monoxide, methane and _ nitrous 
oxide. Originally the research was directed at obtaining information on the 
altitude distribution of the polyatomic molecules (which is of importance in 
connection with studies of the thermo-balance in the upper atmosphere). During 
the course of it the existence of several interesting cycles was realized, and these 
are also discussed. 

2. Model atmosphere.—Table | gives the current estimates of the particle 
concentrations, n(M), and temperatures, 7, at various altitudes, z, up to 130km 
(3). They are significantly different from those accepted at the time of the main 
work on the O-O,—O, equilibrium. Since this equilibrium is of importance 
in the study of the minor constituents it had therefore to be re-investigated. 
‘The processes involved have been discussed fully by earlier writers (1), so that 
we need only indicate them briefly, and record the assumptions made in calculating 
the rates at which they occur. 

TABLE I 
Structure of Atmosphere 
Altitude Particle Concentration ‘Temperature 
n(M) zr 
(cm~*) (deg. K) 
10!” 450 
a 360 
x 101? 305 
XK r0"* 270 
K 16" 245 
x 104 215 
x 1918 190 
K fo” 275 
Ko" 290 
K 20** 250 
< 101? 225 
x 1018 220 


IN 


Sse 


I 

7 

2 
8: 
3°9 
2° 
g- 
2° 
6: 


x 1018 220 
< 101° 285 


Photo-dissociation of molecular oxygen is produced by the absorption of 
solar radiation in the strong Schumann-Runge continuum beginning at A1759: 
O, + Av OP) + OCD), (1) 
and in the weak Herzberg continuum beginning at A 2421: 
O, + Av > 20(°P). (2) 
In the presence of a third body the atoms thus liberated may recombine 
0+0+M-+0,+M, 
or they may unite with oxygen molecules 
0+0,+M+0,+M, 
forming ozone. ‘This substance may be destroyed by the collision process 
O+0;>20,, (5) 


and by photo-dissociation in the strong Hartley continuum beginning about 
A 3500: 


O,+ A4v> O,+ O('D), (0) 
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and in the weak Chappuis continuum beginning about A 8000: 
O,+hv - O, + O(PP). (7) 


Following the usage of other authors the reciprocals of the life-times of oxygen 
and ozone molecules towards photo-dissociation by (1) and (2), and by (5) and 
(6), will be denoted by J, and J;, and the rate coefficients associated with (3), 
(4) and (5) by &,, k, and k, respectively. 

‘The intensity of the incident radiation down to A2200 has been determined 
by means ot rocket-borne instruments {cf. Durant (4)}, and in the spectral region 
of interest beyond it probably approximates to that which would ensue if the Sun 
were a 5000deg.K black body {cf. Nicolet (§), Greenstein (6)}. Since the 
relevant absorption cross-sections are also known* the values of J, and J, at 
various optical depths can hence be computed. ‘The position regarding the 
rate coefficients is much less satisfactory. For reasons discussed elsewhere 
(3) it was tentatively assumed that in ordinary air 


,=5 x 10-471? cm®/sec, (8) 
2=5 x 107-72 cm® /sec, (9) 


where 7’, as usual, is the temperature; and that the relative efficiencies of O,, 
N, and O as third bodies are as 1-0:0-6:0-1 approximately {Penndorf (8)}. 
Since Eucken and Patat (g) have measured the ratio of k, to k, this latter is then 
also determined. 

At altitudes down to about 70 km it is usually supposed that the concentrations 
n(Q,) and n(Q) are controlled by processes (1), (2) and (3), so that 


J n(Ox) = k,n(O)?n(X), (10) 


where n(X) is the effective number of third bodies found by weighting O,, Ny 
and O in the manner indicated above. In addition, if m(M) is the total particle 
concentration (given in ‘Table I) then 


n(O,) + n(N,)+n(O)=n(M), (11) 


and since the relative abundance of oxygen and nitrogen is presumably the same 
as at ground level, 


n(O) + 2n(O,) = 0°53n(N,). (12) 


These equations are somewhat troublesome. They were solved by commencing 
at the 130km level (which is at essentially zero optical depth) and proceeding 
downwards to the 80 km level in altitude intervals so narrow that the attenuation 
in each is small {ef. Penndorf (8)}. ‘he labour involved was however scarcely 
justified, since eqn. (10) over-simplifies the situation. Atomic hydrogen and 
certain radicals, which are present i: the upper atmosphere, catalyse the formation 
of molecular oxygen. Bates and Nicolet (3) have investigated the extent to which 
this is likely to influence the equilibrium. They found that the effect may be 
significant below &85km, and gave crude estimates of m(Q) and n(Q,) down to 
the 65km level. ‘These were adopted, as there are insufficient basic data to 
improve upon them. At low altitudes catalysis becomes unimportant and 
(2), (4), (5), (6) and (7) are the dominant reactionse In principle the equilibrium 

* Reference (3) lists the principal original papers on the subject and discusses the extrapolation 
of some of the experimental results. Its recommendations were followed, except that in the case 


of the Herzberg continuum the absorption cross-section curve obtained by Buisson, Jausseran and 
Rouard (7) was adopted. 
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Fic. 1.— Distribution of n(O) and n(O,). (Day equilibrium.*) 





Ordinates : Altitude in kilometres. 
Abscissae : Log, (particle concentration). 
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—Distribution of n(O) and n(O,) and n(O;). (Day equilibrium.) 


Ordinates : Altitude in kilometres. 
Abscissae : Logy (particle concentration). 


* In the region covered by Fig. 1, #(O;) is too small to be of interest. (It should be noted that 
the moderate rate of fall-off shown in Fig. 2 does not continue to greater altitudes, it becomes much 
more rapid at the O,-O transition layer.) 
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can be determined using theory alone, but various difficulties enter and the full 
calculations were carried out only near the 50km level. For altitudes up to 
40km it was judged best to accept the measurements on m(QO,) {Fabry (7)} and 
to deduce n(O) from them by means of the relation 


n(O) =n(O,){ J, kan(O.)n(M)} (13) 
with 
n(O,) =0-21n(M) (14) 


{ef. (3)}. As has already been mentioned &, is not known reliably. It is possible 
that steric hindrance is not as severe as assumed {cf. Bamtord(10)}, so that the 
use of (13) may cause (QO) to be overestimated. 

The distributions in the different regions were finally combined together as 
shown in Figs. 1 and 2. It may perhaps be mentioned that the main changes in 
the slopes of the n(O) and n(O,) curves are due to the fact that the coefficient 
k, is extremely sensitive to the temperature. While such qualitative features 
are likely to be real the general accuracy of the results is unfortunately very poor. 

Some of the reactions of interest involve atomic hydrogen, hydroxyl and 
perhydroxyl. Estimates of n(H), 2(OQH) and n(HO.,) are available from the work 
of Bates and Nicolet (3). 

3. Carbon dioxide and carbon monoxide.— lhe fractional volume concentration 
of carbon dioxide is not the same all over the Earth’s surface, even apart from 
purely local effects {cf. Rankama and Sahama (1x1)}. However, for present 
purposes it may with sufficient accuracy be taken as constant and equal to 
310 {Paneth (12)}. In the altitude range concerned diffusive separation 
is unimportant. 

In 1949 Migeotte (13), then working at Columbus, Ohio, announced that 
certain lines in the 4-7 region of the solar spectrum are due to atmospheric 
carbon monoxide. Later in Switzerland he and Neven (14) succeeded in 
confirming this. Adel (1§), in Flagstaff, Arizona, failed to detect the lines, but 
as they appear in a spectrum recently obtained by Shaw, Chapman, Howard and 
Oxholm (16) (also working at Columbus, Ohio), their reality is not in doubt. 
The lines observed are members of the fundamental vibrational band. Matheson, 
and Penner and Weber (17) have carried out laboratory measurements on them. 
The former found the Lorentz half-breadth in air at atmospheric pressure to 
be o-10cm ', and the latter found the associated Einstein spontaneous transition 
coefficient to be 33sec. Combining these data with equivalent widths obtained 
from the published spectrograms (in particular from those of Shaw and his 
collaborators, as they include some weak lines attributed to the isotopic molecule 
C0 which probably lie on the linear part of the curve of growth), it is estimated 
that there is about 0-2 atm-cm of cardon monoxide in the atmosphere. ‘The 
figure is of necessity rather crude—an error by a factor of perhaps 3 either way 
is possible. Much more reliable determinations have recently been made by 
Migeotte and Neven*, who on two different days obtained about 0-4 atm-cm 
and 0-1 atm-cm at the Jungfraujoch (showing that the abundance is variable). 
The lower value (which is presumably the least affected by any local contamina- 
tion) will be tentatively adopted; it corresponds approximately to a fractional 
volume concentration of 1 x 10 7 or to 3 x 10!8 CO molecules cm? column. 


* We wish to thank Dr Migeotte for informing us of the results prior to their publication. 
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‘The photo-chemistry of the two gases is clearly closely related. Radiation to 
the short wave-length side of A1690 dissociates carbon dioxide into carbon 
monoxide and excited ('D) atomic oxygen 


CO, +hv->CO+ O('D). (15) 


‘The cross-section for the process has been measured by Wilkinson and Johnston 
(18). Using their results, taking the Sun as usual to be a 5 000deg. K black 
body (cf. Section 2) and averaging over day and night, it may be shown that in 
the region just above the 100km level 7(CO,) the life-time of a CO, molecule 
towards photo-dissociation through (15) is about 5 x 10°sec. Owing to absorption 
in the Schumann—Runge continuum of molecular oxygen the intensity of the 
active radiation falls off rapidly below this level. Consequently 7(CO,) becomes 
extremely long, being very many years even at 95km, and essentially infinite 
at 90km. 
Carbon dioxide may be re-formed through 


CO+0+M+CO,+M. (16) 


This process has been studied by Groth (19) and by Jackson (20), who find that 
the probability that it occurs in a normal three-body collision is some 100 times 
less than the probability of (4), which itself is small. If the estimate given in (9) 
is correct, then the rate coefficient for (16) would only be 5 x 10-87"? cm®/sec. 
As is well known (cf. 21), the oxidation of carbon monoxide is accelerated by the 
presence of moisture. ‘This is generally attributed to the reactions 


CO+OH—>CO,+H, (17) 
CO+ HO, > CO,+OH. (18) 


These must also be taken into account in the upper atmosphere. ‘The rate 
coefficients may be expressed in the form 


k=nPexp(-—E/RT), (19) 


where R is the gas constant (1-987 cal/deg-mole), is the collision frequency for 
unit concentrations, P is the steric hindrance factor and E is the activation 
energy. Unless there is evidence to the contrary it is usual in order-of-magnitude 
calculations to take the product »P to be about 1°5 x 10 ''7"? cm/sec or 
rather less (cf. 22). At the temperatures prevailing in the upper atmospheric 
region of interest the rate coefficients depend critically on the activation energy. 
Unfortunately the values appropriate to (17) and (18) do not appear to have been 
determined. ‘To explore the situation two alternative arbitrary assumptions 
were made: (a) that for both reactions E is zero, (6) that for both it is 5 kcal. 
Calculations were then carried out on 7(CQO), the mean life-time of a CO molecule 
towards oxidation by (16), (17) and (18) in the model atmosphere already described. 
In spite of the uncertainties some general conclusions can be reached. 

Even with the highest possible rate coefficients it is found that 7(CO) is of 
the order of centuries near the 110km level and is longer at greater altitudes. 
Life-times so extended have of course little real meaning. Below 100 km however 
7(CO) may be moderate: at the go, 80 and 70km levels assumption (a) gives it 
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to be about 105, 10% and 10sec and assumption (b) gives it to be about 10%, 10% 
and 10’ sec respectively. ‘Thus 

7(CO,) <7(CO) for z>110km (20) 
and 7(CO,) >7(CO) for z<90km. (21) 
In equilibrium 

n(CO)/n(CO,) = 7(CO)/7(CO,), (22) 
and therefore the carbon would be almost entirely in the dioxide form below 
the 90 km level and in the monoxide* form above the 110 km level. ‘The altitude 
of the transition level is difficult to estimate precisely but it is certainly very close 
to 100km. For reasons which will be apparent later the significance to be attached 
to the derived relative abundances does not justify giving them in full. To 
demonstrate their trend it is sufficient to mention that at 130km (where (16) is 
the only oxidation process of importance) the calculated value of (CO) n(CO,) 
is 10° approximately. ‘ 

It is natural to enquire at this stage whether the photo-dissociation would 
give rise ‘to enough carbon monoxide to account for the observed infra-red 
absorption lines. At first it might be thought obvious that this is utterly 
impossible ; for the fractional volume concentration of carbon dioxide in the 
lower atmosphere is 310 4, and the total particle concentration and local 
scale-height at 100 km are 2-3 x 10'* cm* and 10km respectively. Hence there 
would only be about 7 x 10'!° CO molecules;cm? column above the transition 
level, which is much too small a number. However, such an argument is faulty, 
in that it assumes that equilibrium conditions prevail. ‘This assumption is 
only justified if 7(CO,) and 7(CO) are short compared with a characteristic 
atmospheric parameter, 7(/), which is a measure of the time required for winds 
and diffusion to produce mixing. ‘The value of 7(M) is unfortunately unknown, 
but it must be considerable in view of the stabilizing effect of the positive tempera- 
ture gradient of the atmosphere in the region of the transition layer. As has 
already been noted, 7(CO,) is about 5 x 10° sec at altitudes above tookm. This 
is not unduly long, so that, in general, the degree of dissociation may well be 
high (though perhaps not as high as simple equilibrium calculations would 
indicate). At altitudes below 1tookm, 7(CQO) is the life-time of importance. 
Its value is here mainly controlled by processes (17) and (18). Or, referring to 
the figures already given, it can be that if the activation energy E is zero then 
7(CO) would be quite short. In this case the mean concentration of carbon 
monoxide would be kept low by oxidation. However, if, as is more probable, 
E is 5 kcal or above, the position would be quite different, since 7(CO) would be 
several months or more even at 70km. Consequently the equilibrium calcula- 
tions would be unlikely to be correct. Instead of carbon monoxide being mainly 
confined to the region above the transition level it would also be prevalent in the 
region below, for any reaching there would remain for a very long time before 
becoming oxidized. In extreme circumstances the altitude at which photo- 
dissociation ceases would not be apparent from the carbon monoxide distribution 
curve, which would closely follow that of the main atmospheric gases down to 
perhaps 70km or even lower. The total abundance resulting from photo- 
dissociation cannot be calculated until 7(CO) and 7(M) are known. ‘Though the 
departure from chemical equilibrium may well be very pronounced, the associated 


* Further breakdown of the carbon oxides is unlikely to occur to any marked extent, since CO 
(like the isoelectronic N,) is difficult to photo-dissociate {cf. Gaydon (23)}. 
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increase in the abundance is unlikely to be sufficient to explain the spectroscopic 
observations; for this would only be possible if the monoxide were as prevalent 
as the dioxide down to the 30 or 40 km level. A measurement of the altitude 
of the region in which the CO lines are formed would be of considerable interest. 

It must not be forgotten that carbon monoxide is a major by-product of 
modern civilization. Vast quantities are vented into the atmosphere both by 
internal combustion engines and by solid-fuel furnaces. 

The first source is the easier to assess quantitatively. Petrol consists of about 
85 parts of carbon to 15 parts of other elements (24). ‘The annual consumption 
of it is about 1-5 x 10'g. From an analysis of the exhaust gases of motor 
vehicles given by Dunstan (25) it appears that as much as 40 per cent of the carbon 
is emitted as the monoxide. Hence the amount of carbon monoxide released 
into the atmosphere is approximately 1-2 x 10! g/year. Since the surface area 
of the Earth is 5-1 x 10'* cm?, and the mass of a CO molecule is 4-7 x 10°*8 g, this 
is equivalent to a mean of 5 x 10'? CO molecules/cm?/year. Diesel engines give 
but a minor contribution, as they are highly efficient and as the quantity of heavy 
oil used is relatively small. 

Coal is by far the most important solid fuel. Its composition varies widely. 
The carbon content may be as low as 60 per cent in lignites and as high as 95 per 
cent in anthracites (26). On an average it may be taken to be perhaps 80 per cent. 
‘The amount of coal used in a year is about 1-7 x 10! g. Unfortunately it is difficult 
to obtain representative data on the proportion of carbon monoxide in the 
combustion products, as this is critically affected by the conditions under which 
burning takes place. ‘hus Hurley and Sparks (27) have analysed the flue 
gases from boiler furnaces and find that the CO/CO, ratio ranges from 0-005 
to 0:20, depending on such factors as the air draft. Fuel engineers of course try 
to keep the ratio as low as possible in the interests of thermal efficiency. When it 
equals 0-02 the heat loss is only 1 per cent (26), which suggests that little effort 
would be made to ensure a more complete combustion (especially as the use of too 
much excess air also causes losses). Forest and prairie fires are estimated by Kalle 
(cf. 11) to add as much carbon dioxide to the atmosphere as does the burning of 
coal, but no information is available on the carbon monoxide. From the figures 
given above it would seem not unreasonable to suppose that the mean value of 
the CO/CO, ratio for all forms of non-explosive combustion is at least 0-01 and 
is probably rather more. On this assumption the contribution to the production 
rate (averaged over the Earth) is 3 x 10'7 CO molecules/cm?/year or greater. 

Volcanic emanations contain some carbon monoxide, but from data given by 
Rankama and Sahama (11) it is at once apparent that the quantity released in this 
way is small. 

Summation of the sources mentioned gives that the total production rate is 
at least some & x 10!’ CO molecules/cm*/year, which is sufficient to provide the 
estimated carbon monoxide content of the atmosphere within four years or less. 
Unless some destruction process occurs on the ground itself or at low altitudes, 
such a production rate must be of importance in determining the equilibrium 
abundance. For comparison it may be mentioned that the action of solar 
radiation on the upper atmosphere would yield fewer than 1 x 10/6 CO mole- 
cules/cm® year even if the fractional volume concentration of carbon dioxide in 
the region concerned were kept as high as one-fifth of what it would be if photo- 
dissociation did not occur. 





No. 1, 1952 Photo-chemistry of minor constituents of the Earth’s atmosphere — Log 


Combustion processes certainly lead to considerable localized contamination : 
for example, while the presumed global average of the fractional volume concen- 
tration of carbon monoxide is I x 10~?, it appears from measurements made by 
various investigators that the value in busy streets is frequently 1 x 10-4, that 
in cities 5 x 10°* and that in urban areas 2 x 1077 or less (28).* Clearly winds 
would ensure a widespread dispersal of the contamination even if the mean life 
of the CO molecules (which has yet to be discussed) should prove to be only 
a few days. In spite of its remoteness the Bernese Oberland (where Migeotte 
and Neven conducted their experiments) would be affected. On the other 
hand, the fact that the combustion occurring in the northern hemisphere is 
far greater than that occurring in the southern hemisphere might conceivably 
cause carbon monoxide to be more abundant in the former than in the latter. 
‘Thus the meridional circulation of the atmosphere (30) is not such as to 
facilitate the transfer of contamination from one hemisphere to the other; and 
in any case even a steady breeze having a N.-S. component of 20km/hour 
would take about a month to travel go” of latitude. It would be interesting in 
this connection to have data on the amount of carbon monoxide above some 
southern observatory (say in New Zealand). 

The carbon monoxide molecules resulting from combustion apparently 
cannot have a life exceeding four years (v. supra). It would appear therefore that 
oxidation processes occurring in the upper atmosphere can be ignored. Thus the 
particle concentration at the 50 km level is more than 10° times smaller than that 
at ground level. Hence if a molecule were normally to reach the 50km level 
within the time stated it would follow that a molecule at this level would 
normally take only a day to fall to ground level, which seems unlikely. 

It is difficult to reach any final conclusion about the rate of oxidation of carbon 
monoxide in the dense lower atmosphere. For a process might restrict r(CO) to 
less than four years and yet be soslow as to have escaped detection in the laboratory. 
Some relevant information is however available. 

The reaction between carbon monoxide and molecular oxygen is exothermic, 
but it can probably be dismissed, since Mellor (31) reports that Berthelot could 
find no evidence for its occurrence when he exposed dry and moist mixtures of 
the two gases to sunlight for seven years. As can readily be verified from the 
estimated rate coefficient of (16), and the concentrations given in Fig. 2, atomic 
oxygen cannot cause appreciable loss. Ozone also is a very inefficient oxidizer of 
carbon monoxide. Zatsiorskii, Kondrateev and Solnishkova (32) find that the 
activation energy of 

O,+CO-O,+CO, (23) 


is about 20kcal. On substituting in (19) it can hence be seen that the process 
is unimportant. A prohibitively high activation energy appears indeed to be 
characteristic of the reaction with many oxygen compounds. For example, 
that associated with 

NO, + CO-+ NO+CO, (24) 





* Abnormally high concentrations of carbon dioxide of course also occur but the departure 
from the mean is less pronounced. ‘Thus Haldane (29) states the fractional volume concentration 
in the air of London is about 6°5 ~ 10°* which is little over twice the global average of 3 x 107‘. 
Judging from the ratio of the excess CO to excess CO, in town air it would seem unlikely that we 
have overestimated the rate of production of carbon monoxide. 

In view of the spectroscopic estimates, the figure for urban areas is surprisingly small. 
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is thought to be 28 kcal {Brown and Crist (33)}. However, the atmosphere may 
contain minute traces of some highly active substances and it is not easy to exclude 
the possibility of one of them oxidizing carbon monoxide through a chain such as 
CO+X+M-->-COX+ M, 
COX + O,-> CO, + OX, 
| OX+hv+O+X, 
| OX+CO>CO, +X. 
Quite a moderate collision efficiency would suffice even if the abundance of the 
catalyst at low levels were only 10- times that of all other constituents. 

Reactions on the surface of the Earth may also occur. (i) Carbon monoxide 
unites readily with haemoglobin. ‘The process is however reversible and in fact 
very little of the gas is oxidized in the body tissues. Clark, Stannard and Fenn 
(34) have recently investigated the effect. From their results it may easily be 
demonstrated that the global average of the rate of oxidation is negligible. (ii) 
Plants are affected by the presence of carbon monoxide (35) which may therefore 
be undergoing destruction by them (though the chlorophyll reaction may also be 
reversible). Quantitative data appear to be lacking. But noting that the life-time 
of a CO, molecule towards assimilation is about ten years {cf. Rankama and 
Sahama (11)} it would seem rather improbable (though not impossible) that 
the corresponding life-time of a CO molecule is four years or less. (iii) As has 
long been known, carbon monoxide in contact with soil is rapidly oxidized by 
various micro-organisms, in particular B. oligocarbophilus (36). The effect 
has been investigated by a number of workers but it is difficult to make quantitative 
deductions from their results. Most of the experiments consisted of introducing 
an air—carbon monoxide mixture into a flask containing some form of organic 
matter and studying how its composition changed with time. Usually there was 
little change during an initial period, but after one or two weeks (during which the 
micro-organisms presumably multiplied) oxidation began, and after another one 
or two weeks was complete. The interpretation of all such laboratory studies 
is of course complicated by the fact that both the carbon monoxide content of 
the air used and the surface to volume ratio of the apparatus had values far in 
excess of those typical of the lower atmosphere. Nevertheless it is clear from them 
that B. oligocarbophilus is a very powerful agency for the destruction of carbon 
monoxide. Some suggestive field observations have also been reported. ‘Thus 
after an extensive fire in an anthracite mine Jones and Scott (37) carried out 
periodic analysis of the gases in some sealed galleries of volume 2 x 10! cm’. 
They found that the carbon monoxide disappeared within 20 days and attributed 
this to the oxidizing action of micro-organisms. 

From the foregoing discussion it is apparent that the mechanism responsible 
for the removal of carbon monoxide from the air near ground level cannot yet be 
identified with certainty; and that only a lower limit can be given to the rate of 
the process concerned. In the circumstances no final decision can be reached 
as to the extent to which combustion influences the equilibrium. But so much 
of the gas is produced in this way that considerable quantities must get carried 
into the stratosphere in spite of the possible rapid destruction near ground level. 
Bearing in mind that the yield from the photo-dissociation of carbon dioxide is 
relatively very small, it would seem likely that combustion is the main source of 
the carbon monoxide observed spectroscopically. 
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4. Methane.—The only hydrocarbon so far detected in the Earth’s atmosphere 
is methane.* Migeotte (40) identified some of its characteristic lines in the 
3°4 region of a solar spectrogram taken in Columbus, Ohio, and later (41) 
identified further lines in the 7-7. region of a spectrogram obtained by Adel in 
Flagstaff, Arizona. ‘The discovery was confirmed by McMath, Mohler and 
Goldberg (42) working in Pontiac, Michigan, on the 2-0, region. Later 
measurements at the same observatory (43) give the abundance as 1:2 atm-cm 
or 3:2 x 10'* CH, molecules/cm? column. 

Photo-dissociation of methane, 


CH, + hv > CH, +H, (26) 


begins about A1450, and reaches its peak between A1300 and A1200. The 
absorption cross-sections in the region up to A1365 have been determined by 
Wilkinson and Johnston (18). ‘Those in the region beyond are relatively 
unimportant because of the fall-off in the intensity of the solar emission. ‘They 
were estimated with sufficient accuracy by extrapolation using as a guide crude 
data provided by the work of Duncan and Howe (44). On carrying out the 
usual computations the rate of (26) at zero optical depth was found to be only 
about 4 x 10°7/sec per CH, molecule. It will be seen later that photo-dissociation 
is not the major disruptive process. 
Recombination of methyl radicals and hydrogen atoms requires the presence 
of a third body 
CH,+H+X~+>CH,+X (27) 


{Steacie (45)} and consequently is very slow at the low gas densities prevailing 
in the upper atmosphere. Indeed a collision between a methyl radical and a 
hydrogen atom is more likely to lead to further break-down through 


CH,+H-> CH, + H,. (28) 


For this is a two-body process and its activation energy is probably less than 
5 kcal {Gorin, Kauzmann, Walter and Eyring (46)}. Moreover the methyl 
radicals are destroyed by numerous other collision reactions {cf. (45)} which 
convert them finally to carbon monoxide or carbon dioxide. Amongst the 
primary processes likely to be important is 
CH, + O, +CH,0+ OH. (29) 
Van Tiggelen (47) has given evidence for supposing that this proceeds very 
rapidly (the estimated activation energy being only some 1:5kcal). It is apparent 
therefore that the removal of a single hydrogen atom from a methane molecule 
leads to the essentially permanent destruction of that molecule. ‘The concen- 
tration of methane at eny particular altitude is thus not determined by the local 
equilibrium between dissociation and recombination, but rather by that between 
dissociation and atmospheric mixing effects. As a result the distribution of the 
gas cannot be calculated at present. But some general information on it may be 
obtained from the consideration of 7(CH,), the mean life-time of methane 
molecules at various levels. 
In addition to the direct photo action (already discussed) 7(CH,) is influenced 
by several collision processes. Methane reacts rather slowly with atomic oxygen, 
* The suggestion that ethylene is present {Sutherland and Callendar (38)} is not supported by 


the recent observations of Shaw and Claassen (39) in Columbus, Ohio, which were sufficiently 
sensitive to have revealed much less than o-o1 atm-cm of the gas. 
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atomic hydrogen and hydroxyl. ‘Thus the activation energy of the more probable 
of 

CH,+0O-CH,+H,O (30) 
and 

CH,+O-CH,0+ H, (31) 
is about 8 kcal {Steacie and Parlee (48)}; that of 

CH,+H-~>CH,+H, (32) 

is about 11 kcal {Geib and Steacie (49)}; and that of 
: CH,+OH—> CH, + H,O (33) 
is also about 8kcal {van Tiggelen (47)}. However, Geib and Harteck (50) 
tind that when perhydroxy]l radicals are present methane is rapidly oxidized even 
at temperatures below 100 deg. K, so that with these there must be a reaction 
with a very low activation energy (probably not exceeding 2 kcal). 

The activation energies indicated were substituted in (19) to obtain crude 
estimates of the rate coefficients. ‘These estimates were used in the compilation 
of Table II, which gives the approximate values of 7(CH,) at various altitudes, 
and records which of the reactants is the most important. 


Tase II 
Oxidation of Methane by Collision Processes* 


Altitude 7(CH,) 

(km) (sec) 

120 oO 10? 

100 O 103 

go O, HO, 10° 

80 HO, 10° 

70 HO, 104 
* In the region below 70 km the calculated values of 7(CH,) are too unreliable to be 
worth giving. One major source of error is the uncertainty in the perhydroxyl concentration 
(3). This concentration probably decreases at altitudes much less than 70 km so that 
7(CH,) increases. However, because of the rise in temperature the atomic oxygen reaction 
again becomes mmportant and between 60km and 50 km it limits 7(CH,) to perhaps 

10°-10" sec. 


Principal Reactant 


Clearly methane must be almost completely absent from the region above the 
100km level, since 7(CH,) is there very short. Remembering that a thermo- 
incline hinders mixing, it seems probable indeed that oxidation keeps its abundance 
low, down to at least near the temperature minimum at 70km. Little can be 
said about the equilibrium in the thermo-decline immediately below, but it should 
be noted that there is likely to be a significant loss of methane between the 60 km 
and 50km levels where the temperature is a maximum. Owing to the nature 
of the reaction products the destruction of CH, molecules at such moderate 
altitudes may be a very significant factor in the chemistry of the upper atmosphere. 

In view of the immense amount of fuel won from the Earth’s crust during 
the past century it is natural to enquire if the methane now in the atmosphere 
could have escaped from oil wells and coal mines and thus be of comparatively 
recent origin.* ‘The following figures are of relevance. 

* We wish to thank Mr P. Kerr of the Institute of Petroleum, London, and Mr F. S. Hartwell 


of the Safety in Mines Research and Testing Establishment, Sheffield, for helpful correspondence 
on this topic. 
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Volume of atmospheric methane 
(v. supra) =6 x 10!8 cm? (34) 


[ Amount of crude oil raised to date =6 x 10" barrels (35) 


Average volume of methane associated 
with one barrel of crude oil =7 x 10’ cm*® (36) 


Volume of methane associated with 
total crude oil raised =4 x 10/8 cm% (37) 


_ Amount of coal raised to date =I x 10" tons (38) 





Average volume of methane associated 
with one ton of coal =5 x 10°cm*® (39) 


Volume of methane associated with 
total coal raised =5 x 10!7cm? (40) 


The main uncertainty clearly lies in the assumed quantitative relations between 
the methane and the crude oil and coal. ‘Those quoted are based on information 
given in the standard treatise of Bone and Himus (26). Confirmation of relation- 
ship (36) was obtained from published statistical data on the production of 
natural gas (which is of course mainly methane) and crude oil {cf. ‘Thorp (28)}. 
High accuracy is not required, since the methane released by oil-field operations 
is normally not deliberately vented into the atmosphere because of the risk of an 
explosion. Some is pumped back underground (to increase the oil output) ; 
some is used for local power generation ; the remainder is either piped away from 
the field for general industrial or household purposes, or (in undeveloped countries 
where this is not practicable) is simply burnt in flares. Consequently the volume 
of methane which has actually escaped must be much less than 4 x 10'%cm!, 
As almost all the methane released in mining operations passes into the atmosphere 
relationship (39) is rather critical. Fortunately most authorities seem to make 
much the same estimate as do Bone and Himus. ‘Thus by averaging the methane 
content of the upcast air of 500 pits Hinsley (§1) deduces that the annual discharge 
of the gas from British collieries is some 2 x 10° cm®, which corresponds to I x 107 
cm*/ton* ; and this is consistent with the generally accepted ventilation require- 
ment {cf. Statham (§2)}. To be sure, many of the mines described in the 
literature are much richer in methane; for example, in the Lancashire coalfield 
some have yields of about 1 x 10°cm/ton {Fraser (53)} and in the N. Stafford- 
shire coalfield some have yields of almost 2 x 10° cm*/ton {Jones (§4)}. But such 
cases are reported because of the exceptional fire-damp hazard and should not be 
taken as representative. 

The conclusion would hence appear to be that though the quantity of methane 
which has escaped from oil wells and coal mines is vast it is nevertheless probably 
smaller than the amount now in the atmosphere. 

There are only two obvious naturally occurring sources of methane: seepage 
from fuel beds and vegetable matter decaying under anaerobic conditions. No 
direct quantitative information on the yield could be found but a limit to it can 
at least be set. 

* Since this was written we have received an almost identical (but slightly smaller) estimate 


from Mr G. H. Watson of the National Coal Board, London, who also informs us that the figure - 
for continental Europe is about the same, and that for America is probably less. 
8 
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Assimilation deprives the atmosphere of about 4 x 10?° C-atoms/cm?/year 
{cf. Rankama and Sahama (11)}. ‘The gain of CH, molecules must certainly 
be far less: that it should be more than 0-01 times as much is indeed scarcely 
credible. Consequently the time required to produce the observed 3-2 x 10!” 
CH, molecules/cm* column must be quite long—probably well in excess of ten 
years. Far too little is known about the loss processes in the lower atmosphere 
to make a discussion of them profitable at present, but it may be remarked that 
methane (like carbon monoxide) is oxidized by certain chemically active sub- 
stances and by various micro-organisms (36). 

Slobad and Krogh (§5§) have carried out some mass spectrographic analyses 
of the air near ground level. ‘They found the fractional volume of hydrocarbons 
to be 5 x 10°* tor x 107’, which is some 40 to 20 times smaller than the fractional 
volume concentration methane would have if it were uniformly distributed 
through the atmosphere.* ‘The natural inference is that the gas is largely confined 
to a layer at moderate or high altitudes. Migeotte (56), too, claims that the 
absorption lines are formed in a region where the temperature is 236 deg. K. 
‘This temperature is encountered at only four levels (cf. Table I). Only the two 
lower, one at about 8km and another at about 35 km, need be considered, since 
any methane formed near the two upper would not have a sufficiently prolonged 
existence (v. supra), and since in any case if it were confined to a layer near these 
upper levels its fractional volume concentration would have to be improbably 
high in order to provide the observed total content in a vertical column. 

Opposing evidence on the altitude distribution has recently been presented 
by Goldberg (57), who considers that it is essentially the same as that of the main 
atmospheric constituents. He reached this conclusion from a study of curves 
of growth of the absorption lines obtained at Mt. Wilson, California, and Lake 
Angelus, Michigan, and from a comparison of the total abundance of methane 
over these two stations (which are at altitudes of 1-7 km and 0-3 km respectively). 
In addition he pointed out that the true temperature of the methane is appreciably 
higher than the effective temperature obtained from the spectroscopic measure- 
ments, and on applying the necessary correction deduced a value only slightly 
lower than that at ground level. ‘This disposes of one of the original arguments 
for a methane layer at moderate or high altitude. But the positive evidence on the 
methane distribution is hardly final, since the curve of growth test is not yet very 
sensitive and since (as Goldberg himself realized) the comparison of abundances 
test involves the implicit assumption that the methane content of the air depends 
only on the altitude and not on the geographical position.t There is clearly a 
grievous conflict with the measurements of Slobad and Krogh. A check on 
these is urgently required. As data up to 25km would be useful, the facilities 
offered by balloon flights should be exploited. ‘The altitude of the centre of 
mass of the methane is of considerable interest, since it determines the fractional 
volume concentration in the mesosphere, and hence the importance of the gas 
as a source of hydrogen (3). 


* In view of possibility of the escape of hydrocarbons from oil fields it is interesting to note 
that the air analysed included samples from different parts of Texas. It is desirable that the 
investigations be extended to cover tropical regions (where anaerobic decay is particularly important). 

t To emphasize the importance of this assumption it may be noted that there is only a 20 per cent 
difference in the observed abundances over Lake Angelus and over Mt. Wilson. Methane 
production by decay would be greater near the former station than near the latter—though whether 
\t is appreciable in either is uncertain. 





No. 1, 1952 Photo-chemistry of minor constituents of the Earth's atmosphere 115 


A layer such as that apparently indicated by the results of Slobad and Krogh 
would scarcely be anticipated. It is of course conceivable that the main pro- 
duction takes place in the tropics and that the methane is carried rapidly to the 
stratosphere by the rising air currents there prevalent {cf. Rossby (30)}; further- 
more, the lower atmosphere is probably being depleted of methane (for example, 
through destruction on the Earth’s surface caused by micro-organisms). But if 
the characteristic time parameter associated with the gas is ten years or more, as 
estimated, it would be expected that mixing effects would be sufficient to keep 
the distribution fairly uniform. Goldberg’s results are thus favoured. 

5. Nitrous oxide.—The first evidence for the presence of nitrous oxide in 
the atmosphere was obtained by Adel (58) from observations on the 7-6 region 
of the solar spectrum. Confirmatory evidence has since been obtained by other 
infra-red spectroscopists (§9), and the abundance has been estimated to be 
o-8atm-cm or 2 x 10!® molecules/cm? column. ‘The location has not yet been 
properly determined. However, if the nitrous oxide were mainly confined to a 
10km thick layer* the mean concentration would be as much as 2 x 10!3/em?; 
and since the total particle concentration at an altitude of 7okm is only 
2°5 x 10!°/cm* it would seem most unlikely that any such layer could lie so high. 
More definite evidence will be reported later. 

The long wave-length limits of the nitrous oxide continua are rather ill-defined, 
but it appears that a very weak absorption begins at about A 3070 and that the main 
absorption begins at about A 2400 or slightly above. ‘The products of dissociation 
have not yet been identified with certainty but there is little doubt that in the low- 
energy region the process taking place is 

N,O+h/v>N,+ O(?P or 'D), (41) 
and there is considerable support for the view that that in the high-energy region is 
N,O + hv > NO+ N(4S) (42) 
(cf. Bamford (10), Gaydon (23)}. Dutta, and Romand and Mayence (60), 
have determined the absorption cross-sections to the short wave-length side of 
42700. ‘Their results are in harmony. Unfortunately quantitative measure- 
ments have not been made in the remaining spectral region of interest. A 
qualitative study has however been carried out by Sponer and Bonner (61), who 
find that the absorption diminishes as the wave-length is increased up to about 
2820, and that it then rises to a broad maximum near A29g00. ‘The variation 
does not appear to be very pronounced between A 2800 and the limit at 43070; 
and from an estimate given by Sponer and Bonner the cross-section in this region 
appears to be of order 1 x 10°*8cm?. Sufficient data are thus available to calculate 
the approximate photo-dissociation rate at various altitudes. ‘The values 
obtained are displayed in Table III. Those at very low levels may be consider- 
ably in error because of uncertainty regarding the extent of the attenuation of 
the incoming solar radiation. 

The destruction of nitrous oxide by collision processes is relatively unim- 

portant. ‘Thus 
N,0+0O-+N,+0, (43) 


and N,O0+O0-+2NO (44) 





* jo km is of course the approximate scale height in the lower and middle atmosphere. 
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can be neglected, since the activation energy of both exceeds 13 kcal {Henrique, 
Duncan and Noyes (62)}; and the experiments of Willey and Rideal (63) on the 
chemical action of active nitrogen on various substances indicate that 

N,O+N>N,+NO (45) 
is also slow. 

As it is by no means obvious that there is any upper atmospheric reaction 
capable of yielding the replenishment necessary to maintain the observed 
abundance, it is worth examining whether the nitrous oxide could originate on, 
or near, the surface of the Earth. 


Taste III 
Photo-dissociation Rate of Nitrous Oxide 
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(/N,O molecule/sec) 
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The gas is known to be produced by a variety of widely distributed micro- 
organisms (36), and Adel (64) advocates the view that these are responsible tor 
its presence in the atmosphere. Proper quantitative data to test the hypothesis 
are lacking. Bearing in mind other loss processes {such as the leaking of nitrates 
into the subsoil or into the drainage water (65§)}, it is reasonable to assume however 
that the rate of formation of nitrous oxide by soil micro-organisms is less than the 
rate of fixation of nitrogen over the land which, according to Hutchinson (66), 
averages only some 2 x 10!” to 1 x 10'* N, molecules/cm?/year.* ‘The reliability 
of this figure is difficult to assess.t As would be expected, it is much less than the 
rate of assimilation of carbon (cf. Section 4): but the difference is not so great as 
to suggest that it is a gross underestimate. A yield of perhaps 5 x 10! N,O 
molecules/cm?/year would thus seem the maximum possible ; and a much smaller 
yield would seem more probable. If this is correct, then 40 years or more would 
be required to produce the amount of nitrous oxide indicated by the spectroscopic 
observations. It will be noted that 5 x 10!” N,O molecules/cm?/year is equivalent 
to a production rate of 3 x 10 N,O molecules/cm*/sec distributed over an altitude 
range of 6km, which is of the same order as the calculated destruction rate in the 
troposphere (cf. Table III). Thus if the upper limit to the estimated possible 
yield from them is approached, micro-organisms might well be the source of the 
nitrous oxide. Research is required to determine if indeed the gas is not merely 
an important intermediary, but is a major end-product of denitrification (as is 
necessary). It would be surprising if this should prove to be the case. 

* It may be remarked that lightning discharges give only a minor contribution to the fixation 
rate and moreover lead mainly to compounds other than nitrous oxide. As may easily be verified, 
they are therefore unimportant in the present connection. 


t In a discussion of the nitrogen cycle in agriculture Russell (65) quotes rates which are more 
than an order of magnitude higher; but how representative these are he does not indicate, 
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Nitrous oxide seems to occur in the air near ground level. ‘Ihe mass spectro- 
scopic measurements of Slobad and Krogh (§5§) indicate that its fractional volume 
concentration there is about 5 x 10°’, which is not significantly different from the 
mean (1 x 10°) deduced from the total amount in a vertical column. Con- 
firmatory evidence for its presence was obtained by Gebbie ef al. (67) during 
the course of some work on the atmospheric transmission in the infra-red. It has 
been argued (64) that this strongly favours the suggestion that the nitrous oxide 
originates in the soil. But as will be seen, homogeneous gas reactions may also 
provide a source at low levels. 

The various possibilities can easily be listed In examining them it must be 
remembered that in general the parent particles concerned must themselves be 
appearing at a rate at least equal to the rate at which they are reacting; and 
this in turn must be at least equal to the rate of photo-dissociation of the nitrous 
oxide. Many of the mechanisms which have been discussed by earlier workers 
fail to satisfy the requirement. 

Bamford (10) has suggested that N.O may be formed from O* ions in some 
unspecified way, Nicolet (68) that it may be formed from O- ions through 
associative detachment 

O-+N,>N,O+e (46) 
(though this is endothermic). Both suggestions are unacceptable, since the rate 
of production of neither ion {cf. Bates and Massey (69)} even approaches the 
necessary value. 

Amongst Bamford’s other proposals is that metastable nitrogen molecules 
N,(A%Z,;) may be involved. ‘The reactions envisaged are presumably 

N(A®L})+0+M+>N,0+4+M, (47) 

N,(A3E,') + O, > N,0+0., (48) 

N,(A®Z})+ 0, > N,0+0. (49) 
These are not at all attractive. Under upper atmospheric conditions the 
formation of N,(A*=,)) is difficult. The rate of direct excitation by solar radiation 
is utterly negligible. In connection with the theory of the airglow much abortive 
effort has been directed to finding a plausible collision process. ‘The only ones 
which provide sufficient energy appear to be the three-body association of nitrogen 
atoms through 

N+N+M-N,(A*Z))+M (50) 
or 


N+N+N,>N,+4+N,(A 22). (51) 


Even if the rate coefficients of these are high (which is unlikely in view of the 


electronic transition incurred). there remains the necessity of providing an adequate 


source of nitrogen atoms. Photo-dissociation of molecular nitrogen seems to 
be the sole possibility. Laboratory experiments prove that this does not occur 
readily {cf. (70)}; and reference to the figures given above shows that it should 
actually be taking place in the upper atmosphere at a rate which is comparable 
with the rate of photo-dissociation of molecular oxygen. ‘There can therefore be 
little doubt * that the supply of atomic nitrogen (and so of metastable nitrogen 
molecules) is quite insufficient. Moreover the attenuation of the incoming 

* The fact that some of the nitrogen atoms may be regenerated through (42) does not appreciably 


weaken the argument, since only a small fraction of them would follow the reaction path leading to 
nitrous oxide, 
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radiation is such that any supply there is would be well above the probable location 
of the nitrous oxide layer. It may also be urged that (48) is unimportant in this 
region. Again, the concentration, m(N,A*X,'), is unlikely to exceed 1 x 10°/cm*® 
(cf. Appendix). On taking * n(O) to be 1 x 10!3/cm* and n(M) to be 3 x 10!/cm? 
at the 100 km level, it can be seen that in order that (47) should yield the required 
4x10’ N,O molecules/cm*/sec its rate coefficient would have to be some 
10° *4cm®/sec, which is of course impossible. In the case of (49) this alternative 
approach is also useful. For since m(O,) does not exceed 5 x 10!*/cm® a rate 
coefficient of at least 10°!’ cm*/sec would be needed. Similarly 


N+N+0,>N,0+0 (52) 


is most unpromising, as the rate coefficient would have to be more than 10-*! cm®/sec 
even if m(N) were 1 x 10!2/cm*. ‘Though these last two rate coefficients do not 
imply collision cross-sections exceeding gas kinetic values, they are improbably 
large. It is to be noted in this connection that the formation of nitrous oxide from 
active nitrogen and molecular oxygen does not appear to have been reported 
{ef. (73)}. 
No process involving metastable oxygen molecules, O,(a'A, or b12}), is at 
all plausible. For example, 
O,(a'A, or 6'X7)+N,>N,0+0 (53) 


is strongly endothermic, so that it can be dismissed without detailed consideration. 
The low heat of formation of N,O also prevents most oxides from interacting 
with molecular nitrogen. ‘To be sure 


HO, + N,>N,O0+0OH (54) 


may be almost thermochemically neutral {cf. Bodenstein, Schenk and Kornfeld 
(72) and Gaydon (23)}. However, its rate is probably slow, since the influence 
of molecular nitrogen on the hydrogen—oxygen reaction seems to be merely that 
of an inert gas {Gibson and Hinshelwood (73)}. Furthermore it may be shown 
from the work of Bates and Nicolet (3) that the perhydroxyl radicals in the upper 
atmosphere are not replenished sufficiently rapidly. 

The main processes thus far discussed (that is, processes involving positive 
or negative ions, metastable molecules or oxides) all fail, in that the supply of 
the parent particles is inadequate to balance the loss that would be incurred 
in producing the nitrous oxide at as great a rate as it is being destroyed. This rate 
is indeed so great that it was clear from the onset that the only mechanism likely 
to give sufficient parent particles is the photo-dissociation of molecular oxygen. 
The reactions following such photo-dissociation would have been considered 
earlier, but that it was judged important to demonstrate the need for invoking 
them as conclusively as possible by first disposing of the alternatives that had 
been proposed. 

The altitude distribution of the rate of photo-dissociation of molecular 
oxygen was obtained in the course of the calculations outlined earlier. It is 
characterized by two maxima, one due to the Schumann—Runge continuum and 
another due to the Herzberg continuum, located at altitudes of approximately 
100km and 30km respectively. The yield of oxygen atoms in the first is some 
I x 107/cm*/sec; that in the second is some 5 x 107/cm*/sec. When the particular 
reactions which might lead to the formation of nitrous oxide are described it 


* Rather high values of the concentrations are deliberately chosen. 
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will be apparent that they could not occur to a sufficient extent at the altitude of 
the Schumann-Runge maximum. For this and other reasons attention will 
be confined to the Herzberg maximum. 

In the region near 30km the atomic oxygen initially formed reacts rapidly 
to give ozone, from which it is rapidly re-liberated by further photo-dissociation 
(cf. Section 2). As a result the supply of either substance is sufficient to meet 
the nitrous oxide requirement. With even normal oxygen atoms and ozone 
molecules there is the possibility of 

O+N,+M->N,0+M (55) 
and 
O,+ N,>N,0+ O,. (56) 
‘hese processes certainly do not occur readily {cf. Bamford (10)}, but as the 
concentrations of the reactants are considerable it does not follow that they are 
unimportant. In addition both types of particle will make collisions while still 
retaining some of their energy of formation. ‘The O('D) atoms* provided by (6), 
and the activated O,* molecules provided by (4), may yield N,O by reactions which 
take the same course as (55) and (56) and will be referred to as (55’) and (50’) 
respectively. 

Just above ground level the temperature is approximately 285 deg. K, n(N,) 
is 2x 10!/cem$; n(O,) is 5 x 10'*/em*; n(O,) is perhaps some 5 x 10"/cm?* and 
n(O) is perhaps of order 1 x 10/cm*: further O,* molecules and O('D) atoms 
are being formed at rates of about 2 x 105/cm*®/sec and 4 x 10’,cm*;sec respectively. 
Thus there would be a yield of say 5 x 10 N,O molecules/cm* sect provided 
either 

(i) the rate coefficient of (55) is I x 10 **cm®/sec; 

(ii) the rate coefficient of (56) is 5 x 10 *’ cm/sec; 

(iii) there is a I in I x 10% chance that an O('D) atom interacts through (55’) 
before deactivation (as would be the case if, for example, the rate coefficient 
of (55') were I x 10° *6cm®/sec and that associated with deactivation were 
2 x 10° !4cm*/sec); or 
there is a I in 4 x 10% chance that an O,* molecule interacts through (56’) 
before deactivation (as would be the case if, for example, the rate coefficient 
of (56’) were 6 x 10° !4cm*/sec and that associated with deactivation were 
1 x 10! cm?/sec). 

It is difficult to determine whether any one of these requirements is actually met. 
‘The following comments are however offered. 

If (55) leads to normal X!X N,O molecules the rate coefficient would scarcely 
be expected to be gveat as I x 10-** cm®/sec, since a reversal of electron spin would 
be involved and since the activation energy would probably be appreciable 
{cf. Stearn and Eyring (75)}. However, a final decision cannot yet be reached 
theoretically. Laboratory investigations, too, provide but inconclusive evidence 
on the matter. The most significant experiments are those carried out by 
Leighton and Steiner (76), who passed mixtures of molecular nitrogen and 
oxygen through a cell illuminated by radiation in the Schumann—Runge 
continuum so that oxygen atoms (half of them in the 'D state) were freed and 


* There does not appear to be an adequate source of O('S) atoms {cf. Bates and Seaton (74)}. 

t This is greater than the rate of photo-dissociation of molecular oxygen at ground level, but 
no difficulty arises, since the time scale is such that mixing with the air at higher levels ensures an 
ample supply of odd oxygen atoms and since moreover (41) tends to replace those lost. 
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reactions involving them could proceed. In one run the n(N,)/n(O,) ratio 
was as high as 10%; no nitrous oxide could be detected and apparently only ozone 
was formed. However, on referring to the detailed results it will be seen that it 
cannot be inferred with complete certainty that condition (i) is not satisfied. 

‘There appears no reason to suppose that the rate coefficient for (56) could 
not be 5 x 10°*7cm*/sec: for such a rate coefficient would be possible even with 
an activation energy of 20kcal; and, as it would allow an ozone molecule to have 
a life of several months in nitrogen at room temperature and atmospheric pressure, 
the process would not have been detected during any of the laboratory investiga- 
tions. 

Conditions (iii) and (iv) also do not appear to be excluded by the published 
experimental data. Neither is particularly severe. ‘Thus judging from the 
potential energy curve along which an O('D) atom and an N,(X!2>) molecule 
approach each other {cf. Stearn and Eyring (75)} (55’) should not be abnormally 
slow; since (56’) involves an O,* molecule its activation energy might well be 
small, as is implied by the suggested rate coefficient; and the assumptions 
regarding the deactivation of O('D) and O,* are not unreasonable. 

As the high activation energy makes the rate of (56) extremely sensitive to 
the temperature the yield of nitrous oxide from it must decrease through the 
troposphere. In view of the increase in the photo-dissociation rate (cf. ‘Table 111) 
it would clearly lead to a low-lying bank. Examination of the other processes 
shows that with them also the greatest equilibrium concentration would occur 
near ground level. However, the fall-off with altitude would be comparatively 
slow. Consequently (42) would give rise to a considerable number of odd 
nitrogen atoms. Now in the lower mesosphere their re-association is hindered 
by the rapid oxidation occurring {Bates (77)}. ‘here is therefore the possi- 
bility that (55), (55) or (56’) would lead to an unacceptable accumulation of 
nitric and higher oxides if they were as effective as required. With (56) no such 
difficulty arises. 


APPENDIX 
The concentration of N.(A*X,;) molecules in the upper atmosphere 


Illuminated metastable nitrogen molecules readily absorb and re-emit through 
the First Positive band system 


N,(A®E$) +Av>N,(B3I1,), 


N,(B°Il,) > N,(A23¢) + hy. (1) 


Resonance scattering at twilight should therefore provide a sensitive means of 
detecting their presence in the upper atmosphere. ‘To investigate the effect it 
is necessary first to determine the transition probability. This can be done by 
standard quantal methods. 

The oscillator strength is given by the formula 


f= (1/3Ray’A)D, (2) 


where R is Rydberg’s constant, a@» is the radius of the first Bohr orbit, A is the 
wave-length of the radiation concerned and D is 


| J X(o,2p |r)rX*(7,2p| ) dr|? (3) 
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with X(o,2p|r) and X(z,,2p|r) the normalized wave functions of the active electron 
in the A®Z and B§II, states respectively and r its position vector (ef. Mulliken 
(78)}. On the crude LCAO approximation 
X(o,2p) = No(Y4o" — Ho"), (4) 
X(7 2p) = Ny(i" +f"), (5) 
where N, and N, are normalizing factors, a and 6 denote the atom on which the 
orbital is based, and 
thy = f(r) cos 8, (6) 
ree gi cos 3 
dn =fosino {orf (7) 
Representing f(r) by the function 
(y°/7)"*r exp (—yr) (8) 
and substituting in (3) it can be shown that when the internuclear distance is 
Ra, then 
D=T?/(1—So)(1 + Sj), (9) 
where 


a ) 
T= 5 | 541—64,+ Ag, (10) 
5A, 


o5 ' 
S, = 2 {5Ay— 184, 4: p (11) 


S,=2T7/R, (12) 
a=yR, (13) 


A,= [exp —ax)x" dx, (14) 


Spectroscopic measurements show that R is approximately 24a, {cf. Herzberg 
(79)}. As usual « is poorly determined: a value of about 4 is indicated by the 
self-consistent field wave functions of Brown, Bartlett and Dunn (80) and a 
value of about 5 by Slater’s rules (which, though crude, have proved remarkably 
successful in molecular calculations). On carrying out the necessary computation 
it was found that D is either 0-068 or 0-026, depending on whether the smaller or 
larger « is used. Adopting the mean, 0-05, and taking the wave-length to be 
6000 A (v. infra), it can be seen from (2) that the oscillator strength is of order 
2x10-%, The corresponding Einstein coefficient A is 2 x 105/sec. 

If the effect of vibration in producing energy spread is neglected, then the 
photon emission due to (1) is 

P=sQA exp(—hc/ART)/N,(A*®={) molecule/sec. (15) 

where s is the ratio of the statistical weight of upper state to that of the lower (2), 
Q is the dilution (5-4 x 10~*) and T is the temperature (6000 deg. K) of the solar 
radiation, k is Boltzmann’s constant, and the other symbols are as already defined 
(81). Numerical substitution yields 


P=4x10-*/N,(A*X}) molecule/sec. (16) 
Montgomery and Nich@lls (82) have recently evaluated the nuclear overlap 


integrals associated with the individual bands. The intensity distribution 
through the system can thus be deduced (81). This need not be given in detail. 
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It is sufficient to mention that approximate calculations indicate that AA 6705, 
6624, 6545 and 6468 of the v’ — v” =3 sequence should be prominent, the photon 
emission in each being at least 0-05 times the total. Now resonance scattering of 
the First Positive system at twilight has never been observed.* If the photon 
emission in any of the bands mentioned were as great as 2 x 107/cm*/sec it could 
scarcely have escaped detection. From the estimated efficiency of the scattering 
process it is hence apparent that the concentration of metastable nitrogen mole- 
cules does not exceed 10°/cm*. ‘The assumed upper limit to the twilight intensity 
of the First Positive bands is considered to be conservative. A precise spectro- 
scopic determination of this limit would be very useful in that it would enable a 
more severe restriction to be placed on the concentration of metastable molecules. 


Dept. of Applied Mathematics, Dept. of Physics, 
Queen’s University, University College, 
Belfast: London. 
1951 October 2. 
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